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"Alguna vez, los senderos de ese laberinto convergen: por ejemplo, usted llega a esta 
casa, pero en uno de los pasados posibles usted es mi enemigo, en otro mi amigo. Si se 
resigna usted a mi pronunciacion incurable, leeremos unas pdginas. 7 |f| 

Anabella Araudo 



3 J.L. Borges, Extracto del cuento "El jardin de los senderos que se bifurcan", Ficciones. 



Indice general 



Resumen 
Abstract 

1. Introduccion general 

I Procesos no termicos 

2. Ondas de choque y aceleracion de particulas 

2.1. Introduccion 

2.2. Definition y estructura de las ondas de choque 

2.3. Ondas de choque no relativist as 

2.3.1. Ondas de choque hidrodinamicas 

2.3.2. Ondas de choque magnetohidrodinamicas 

2.4. Ondas de choque relativistas 

2.5. Production de ondas de choque 

2.5.1. Ondas de choque por perturbaciones locales . . . 

2.5.2. Ondas de choque por interacciones globales . . . . 

2.6. Aceleracion de particulas en choques 

2.6.1. Poblacion de particulas no termicas 

3. Procesos radiativos 

3.1. Conceptos basicos 

3.2. Procesos radiativos 

3.2.1. Radiation sincrotron 

3.2.2. Radiation Compton inversa 

3.2.3. Bremsstrahlung relativista 

3.2.4. Interacciones proton-proton 

3.3. Absorcion 

3.3.1. Absorcion por creation de pares electron-positron 



Estudios de fuentes a diferentes escalas 



Objetos estelares jovenes 

4.1. Introduction 

4.2. La fuente IRAS 16547-4247 

4.3. Poblacion de particulas relativistas 

4.3.1. Choques terminales 

4.3.2. Aceleracion de particulas y perdidas radiativas 

4.4. Distribuciones espectrales de energia 

4.4.1. Interacciones leptonicas 

4.4.2. Interacciones hadronicas 

4.5. Discusion 

Micro cuasares 

5.1. Introduccion 

5.2. Escenario 

5.2.1. Interaction jet- clump 

5.2.2. Emision termica del clump 

5.3. Poblacion de particulas relativistas 

5.3.1. Aceleracion de particulas y perdidas radiativas 

5.4. Distribuciones espectrales de energia 

5.4.1. Emision asociada al bow shock 

5.4.2. Emision asociada al clump 

5.5. Interacciones simultaneas 

5.6. Discusion 

Nucleos de galaxias activas 

6.1. Introduccion 

6.2. Escenario 

6.2.1. Interaction jet-nube 

6.3. Poblacion de particulas relativistas 

6.3.1. Aceleracion de particulas y perdidas radiativas 

6.4. Distribuciones espectrales de energia 

6.4.1. Emision asociada al bow shock 

6.4.2. Emision asociada a la nube 

6.5. Interacciones multiples 

6.6. Aplicaciones 

6.6.1. Galaxias FR I: Cen A 

6.6.2. Galaxias FR II: 3C 273 

6.7. Discusion 



7. Cumulos de galaxias Ill9 

7.1. Introduction 119 

7.2. El cumulo Abell 3376 121 

7.3. Poblacion de partfculas relativistas 125 

7.3.1. Aceleracion de partfculas y perdidas radiativas 126 

7.4. Distribuciones espectrales de energfa 131 

7.4.1. Interacciones leptonicas 132 

7.4.2. Interacciones hadronicas 132 

7.5. Discusion Il34 

8. Conclusiones y perspectivas Il35 
Referencias Il40 
A. Lista de acronimos 147 



B. Lista de publicaciones 



1491 



Resumen 



El objetivo principal de esta tesis es investigar los procesos fisicos que dan lugar a emi- 
sion no-termica a altas energfas en objetos astrofisicos capaces de acelerar particulas hasta 
velocidades relativistas. En particular, se ha estudiado la emision de rayos gamma produ- 
cida en fuentes cosmicas con diferentes escalas espaciales, desde objetos estelares jovenes 
hasta cumulos de galaxias, pasando por microcuasares y micleos de galaxias activas. En los 
dos primeros tipos de objetos se ha modelado la emision de rayos gamma a partir de los 
datos obtenidos en frecuencias radio de las fuentes IRAS 16547-4247 y Abell 3376. En los 
dos ultimos, se ha desarrollado un modelo especffico de emision basado en la interaction 
de inhomogeneidades del medio externo con los jets producidos por el objeto compacto. 
Especificamente, se han considerado clumps o grumos del viento de la estrella companera 
en los microcuasares y nubes de la region de formacion de lineas anchas en las galaxias 
activas, interactuando con los jets de las fuentes. En todos los casos, los modelos desarrolla- 
dos permiten realizar predicciones contrastables por la nueva generation de instrumentos 
que operan en altas energias, tales como los satelites Fermi y AGILE y los telescopios 
Cherenkov HESS, MAGIC y el planeado CTA. 



Palabras claves: 

■ Rayos gamma: general 

■ Ondas de choque: general 

■ Ondas de choque: aceleracion de particulas 

■ Procesos radiativos: no termicos 

■ Estrellas de gran masa: formacion: emision no termica 

■ Estrellas de gran masa: vientos 

■ Sistemas binarios: microcuasares: general 

■ Galaxias activas: general 

■ Cumulos de galaxias: Abell 3376 



Abstract 



The main goal of this thesis is to study the physical processes that can produce non- 
thermal emission at high energies in astrophysical objects capable to accelerate particles 
up to relativistic velocities. In particular, we have studied the gamma-ray emission produ- 
ced in cosmic sources with different spatial scales, from young stellar objects to clusters 
of galaxies, going through microquasars and active galactic nuclei. In the former cases, we 
have modeled the gamma-ray emission using the radio data from the sources IRAS 16547- 
4247 and Abell 3376. In the latter, we have developed a specific radiation model based on 
the interaction of the inhomogeneities of the external medium with the jets generated by 
the compact object. Specifically, we have considered clumps of the massive stellar wind 
in microquasars, and clouds of the broad line region in active galactic nuclei, interacting 
with the jets of the sources. In all cases, the developed models allow us to make predic- 
tions testables with the new generation of instruments operating at high energies, such as 
the satellites Fermi and AGILE, and the Cherenkov telescopes HESS, MAGIC, and the 
forthcoming CTA. 



Key words: 

■ Gamma-rays: general 

■ Shock waves: general 

■ Shock waves: particle acceleration 

■ Radiative processes: non-thermal 

■ Massive stars: formation: non-thermal emission 

■ Massive stars: winds 

■ Binary systems: microquasars: general 

■ Active galaxies: general 

■ Clusters of galaxies: Abell 3376 
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Capitulo 1 

Introduction general 



La astronomfa es la ciencia que estudia los fenomenos que ocurren fuera de nuestro 
planeta. La tecnologfa actual ha permitido estudiar in situ los planetas y satelites naturales 
cercanos, accediendo a ellos a traves de sondas y naves tripuladas. Sin embargo, para 
los procesos que ocurren fuera del sistema solar, la astronomfa aiin sigue investigando 
casi como en la antiguedad, es decir, a traves de la luz que llega a nuestros detectores 
desde los objetos celestes. Esta radiation, en cada banda de frecuencia, nos ofrece una 
fenomenologfa distinta, ya que los procesos ffsicos subyacentes pueden ser muy variados. 
En particular, la emision de rayos gamma da cuenta de los procesos no termicos, es decir, 
fuera del equilibrio termodinamico, que puedan tener lugar en la fuente. El Universo en 
rayos gamma es puramente no termico, ya que la temperatura requerida para emitir fotones 
gamma termicamente es extremadamente alta (~ 10 13 K) y diffcilmente puedan tenerla 
los sistemas fisicos conocidod^. 

Los rayos gamma forman la ultima banda del espectro electromagnetico, abarcando 
mas de 14 ordenes de magnitud en energfa: 



El lfmite inferior, E 1 ~ m e c 2 ~ 5 x 10 5 eV, corresponde a la emision de lfneas, como 
la de aniquilacion de los pares electron-positron e , mientras que el valor superior, E 1 ~ 
10 20 eV, corresponde a los rayos cosmicos mas energeticos que han sido detectados. Dada la 
gran amplitud de esta banda de energfa, resulta conveniente subdividirla en las siguientes 



■ E 1 < 30 MeV: Baja energfa (LE, por Low Energy). 

m 30 MeV < £ 7 < 30 GeV: Alta energfa (HE, por High Energy). 

■ 30 GeV < E y < 30 TeV: Muy alta energfa (VHE, por Very High Energy). 

1 Temper aturas ~ 10 13 K podrfan alcanzarse por periodos de tiempo cortos en eventos explosivos muy 
energeticos como el Big Bang o los eruptores de rayos gamma (GRBs, por Gamma Ray Bursts). 
2 1 TeV = 10 12 eV, 1 PeV = 10 15 eV y 1 EeV = 10 18 eV. 



5 x 10 5 eV < E~ < 10 20 eV. 
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- 30TeV < £ 7 < 30 PeV: Ultra alta energia (UHE, por Ultra High Energy). 

■ E 1 > 30 PeV: Extremadamente alta energia (EHE, por Extremely High Energy). 

La observation en las bandas LE y HE debe realizarse desde satelites espaciales, ya que 
la atmosfera terrestre es opaca para las frecuencias correspondientes. Por otro lado, los 
rayos gamma de energias > 30 GeV pueden detectarse indirectamente desde la superflcie 
de la Tierra mediante telescopios Cherenkov. Los rayos gamma, al penetrar en la atmosfera 
interactiian con los campos alii presentes (magneticos, de materia y de fotones) produciendo 
pares e ± . Estos leptones, interactuando con los mismos campos, pueden crear mas fotones 
y pares menos energeticos, desarrollando asi una cascada electromagnetica. La luz (visible) 
Cherenkov producida por los leptones relativistas en la atmosfera puede detectarse desde 
la superficie terrestre. Reconstruyendo la cascada se puede determinar la energia del foton 
gamma original y la direction de arribo, que nos indica la localization de la fuente. Sin 
embargo, solo hay detection de fuentes hasta la region de las UHE y por lo tanto solo para 
< PeV puede hablarse de una astronomia de rayos gamma. 

Los telescopios de rayos gamma que actualmente estan en funcionamiento son los si- 
guientes: 

■ Telescopios espaciales (E 1 < 30 GeV) 

• Fermi: Es el instrumento mas nuevo que se ha puesto en funcionamiento. Fue 
lanzado por la NASA en el aho 2008 y puede detectar fotones con energias 
20 MeV < £ 7 < 100 GeV. 

• AGILE: Es un satelite italiano que funciona en el rango de energias 0.1 GeV 
< E, < 30 GeV. 

• Swift: Lanzado por la NASA en el ano 2004, puede detectar fotones gamma con 
energias de ~ 150 keV. Ademas del instrumento BAT de rayos X duros (pero 
que llega a detectar algunos fotones gamma blandos), a bordo del satelite hay 
un telescopio ultravioleta (UV) y otro optico, con el fin de poder detectar GRBs 
en diferentes longitudes de onda. 

■ Telescopios Cherenkov {E 1 > 30 GeV) 

• HESS: Es un arreglo de 4 telescopios franco-alemanes de 12 m de diametro cada 
uno ubicados en Gamsberg, Namibia. Opera en el rango de energias 50 GeV 
<£ 7 <10TeV. 

• MAGIC: Es un unico receptor germano-espanol de 17 m de diametro ubicado 
en La Palma y detecta fotones gamma de energias 50 GeV < E 1 < 10 TeV. 

• VERITAS: Es un arreglo de 4 telescopios estadounidenses de 12 m de diametro 
cada uno, ubicados en Arizona, y que opera en la banda de energias 100 GeV 
<£ 7 <50TeV. 
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• Cangaroo III: Es una colaboracion japonesa-australiana que opera un arreglo 
de 4 telescopios de 10 m cada uno ubicados en Woomera, en el sur de Australia. 
Detecta rayos gamma de E 7 ~ 100 GeV. 

Todos estos instrumentos detectan fuentes, y de ahi que los llamemos telescopios. Estas 
fuentes pueden ser puntuales o extendidas, como por ejemplo la region central de la Galaxia. 
Entre las fuentes detectadas por los telescopios listados anteriormente, hay nucleos de 
galaxias activas (AGN, por Active Galactic Nuclei), pulsares, remanentes de supernovas, 
binarias de rayos X, GRBs y algunas regiones de formacion estelar. Hay otras fuentes que 
han sido observadas, pero cuya deteccion H no ha sido clamada aiin como ocurre con los 
cumulos de galaxias y algunas regiones de formacion estelar. Por otro lado, hay muchas 
fuentes que aiin no han sido identificadas. En la Figura 11.11 se muestra un mapa del cielo 
con las fuentes detectadas en el rango de las VHE por diferentes telescopios. 

La deteccion de fuentes gamma es importante porque nos provee information sobre 
procesos fisicos extremos. La astroffsica de rayos gamma es una disciplina unica para es- 
tudiar los fenomenos fuera del equilibrio que involucran las enormes cantidades de energfa 
requeridas para emitir los fotones mas energeticos del espectro electromagnetico. Asi, los 
emisores de rayos gamma deben ser regiones del espacio en las cuales las parti'culas rela- 
tivistas alh presentes pueden enfriarse eficientemente debido a interacciones con grandes 
densidades de materia y/o radiation. 

En esta tesis nos hemos propuesto estudiar los procesos no termicos que puedan desa- 
rrollarse en una variedad de fuentes astrofisicas con diferentes escalas espaciales, desde 
objetos estelares jovenes (YSOs, por Young Stellar Objects) hasta cumulos de galaxias, 
pasando por microcuasares (MQs, por Microquasars) y AGNs. En particular, nos hemos 
concentrado en estos cuatro tipos de objetos para analizar como las ondas de choque pue- 
den acelerar parti'culas cargadas hasta velocidades relativistas y como luego estas particulas 
pueden interactuar con el medio y producir fotones mediante procesos radiativos no termi- 
cos. Hemos calculado el espectro de fotones producido, haciendo enfasis en el analisis de la 
emision de rayos gamma. 

■ En los YSOs masivos, estudiamos la aceleracion de particulas en los choques termina- 
les de los jets, producidos cuando estos son frenados por el medio circundante. Hemos 
calculado el espectro de emision de rayos gamma de la fuente IRAS 16547-4247 obte- 
niendo niveles de luminosidad detectables con los telescopios que actualmente estan 
funcionando. 

■ En los microcuasares de alta masa (HMMQ, por High Mass Microquasar), el viento 
de la estrella companera puede tener estructura (inhomogeneidades en la densidad). 
Estas inhomogeneidades al llegar al jet pueden penetrar en el y producir ondas de 
choque. Las particulas alh aceleradas pueden producir emision variable y, en rayos 
gamma detectable por los telescopios actuales. Este mecanismo podria explicar cierta 
variabilidad observada en fuentes como Cygnus X-l, LS 5039 y LSI +31 603. 

3 En frecuencias gamma, para que una deteccion sea confiable se requiere una relation se nal (S) ruido 
(N) tal que S/N > 6. 
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VHE y-ray Sky Map 

(E t >100 GeV) 



+180°i 




VHE --ray sources 
# Blazar (HBL) 
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-180° 



2010-03-25 - Up-to-date plot available al http://www.mpprnu.mpg.de/~rwagner/sourGeB/ 



Figura 1.1: Mapa del cielo en rayos gamma de VHE. Se muestran las fuentes puntuales 
detectadas hasta la actualidad por diferentes telescopios. Entre las fuentes identificadas, las 
mas comunes son los AGN, en particular, los blazares (http:/ /www. mppmu.mpg.de/ rwag- 
ner /sources/). 



■ En los AGNs estudiamos una situation similar a la analizada en los HMMQs. Sin 
embargo, en los primeros, seria la interaction de las nubes que se encuentran orbitando 
al agujero negro central con el jet la que daria como resultado la emision de rayos 
gamma. Este mecanismo podria explicar la emision observada en algunas galaxias. 

■ Finalmente, en el cumulo de galaxias Abell 3376 estudiamos el espectro en rayos 
gamma producido por la poblacion de particulas relativistas generada en choques 
detectados en el borde de la fuente. 

Los resultados de nuestra investigation sobre "Radiation no termica asociada a las 
ondas de choque astrofisicas" , tematica que da nombre a esta tesis, seran presentados 
de la siguiente manera: en los Capitulos 2 y 3 haremos una somera description de los 
procesos no termicos: ondas de choque y radiation. Los cuatro capitulos siguientes estan 
dedicados a cada una de las fuentes estudiadas en esta tesis: YSOs (Capitulo 4), HMMQs 
(Capitulo 5), AGNs (Capitulo 6) y ciimulos de galaxias (Capitulo 7). Finalmente, en el 
Capitulo 8, discutiremos los resultados y las conclusiones. A lo largo de toda la tesis, los 
valores de las magnitudes fisicas consideradas estaran dados en unidades cgs y las unidades 
electromagneticas seran Gauss (G) y esu. 



Parte I 
Procesos no termicos 



Capitulo 2 

Ondas de choque y aceleracion de 
particulas 

2.1. Introduccion 

En un plasma que se mueve a una velocidad v, tiene un campo magnetico B, una presion 
P y una densidad p, en el cual pueda aplicarse la aproximacion magnetohidrodinamica0 
(MHD), las ecuaciones que describen el estado del sistema pueden separarse en 3 grupos: 
mecanicas, electromagneticas y termodinamicas. 

■ Las mecanicas incluyen la conservation de la masa 

^ + pV-t; = 0, (2.1) 

y la ecuacion de movimiento 

= -VP + J x B + F', (2.2) 

donde F' incluye a todas las demas fuerzas, ademas del gradiente de P y de la fuerza 
de Lorentz (I x B), que actuan sobre el sistema. La fuerza F" podria ser la gravedad, 
la fuerza viscosa, etc. Por otro lado, la derivada convectiva D/Dt se define como 
D/Dt = d/dt + v ■ V. En la aproximacion MHD, la relation entre I y B esta dada 
por la ecuacion de Ampere-Maxwell, (Air/c 2 ) I = V X B. 

■ Dentro de las ecuaciones electromagneticas incluimos a la de divergencia nula: 

V-5 = 0, (2.3) 



1 La aproximacion magnetohidrodinamica considera que el campo magnetico es relevante en la dinamica 
del plasma. En las ecuaciones de Maxwell, esto se introduce despreciando la corriente de desplazamiento. 
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Ondas de choque y aceleracion de particulas 



y a la ecuacion de induction 

dB - c 2 

donde o"b es la conductividad magnetica. 

Finalmente las ecuaciones termodinamicas son la ecuacion de estado, por ejemplo 
la de gas ideal, 

P = nK B T, (2.5) 
y la ecuacion de conservation de la energia: 



p 7ad d / P 



7 ad - 1 Dt VP 7ad 



-A (2.6) 



donde 7 a( j es el indice adiabatico y toma el valor 5/3 para gases monoatomicos no 
relativistas y 4/3 para los relativistas. Las perdidas de energia C pueden estar debidas 
a la viscosidad, a la difusion termica o a la radiation. En la ecuacion ( 12.5f) n = p/m es 
la densidad numerica de particulas de masa m, T es la temperatura y Kb la constante 
de Boltzmann. 

Si provocamos una pequena perturbation en el movimiento del plasma, las propiedades 
del sistema se veran modificadas de la siguiente manera: 

£ = 6) + a(r- t), (2.7) 

donde £i(r, t) es la perturbation respecto del estado de equilibrio £ - La magnitud £ puede 
ser la presion, la densidad, la velocidad o el campo magnetico. Si la perturbation es chica, 
es decir, |^| -C |£o|; entonces las ecuaciones de la MHD perturbadas se pueden linealizar 
considerando solo los terminos de primer orden en las perturbaciones £i. Si no hay campos 
magneticos, la perturbation satisface la ecuacion de onda: 

V>6 - ^ = 0, (2.8) 

donde C es (el modulo de) la velocidad de propagation de la perturbation^. 

La ecuacion de movimiento ( 12. 21) depende de las fuerzas que actuen sobre el sistema: 
el gradiente de presion, la fuerza de Lorentz u otro tipos de fuerzas F" . Dependiendo de la 
importancia relativa de unas respecto de las otras, tendremos diferentes tipos de ondas: 

■ Si | VP | ^> | J A B\ y | VP | ^> |P'|, y consideramos un fluido inicialmente en reposo, 
entonces tendremos variaciones en P, p y v. Estas son las llamadas ondas de sonido 
y se propagan a la velocidad del sonido C s = a/ 7ad P/ P en el medio no perturbado. 



2 A los modulos de las magnitudes vectoriales los representamos a traves de la letra usada para re- 
presentar a la magnitud pero sin la flecha que indica la naturaleza vectorial de la misma. Por ejemplo, 
\A\ = A. 



2.2 Definicion y estructura de las ondas de choque 
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■ Si | J A B\ | VP | y \I A B\ 3> |P'|, entonces tendremos las llamadas ondas magne- 
tohidrodinamicas o de Alfven. Estas tienen diferentes modos dependiendo de que la 
variation sea en la direccion o en la intensidad de B. Las ondas de Alfven de corte 
son transversales (se propagan perpendicularmente a Bq) y varian solo la direccion 
de B. La velocidad de fase de estas ondas es Vf = ±va cos(6>b) k/k, donde k es el 
vector de onda y 9 B es el angulo entre k y B , y la velocidad de grupo es v g = va- La 
velocidad de Alfven se define de la siguiente manera: va = Bq/ p$. Por otro lado, 
las ondas de Alfven de compresion son longitudinales; modifican tanto la direccion 
como la intensidad de B y se propagan con una velocidad Vf = v g = ±va k/k. 

■ Si \I AB\ ~ | VP | y \F'\ ~ 0, entonces tenemos las llamadas ondas magnetoacusticas. 
La velocidad C s a de propagacion (de fase) de estas ondas es < C s a < a/^a + C s 2 . 
Cuando max{v a, C s } < C s a < \Aa + se dice que la onda esta en el modo rapido, 
mientras que si < C s a < min{t>A, C s } se dice que esta en el modo lento. 

Notemos que en todos estos casos, la velocidad de fase, V{, y por lo tanto la de grupo, 
son constantes. Esto implica que no hay dispersion, es decir, que el perfil de las ondas 
no cambia, contrariamente a lo que ocurre, como veremos luego, en las ondas de choque. 
Hasta aqui hemos realizado un breve resumen de lo que sucede cuando las perturbaciones 
son pequenas. /,Que ocurre cuando las perturbaciones ya no son tan chicas, comenzando a 
hacerse importantes los terminos no lineales? 

2.2. Definicion y estructura de las ondas de choque 

Dada una onda sonora o una onda magnetoacustica (de compresion), cuando la ampli- 
tud es grande, es decir cuando los terminos no lineales de las ecuaciones son importantes, el 
perfil de la onda se modifica con el tiempo. Esto se debe a que las perturbaciones de mayor 
amplitud viajan mas rapidamente que las de menor amplitud, es decir v g no es constante, 
y el perfil de la onda se modifica empinandose cada vez mas a medida que la onda avanza, 
como se muestra en la Figura (Priest 1982, Platzeck 2009). 

Este comportamiento es similar tanto en v como en P y en p. Los gradientes de estas 
cantidades se hacen cada vez mas grandes de tal manera que efectos disipativos como la 
viscosidad o el termino disipativo de la ecuacion de induction no pueden despreciarse. Son 
estos efectos los que impiden que esta situation continue indefinidamente. Cuando la fuerza 
viscosa (no consideraremos disipacion magnetica) se equilibra con los gradientes de v, P y 
p se alcanza el estado estacionario y se frena el empinamiento del perfil de la onda. Luego 
continua propagandose con este perfil constante y empinado. Esto es lo que llamamos 
una onda de choque. Cuando se llega a esta situation, la velocidad de propagacion, que 
llamaremos t> c h, es mayor que la velocidad de la onda linealizada correspondiente, C s o C s a- 

El equilibrio entre la fuerza viscosa y los gradientes ocurre en una region muy peque na 
comparada con el tamano del sistema. El ancho Ax de esta region de transition se estima 
que es del orden de algunos caminos libres medios A (longitud que recorre una particula 
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Figura 2.1: Empinamiento del frente de la onda a medida que los terminos no lineales de 
las ecuaciones perturbadas se hacen mas importantes. 



antes de desviarse significativamentqfl) y decrece a medida que crece. De esta manera, 
si v ch ^> C s el choque se dice fuerte y la region de transition puede considerarse como 
una "superficie" de discontinuidad. El tiempo que tarda el fluido en atravesar la region 
de transition es muy chico comparado con el tiempo de los cambios en las regiones fuera 
de ella. De esta manera podemos considerar que el material que se encuentra fuera de la 
region de transition esta en un estado estacionario, como se muestra en la Figura [2T2l En la 
aproximacion hidrodinamica (HD) y MHD Ax ~ 0. Podemos definir a las ondas de choque 
de la siguiente manera: 

■ Onda de choque: Onda de gran amplitud que produce una discontinuidad en los 
valores de las funciones que representan a las propiedades del plasma. El frente de 
choque divide al medio en dos regiones: chocada y no chocada, en las cuales los valores 
de P, p, T, v y B son muy distintos. 

Denominamos con el submdice a las propiedades en el medio sin perturbar (delante 
del frente de onda) y con el submdice 1 a aquellas en el medio chocado (detras del frente 
de onda). Suponemos que ambos medios son uniformes. En el sistema de referenda (SR) 
con el medio en reposo, el choque se mueve con una velocidad U y el medio perturbado 
con una velocidad U\. En el SR fijo al frente de choque, el medio no perturbado se acerca 
al choque con una velocidad v mientras que el material perturbado se aleja del choque 
con una velocidad v%. En la Figura [2731 se esquematizan ambas situaciones. Notar que 

v = -U y 1 = U 1 -0 O . (2.9) 

Los modulos de las velocidades se relacionan de la siguiente manera: Vq = U y v\ = 
v -U 1 = U-U 1 . 

3 Si el mecanismo de interaction dominante entre las particulas es colisional, entonces A es la distancia 
recorrida por la particula antes de interactuar con otra. En el proximo capitulo (ver ecuacion (|3.6jl ) 
definiremos A de una manera mas apropiada. 
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► Onda de choque 

Ax 



Figura 2.2: Esquema de las variaciones en la presion (P), la densidad (p), y la velocidad (v) 
del fluido. Con los submdices y 1 indicamos las magnitudes en la region no perturbada y 
perturbada, respectivamente. En la figura de la izquierda se muestra la region de transicion 
Ax ampliada con respecto a la figura de la derecha. En HD y MHD, Ax — > 0, es decir, la 
region de transicion puede considerarse como superficie de discontinuidad de las magnitudes 
graficadas. 




Figura 2.3: Velocidades en el SR con el medio no perturbado en reposo (izquierda) y en el 
SR solidario al choque (derecha). 
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2.3. Ondas de choque no relativistas 

Las ondas de choque no relativistas (NR) son aquellas tales que t> c h < 0.1c, siendo c 
la velocidad de la luz en el vatio. Veremos a continuation que en el caso de choques no 
relativistas los valores de la velocidad, la presion, la densidad y el campo magnetico en 
los medios chocado y no chocado pueden relacionarse teniendo en cuenta las ecuaciones 

2.3.1. Ondas de choque hidrodinamicas 

Las ondas de choque hidrodinamicas son aquellas que se producen en un medio en el 
cual la fuerza dinamicamente relevante es VP. Esto es |VP| ^> \I AB\ y |VP| ^> \F'\. Asi, 
las ondas de choque hidrodinamicas se propagan a una velocidad mayor que las ondas de 
sonido, es decir, t> ch > C s . 

Consideremos un medio no conductor y en estado estacionario, con densidad po, presion 
Po, campo magnetico B ~ y un choque propagandose con una velocidad v^. En el 
SR solidario al frente de choque, la velocidad del medio no chocado es v$ y debido a que 
hemos tornado convenientemente el eje x como la direccion de propagation del choque, solo 
habra variaciones espaciales en tal direccion. En este caso de ondas de choque propagandose 
en un medio con B ~ 0, las ecuaciones que describen el estado del sistema son la ecuacion 
( 12. ip de conservation de la masa, la ecuacion ( 12. 2 j) de movimiento y la de conservation de 
la energfa ( 12. 6p . 

Debido a que el sistema esta en estado estacionario, la ecuacion (12. ip se simplifica a 

d(pv) = Q 
dx 

Esta condition debe cumplirse a lo largo de cualquier lmea de fluido, por lo cual, si la 
integramos entre un punto en la region 1 (xi) y otro en la region no perturbada (xq) 
obtenemos: 

Piv 1 = p v . (2.11) 

Esta ecuacion describe que el flujo de masa pv que pasa a traves de la onda de choque se 
conserva. Por otro lado, la ecuacion de movimiento en un medio no conductor (J A B = 0) 
y estacionario se reduce a 

_ dv dp 4 d A dv 
^ dx dx 3 da; \ v dx 

La fuerza viscosa F' = —(4/3) d(A^dv/dx)/da; i, donde X v es una constante, esta presente en 
la region de transition Ax pero es nula tanto en el medio chocado como en el no chocado. 
Por lo tanto, si integramos la ecuacion (I2.12p entre Xi y x y usamos que piv 1 = poV , 
obtenemos 

p v 2 + P = p 1 v 2 1 +P 1 . (2.13) 



(2.10) 



(2.12) 
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La magnitud (pv) ■ v, conocida como "presion cinetica", es la cantidad de impulso trans- 
portado por unidad de tiempo y por unidad de area transversal al movimiento. Por otro 
lado, P es la fuerza que actua sobre dicha area. 

Finalmente, en la ecuacion de conservacion de la energfa consider aremos las perdidas 
por conduccion termica C c = V(— A C VT), donde A c es una constante, y por viscosidad 
£ v = (— 2/3)pA v (V • v) 2 , ambas presentes solo en la region de transition. La ecuacion (12. 6 j) 
escrita en funcion de la energfa interna por unidad de masa e y teniendo en cuenta la 
viscosidad y la conduccion de calor toma la forma 

fde_ p_dp\ = _d_ A dT\ + 2 /dw\ 2 
\ da; p 2 dx J dx \ ° dx / 3 v \ dx / 

Integrando como antes entre xi y xo, donde dT/dx = 0yd-y/dx = 0,y usando las relaciones 
f l2~TT|) y fl2TT3|) . obtenemos: 

^pit> 2 + eipij iTi + PiUt = Qpo^o + e oPo^ + PoVq, (2.15) 

donde ((l/2)pt> 2 + ep) es la cantidad densidad de energfa cinetica e interna transportada 
por unidad de superficie y de tiempo; pv es el trabajo por unidad de tiempo y de area 
realizado por la presion del gas. Si consideramos que piVi = PoVq tenemos 

Pi v 2 Po v 2 , 

Pi 2 Po 2 

Teniendo en cuenta que para un gas ideal e = p/p( 7 — 1), la ecuacion de conservacion de 
la energfa nos da la relation 

7ad Pi + Vf = 7ad Pp + Vq ^ ^ 



7 ad - 1 pi 2 7ad - 1 Po 2 
Asf hemos obtenido las relaciones de Rankine-Hugoniot, o relaciones de salto: 

p v = Pl v t (2.18) 
povl + Po = Piv\ + Pi (2.19) 

7ad^0 1 2 7ad^l 1 2 



+ _^ = — — ^ + (2.20) 

( 7ad - l)p 2 ( 7ad - l) Pl 2 1 y ' 

Resolviendo estas ecuaciones podemos obtener information sobre como la onda de choque 
modifica el medio. De la ecuacion (12.1 8[) obtenemos que 

* = *. (2.21) 

vo Pi 

Por otro lado, de las ecuaciones (12. 19ft y (12.201) podemos hallar expresiones para po/pi y 
para Pq/P\ en funcion del mimero de Mach M en el medio no perturbado: 
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Se obtienen las siquientes relaciones (Landau & Lifshitz, 1959): 

pi = vo _ (7ad + 1)M 2 
Po vx ( 7ad - 1)M 2 + 2 

Pi = 2 7ad M 2 - ( 7ad - 1) 

Po 7ad + 1 

Si ademas consideramos la ecuacion de estado (12.51) . entonces la temperatura del medio 
chocado es T\ = P\f [K^ri\) y en funcion de parametros conocidos el cociente Ti/T resulta: 

Ti [2 7ad M 2 - ( 7ad - 1)] [( 7ad - 1)M 2 + 2] 



(2.23) 
(2.24) 



T ( 7ad + 1) 2 M 2 



(2.25) 



Recordemos que todas estas expresiones (I2.23l) - fl2.25l) han sido deducidas en el SR que 
se mueve con el frente de choque. Teniendo en cuenta que v\ < vo, podemos decir: 

■ La onda de choque comprime el gas: pi > po- 

■ La velocidad de la onda de choque es supersonica: U = vq > C So . 

■ El material chocado es subsonico: Mi < 1. 

■ La onda de choque aumenta la presion del medio: Pi > P . 

■ La entropia aumenta: S\ > So. 

■ La onda de choque calienta el plasma: T\ > To. 

Si bien la presion crece si la velocidad del choque aumenta, el incremento de la densidad 
esta acotado: 

hm a = 2=i±i. (2.26) 

Mo-kx) p 7ad - 1 

Para el caso particular de un gas monoatomico, 7ad = 5/3 y se obtiene pi/po = 4. Cuando 
M ^ 1 se dice que el choque es fuerte y en este caso las relaciones de salto se expresan 
de una manera muy sencilla: 

Pl = 4p y u i = j> ( 2 - 27 ) 

Pi = \povl y Ti~2x 10-% 2 K, (2.28) 

y estas relaciones valen en toda la region chocada. 

Los choques que hemos estudiado en esta seccion se conocen como choques adiabati- 
cos, ya que hemos considerado que no hay perdidas de energia ni en la region chocada ni 
en la no chocada, es decir, C\ = Co = 0. Sin embargo, hay situaciones en las cuales las 
condiciones en la region 1 (originalemente adiabatica) son tales que las perdidas radiativas 
comienzan a hacerse relevantes. En este caso, la propagacion del choque se modifica dismi- 
nuyendo vecv c h y el estado del medio 1 cambia, haciendose mucho mas denso y frfo, como 
veremos a continuacion. 
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Choques radiativos 

Dado un choque que se propaga a una velocidad v^, y cuyas perturbaciones en el 
medio estan descriptas a traves de las ecuaciones (I2.23p . (I2.24p y (I2.25p . estas ultimas 
pueden no valer en toda la region chocada si las perdidas radiativas en dicho medio son 
importantes. Si bien en la region mas cercana al choque el estado del medio es aquel 
descripto por las condiciones de un choque adiabatico, a mayores distancias la densidad 
crece y la temperatura disminuye. La region 1 es entonces muy densa y fria. 

Para estimar la importancia de las perdidas radiativas en el medio chocado debemos 
comparar el tiempo de enfriamiento iter con otros tiempos caracterfsticos de cada sistema. 
Si consideramos las perdidas por emision termica (de lmeas y continuo) £ter ~ n 2 A(T), 
donde A(T) = A t T at (Myasnikov et al. 1998): 

f 7 x 10- 19 T-°- 6 si 10 4 < T < 4 x 10 7 K 
A[I } ~ \ 3 x 1(T 27 T - 5 si T > 4 x 10 7 K, [Z - ZJ) 

entonces t teT = (5/2)P 1 /£ ter . Por otro lado, si consideramos como tiempo caracterfstico 
aquel en el cual el material chocado recorre una distacia caracterfstica X, entonces ~ 
X/Ui y el choque sera radiativo si t tCT < t d . 

En este caso, debido a las perdidas radiativas, la energfa no se conserva a lo largo de 
una hnea de flui'do en la region chocada. La ecuacion ( 12 .6p se escribe ahora como (Zhekov 
& Palla, 2007) 

d < ^ . PlVl + Elf) = (2 .30) 



dx V7 ad - 1 2 
y si consideramos que P\ se mantiene constante, entonces la ecuacion anterior puede rees- 
cribirse como 

-n 1 v 1 K B - £ tcr = —n\ A(T X ). (2.31) 

2 ax 

Resolviendo esta ecuacion obtenemos el siguiente perfil de temperatura en la region cho- 
cada: 

/ O p2 „ \ V(3-at) 

Ux) = — l —zA(3 -a t )x+ (T? d ) 3 - at , (2.32) 

donde P x y Tf son la presion y la temperatura obtenidas con las relaciones de Rankine- 
Hugoniot. En la Figura 12.41 se muestra el perfil de temperatura en la region chocada de 
un choque radiativo obtenido bajo la aproximacion de presion constante en dicha region. 
Conociendo Ti(sc), la densidad puede hallarse a partir de la ecuacion de estado. Para fluidos 
ideales la densidad en la region 1 varfa de la siguiente manera: 

ni (x) = ^ (2.33) 

De las ecuaciones (I2.32p y (I2.33H vemos que a medida que el material se aleja del frente 
de choque, 7\ disminuye drasticamente (ver la Figura 12. 4p y n\ crece. De esta manera, el 
medio chocado se enfrfa en un tiempo menor que el tiempo en el cual el choque recorre una 
distancia del orden del tamano del sistema. La region adiabatica del medio 1 es pequena y 
luego el material se enfria hasta que se llega a un estado de alta densidad; muy compacto. 
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Figura 2.4: Perfil de temperatura en el medio chocado. Se muestran las soluciones exactas 
(lmea punteada) y las calculadas con la aproximacion de presion constante (lmea de rayas). 
Los niimeros 1, 2, 3 y 4 indican que las velocidades de los choques consideradas han sido 
v ch = 500, 600, 750 y 900 km s" 1 (Zhekov & Palla, 2007). 
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2.3.2. Ondas de choque magnetohidrodinamicas 

En el caso en que B ^ y relevante dinamicamente, entonces tendremos ondas de 
choque MHD. Las hay de diferentes tipos dependiendo de la geometria de B. 



Choques transversales 



Los choques transversales son aquellos en los cuales el campo magnetico en la region 
no perturbada, B , es perpendicular a la velocidad de propagation del choque. Es decir, 
Bo _L v^. En este caso, ademas de considerar las ecuaciones de movimiento, de conser- 
vation de la masa y de conservation de la energia, se considera tambien la ecuacion de 
induction. Asi, aparece en las relaciones de Rankine-Hugoniot ( I2.19P y (12.201) el termino 
B 2 /87i correspondiente a la presion y a la densidad de energia magnetica. Se tiene la 
relacion adicional 

t = t- (2 - 34) 

B Po 



Definiendo el parametro ( = pi/po, se tiene que ( satisface una ecuacion cubica. La 
solucion ( = 1 corresponde a una onda de Alfven de corte si v q = v\ , que no son ondas 
de choque. (Las ondas de Alfven de corte son ondas MHD que no se deforman aunque la 
amplitud de la pertubacion sea grande.) La otra solucion es ( < 0, pero se descarta por 
no tener sentido fisico. Finalmente, la solucion ( > corresponde a una onda de choque. 
En esta situacion se tiene un choque propagandose perpendicularmente a B y el efecto de 
este ultimo es reducir el valor de ( respecto del caso hidrodinamico. Es decir, en presencia 
de un campo magnetico, el aumento de la densidad de la region chocada es menor que en 
el caso hidrodinamico. 



Choques oblicuos 

Este caso corresponde a aquel en el cual B es oblicuo a v c h- Aqui, para obtener las 
relaciones de salto, ademas de tener en cuenta las ecuaciones consideradas en el caso de 
los choques transversales, usamos tambien la ecuacion V • B = 0. 

Resolviendo como antes estas ecuaciones se llega nuevamente a una relacion cubica para 
(. Una de las soluciones posibles es la trivial, ( = 1, en la cual no hay onda de choque. 
La solucion £ = 1 corresponde a una onda de Alfven si Vq = v\ que, como mencionamos 
antes, no son ondas de choque. Se las llama ondas intemedias porque se propagan a la 
velocidad de Alfven en la direction x, es decir, v c ^ = va 0x ■ Las otras dos soluciones, dadas 
por ( > 1, corresponden a las ondas de choque asociadas a las ondas magnetoacusticas 
rapida (v 0x > v Aox ) J lenta (v 0x < v Ao J- 
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Figura 2.5: Ondas de choque transversales (izquierda) y oblicuas (derecha). El ancho del 
choque es Ax — > y Bq y B\ son el campo magnetico en la region no perturbada y 
perturbada, respectivamente. 



2.4. Ondas de choque relativistas 

Un fluido que se mueve con una velocidad U > 0.1c es relativista. La information sobre 
su estado esta contenida en el tensor de energia- impulso T^ u . Para un fluido ideal: 

= ^ u u^Uy + P g^, (2.35) 

donde = T(c, u) es la tetra- velocidad, /3 = u/c y T = l/y/l — /3 2 es el factor de Lorentz 
y uj = e t + P es la entalpia. La densidad de energia total e t de un fluido relativista se 
escribe e t = nmc 2 + e, donde nmc 2 es la densidad de energia en reposo y e es la cinetica. 
Finalmente, es la metrica del espacio-tiempo. En esta tesis consideraremos la metrica 
de Minkowski, con g tt = — 1 y g xx = g yy = g zz = 1. En el SR propio, es decir, que se mueve 
con el fluido, T^ u es una matriz diagonal: 
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(2.36) 



El indice adiabatico para gases monoatomicos varia de la siguiente manera: 4/3 < 7 a( j < 
5/3. En el caso no relativista (NR) 7 a d = 5/3 mientras que en el caso ultra- relativista (UR) 
7ad = 4/3. En esta ultima situacion, la ecuacion de estado es P = e/3 y la velocidad del 
sonido resulta C s = c/ \^3. Las ecuaciones de movimiento y de conservacion de la energia 
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estan contenidas en la derivada <9T^/<9x M = mientras que la ecuacion de continuidad de la 
masa, o de conservation de la densidad de partfculas n, se obtiene a partir de d(nu u ) / dx u = 
0. 

La teorfa de las ondas de choque en fluidos relativistas se construye de la misma manera 
que para las NR presentadas en la Section [2731 (Taub 1948). Si consideramos un choque que 
se propaga en la direction x, las condiciones de salto se obtienen de considerar la conser- 
vation del numero de partfculas (nT/3), y las densidades de flujo de momento (uT 2 /3 2 + P) 
y de energia (uT 2 f3). En el SR del choque, estas ecuaciones se escriben de la siguiente 
manera: 

n T f3 = n 1 T 1 f3 1 , (2.37) 

u T 2 ft + P = u 1 Tj Pi + P u (2.38) 

UoTlPo^uxTlPt. (2.39) 

Para resolverlas consideraremos el caso UR, en el cual la ecuacion de estado es P = e/3 
y la entalpia se reduce aw = 4P. Ademas consideraremos que P <C u Tq /3q, es decir, un 
choque fuerte. En el SR del choque, el medio no chocado llega al frente de choque con una 
velocidad v o ~ c, con lo cual /3q ~ 1 y T — > oo. El medio chocado se aleja del choque a una 
velocidad menor que la del medio no chocado, pero que aun es moderadamente relativista: 

h = \ — ► ri = ^7f~ L06 - (2 - 40) 

Por otro lado, la densidad chocada resulta n x = 2\/2 no To- Pero recordemos que no esta me- 
dida en el SR del fluido. Si transformamos esta magnitud al SR del observador a traves de 
la relation n' = Tn, donde las magnitudes primadas son las medidas en este SR, obtenemos 

ni~3n . (2.41) 

Finalmente, la presion en el medio chocado resulta 

P^^mn'uToc 2 . (2.42) 

Al ser la presion un invariante relativista, Pi = P[. 

Si consideramos campos magneticos, entonces, al igual que en el caso no relativista, solo 
debemos hacer los siguientes reemplazos en las ecuaciones (12.381) y (I2.39[) : P — > P+B 2 / (87r) 
y e -^ + 57(871-). 



2.5. Produccion de ondas de choque 

Hay muchas maneras de producir ondas de choque en sistemas astrofisicos. En esta 
tesis, hemos estudiado aquellas que se forman por perturbar un fluido que se mueve con 
una velocidad mayor que la del sonido. Esto es, dado un choque que se propaga por un 
medio, estudiamos de que manera la interaction con el medio externo puede producir otros 
choques. 
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2.5.1. Ondas de choque por perturbaciones locales 

Una de las formas mas comunes de perturbar un fluido en movimiento es interponien- 
dole un obstaculo. Si el fluido es supersonico y el obstaculo es rigido, entonces se produce 
una unica onda de choque que se propaga en el sentido opuesto al movimento del fluido. 
Si el obstaculo no es rigido, se puede producir ademas un choque que se propaga a traves 
del obstaculo, como se muestra en la Figura l^Bl Un escenario posible para esta interaccion 
es un jet al cual se le interpone una inhomogeneidad del medio circundante. En esta tesis 
hemos estudiado dos situaciones similares: en un HMMQ y en un AGN. En el primer caso, 
estudiamos la interaccion de grumos presentes en el viento de la estrella companera con 
el jet del objeto compacto y en el segundo hemos estudiado la interaccion de nubes que 
circundan el agujero negro del centro de la galaxia con el jet del AGN. 

A fin de saber cuales son los procesos ffsicos mas relevantes en este tipo de interacciones, 
haremos un analisis de las escalas de tiempo de los procesos que tienen lugar debido a la 
interaccion de un obstaculo con un jet. 



Escalas de tiempo 

Supongamos un obstaculo esferico, de radio R Q , con una densidad homogenea n Q y que 
se mueve con una velocidad v Q en la direccion perpendicular a la direccion de propagation 
del jet, que designaremos z. Por otro lado supondremos un jet que se propaga con una 
velocidad t>j y con un factor de Lorentz iy Considerando un angulo de abertura del jet 
<f) ~ 6° (valor tfpico de jets colimados de MQs y AGNs), la relacion entre el radio y la 
altura es -Rj ~ 0.1z y Considerando que la luminosidad cinetica del jet, L } , es constante, la 
densidad, n, (en el SR del laboratorio), de este varia con Zj a traves de la relacion 

L i = ( r j - !) m P c2 °j n i v h ( 2 - 43 ) 

donde Oj = irr? es la section del jet a la altura Zj = 10 rj. Un parametro importante en 
nuestro estudio sera el cociente entre la densidad del obstaculo y la del jet: x = n /riy 

El tiempo de penetracion, t Q , del obstaculo en el jet puede estimarse a traves de la 
ecuacion: 

t ~ (2.44) 

Vo 

Por otro lado, si suponemos que dentro del jet el obstaculo se mueve siguiendo su trayectoria 
original, el tiempo que tardana en cruzar el jet (recorriendo su ancho igual a 2Ry) esta dado 
por 

H ~ IS ~ L^M. ( 2 .45) 

Una vez que el obstaculo ha entrado en el jet, la interaccion con el material de este 
ultimo producira dos choques: uno en el jet y otro (o una onda sonica) en el obstaculo, 
como se muestra en la Figura [231 El choque en el primero se mueve a una velocidad ~ uj 
y lo llamaremos bow shock ("choque de proa"). Al cabo de recorrer una distancia Z en el 
jet, este choque alcanza el estado estacionario, es decir, la separacion entre el bow shock y 
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Figura 2.6: Esquema de la interaction de un jet con un obstaculo (verde). Izquierda: el 
obstaculo penetra en el jet moviendose a una velocidad v Q . Derecha: Una vez que el obstacu- 
lo esta completamente adentro del jet, se forma un bow shock simetrico en el jet. Mientras, 
un choque se propaga en el obstaculo a una velocidad v co . El bow shock esta a una distancia 
Z del obstaculo. (El grafico no esta a escala.) 



el obstaculo es fija e igual a Z. Para estimar esta distancia consideramos que el numero 
de particulas del jet que por unidad de tiempo atraviesan el bow shock (~ njVj0~ o , donde 
a = it Rl es la section del obstaculo) y luego escapan siendo advectadas de la region 
chocada a la velocidad Vji del medio chocado del jet (fji ~ fj/4 6 fj/3, para los casos NR y 
UR, respectivamente) se conserva. De esta manera, de la igualdad ^jfj0" o = Wji fji 2tc R Z, 
obtenemos que Z ~ 0.2i? o y ~ 0.3i? o en los casos de choques UR y NR, respectivamente. 
Luego, el tiempo en el cual el bow shock alcanza el estado estacionario es 



Respecto del choque en el obstaculo, la velocidad de propagation puede determinarse 
asumiendo que se establece el equilibrio de presiones en la superficie de separation entre 
la region chocada del jet y del obstaculo y usando las relaciones de salto. Se obtiene la 
siguiente expresion para la velocidad: 



Z 

tbs ~ 




(2.46) 



CO 



(2.47) 



Luego el tiempo en el cual el choque recorre todo el obstaculo es 



t, 



2 R, 



o 




CO 




CO 



Esta magnitud t co puede considerarse un tiempo caracteristico (o de vida) del sistema, ya 
que como veremos a continuation, las inestabilidades pueden destruir el obstaculo en un 
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tiempo ~ t co o este podn'a ser acelerado y comenzar a moverse con el jet antes de escapar 
del mismo (en un tiempo ij). 

La aceleracion g que el material del jet ejerce contra el obstaculo puede determinarse 
considerando conservation de la entalpia y suponiendo que toda la presion cinetica del jet 
se convierte en presion termica del material chocado en la superficie de contacto con el 
obstaculo. Asi obtenemos que 

(Fi - 1) v? 

con lo cual el tiempo de aceleracion del obstaculo hasta la velocidad del jet es 

* g ~-~Vx*co. (2.50) 

Sin embargo, en un tiempo ~ t co el obstaculo se acelera hasta la velocidad v co del choque 
que se propaga en su interior. La aceleracion g ademas de acelerar el obstaculo, hace que 
se desarrollen inestabilidades de Rayleigh- Taylor (RT) en la parte de este que se encuentra 
mas cerca de la base del jet, como se muestra en la Figura 12.71 El tiempo en el cual estas 
inestabilidades crecen hasta una longitud de escala I ~ R esta dado por 




£rt ~ W - ~ W — ~ t co . (2.51) 



Por otro lado, la diferencia relativa v Te \ entre la velocidad del material chocado del jet que 
rodea al obstaculo, v$i, y la velocidad del material chocado del obstaculo, v Q i, desarrolla 
inestabilidades de Kelvin-Helmholtz (KH), con una escala de tiempo dada por: 

tKH „ L£ „ (2.52) 

donde nuevamente hemos considerado I ~ R Q . Notar que como v Te \ ~ v^, ^kh ~ t co . En las 
estimaciones de £rt y de ^kh no hemos tenido en cuenta al campo magnetico, que podria 
estabilizar el sistema (Blake 1972, Romero 1995). 

De las diferentes escalas de tiempo estimadas anteriormente, podemos resumir que una 
vez que el obstaculo penetra en el jet, el bow shock alcanza la distacia Z rapidamente, 
en un tiempo tbs <C t co y luego el obstaculo podria ser destruido por las inestabilidades, 
acelerarse y comenzar a moverse conjuntamente con el jet o escapar de este. Sin embargo, 
para hacer un analisis mas exaustivo de cuales seran los procesos mas relevantes necesitamos 
especificar la localizacion en el jet a la cual se produce la interaccion. La determinacion de 
este parametro, conjuntamente con otras magnitudes del sistema (jet y obstaculo), seran 
expuestas en los Capitulos[5]y[6l donde los escenarios considerados seran HMMQs y AGNs, 
respectivamente. 



2.5.2. Ondas de choque por interacciones globales 



Cuando un fluido supersonico choca con un medio externo, la onda de choque se disipa 
a medida que se propaga si la inyeccion de energfa es puntual en el tiempo, como ocurre 



2.5 Produccion de ondas de choque 



31 




Figura 2.7: Diagrama del obstaculo y del bow shock. Las inestabilidades de RT se desarrollan 
en la superficie frontal del obstaculo (figura izquierda), mientras que las de KH lo hacen 
en las superficies laterales (figura derecha). 



en las explosiones de supernovas. Por otro lado, si la inyeccion de energfa es continua el 
choque se frena, pero mas lentamente que en el caso anterior. Este choque llega a un estado 
estacionario en el caso radiative En las dos situaciones descriptas anteriormente se produce 
un choque reverso que se propaga en el sentido opuesto al movimiento del fluido. En el 
caso del evento puntual, el choque reverso al cabo de un cierto tiempo se disipara. Si la 
inyeccion inicial de energfa es continua, se formara un choque estacionario. 

Situaciones en las cuales pueden producirse choques de las maneras descriptas anterior- 
mente pueden darse cuando los jets son frenados por el medio en el cual se propagan y se 
producen los llamados choques terminales, como ocurre por ejemplo en los jets de YSOs 
masivos y en los hot-spots de las radio-galaxias. Otra situation en la cual se desarrollan 
choques es cuando dos cumulos de galaxias colisionan. Estos choques a gran escala se pro- 
pagan hasta una cierta distancia del centro de gravedad de la fusion de ambos cumulos y 
luego se disipan. 

A diferencia de lo expuesto en la section 12.5.11 aqui estudiaremos dos sistemas con 
propiedades muy diferentes: el jet de un YSO y un cumulo de galaxias. Es por esto que la 
description de las escalas de tiempo que haremos a continuation sera muy general y luego, 
en los capftulos correspondientes a cada fuente particular: YSOs (Capftulo Hj) y cumulos 
de galaxias (Capitulo Ej), haremos un estudio mas detallado de las escalas de tiempo con 
formulas pertinentes a cada fuente. 

Escalas de tiempo 

Supongamos un choque que se propaga a una velocidad en un medio externo. El 
choque, a medida que se propaga empuja y apila material delante de el, y asi el frente de 
choque es ahora la "cabeza" de ese material que se mueve a la velocidad de la onda. Defini- 
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mos el tiempo de propagation t p como aquel en el cual el choque se frena sustancialmente, 
es decir, podemos suponer que durante t p este se propaga "libremente" por el medio una 
distancia Z p en un tiempo: 

t p ~ ^. (2.53) 

Vch 

Sin embargo, mientras mas material se haya acumulado, mas dificil sera moverlo y asi la 
velocidad del choque va disminuyendo. Si la velocidad decrece, entonces la luminosidad 
tambien (de acuerdo a la ecuacion (12.431) ) y es asi como el choque se va debilitando. Si v c h 
disminuye, entonces hay fuertes variaciones de presion en la region chocada y se produce 
una onda. Si esta se propaga a una velocidad mayor que C s (6 C s a) entonces sera una onda 
de choque y la llamaremos choque reverso. Este choque reverso se disipa o llega al estado 
estacionario si la inyeccion de energfa es puntual o continua, respectivamente. En ambos 
casos, el choque reverso recorre una distancia Z T a la velocidad v r en un tiempo 

t r -. (2.54) 

v T 

La velocidad v r puede estimarse a traves del contraste de densidades entre ambos medios 
chocados. Entre las regiones de ambos medios perturbados (el externo y el chocado por el 
choque reverso) se genera una superficie de discontinuidad, como se muestra en la Figu- 
ra 12.81 En esta superficie, las inestabilidades de RT pueden ser importantes y mezclarse 
material de ambos medios chocados. La escala temporal de esta inestabilidad esta determi- 
nada por la ecuacion ( 12. 51ft : £rt ~ (l/g) 1 ^ 2 , pero considerando I ~ Z r y que la aceleracion 
es ejercida o bien por el material chocado del jet o bien por el centro del cumulo de galaxias 
(en los casos particulares que estudiamos en esta tesis). 

Finalmente, el choque reverso puede ser radiativo si l te r < Z T , siendo l teT ~ t t crV T y el 
tiempo de enfriamiento por emision termica 

iter ~ (2.55) 

* *-ter 

La funcion £ ter depende de la temperatura del medio chocado a traves de la funcion A(T), 
como se indica en la ecuacion (12. 29ft . 

2.6. Aceleracion de particulas en choques 

Las ondas de choque, ademas de modificar las propiedades termodinamicas del medio 
por el cual se propagan, pueden tambien acelerar particulas hasta energias relativistas. 
Estas particulas aceleradas son llamadas particulas no termicas, ya que su distribucion 
en energia no es Maxweliana. El numero de estas particulas no termicas es una fraccion 
pequena del numero de particulas termalizadas que forman el medio en el cual se propaga el 
choque. Es decir, el choque acelera solo una fraccion peque na de las particulas termicas del 
medio sacandolas del estado de equilibrio termodinamico y convirtiendolas en no termicas. 
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Z r 




Superficie de discontinuidad 



Figura 2.8: Esquema de un choque que se propaga en un medio externo con una velocidad 
v c h. Se forma un choque reverso que se propaga con una velocidad v T en sentido opuesto a 
v c h, y entre ambos medios chocados se forma una superficie de discontinuidad en la cual 
pueden desarrollarse inestabilidades de RT. 



34 



Ondas de choque y aceleracion de particulas 



La presencia de estas particulas en fuentes astrofisicas es usualmente detectada a traves 
de la emision no termica en radio y de la radiacion gamma (la cual no puede ser termica). 

El mecanismo propuesto originalmente por Fermi (1949) describe como las particulas 
podrian ser aceleradas por rebotes sucesivos entre centros dispersores con una distribution 
aleatoria de velocidades. Este proceso resulta poco eficiente ya que la ganancia de energia 
E cada vez que la particula rebota es (AE)/E oc /3 2 , y como /? <C 1 la ganancia es muy 
poca (/3c es la velocidad de la particula). Por otro lado, ademas de ser poco eficiente, 
la situation en la cual tiene lugar el mecanismo es poco frecuente en la naturaleza. Sin 
embargo, es facilmente adaptable a un escenario astrofisico si reemplazamos los centros 
dispersores por un choque e inhomogeneidades magneticas en el medio. La teoria original 
de Fermi fue modificada en los afios 70 por diversos autores, entre ellos Axford, Lear & 
Skadron (1977) y Bell (1978) para describir un mecanismo mas eficiente ((AE)/E oc 0) y 
que tuviese lugar en sistemas astrofisicos con choques. Para que este proceso tenga lugar, es 
necesario que las particulas en el medio perturbado puedan difundir y alcanzar el choque. A 
continuation describiremos someramente el mecanismo de aceleracion de Fermi modificado 
y que se conoce como mecanismo de Fermi de tipo I. Este proceso tambien se conoce 
como mecanismo de aceleracion difusiva ya que la difusion juega un rol importante en el 
desarrollo del mismo. 

En un plasma en el cual el campo magnetico es uniforme, las particulas cargadas se 
mueven siguiendo trayectorias helicoidales alrededor de las lfneas de B. El radio de giro 
(r g ) esta determinado por la energia E de las particulas y por la intensidad del campo 
magnetico: 

E 

rv = , 2.56 

donde Z e es el mimero atomiccQ y q e es la carga electrica del electron. Sin embargo, en los 
sistemas astrofisicos que nos interesan el campo B presenta irregularidades que perturban 
el movimiento helicoidal de las particulas produciendo una reorientation de las mismas, y 
por lo tanto un movimiento desordenado. Las particulas con velocidad v difunden en el 
medio con un camino libre medio A = 3D/v, donde D es el coeficiente de difusiorjf). Si el 
movimiento de la particula es en una direccion que forma un angulo 9 con B, entonces 
D = D|| cos 2 (#) + D± sen 2 (9), siendo D\\ y D± los coeficientes de difusion en la direccion 
paralela y perpendicular a B, respectivamente. Debido a que D es desconocido en la mayoria 
de los sistemas astrofisicos, se suele considerar que D» es un mimero rj veces el coeficiente 
de difusion minimo o de Bohm, Db = r g c/3, y que significa considerar que el camino libre 
medio en la direccion paralela a B es An ~ r g . Luego resulta D± w E)\\/{1 + rj 2 ) (Jokipii 
1987). 

Con cada ciclo en el cual la particula va del medio no chocado al chocado y vuelve al 
medio original, la ganancia de energia es (AE)/E ~ (4/3)/3. Al cabo de k c ciclos, la energia 



En todas las aplicaciones realizadas durante esta tesis Z e = 1. 
5 En el Capftulo [3] se dara la definition del camino libre medio. 
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de la particula, que inicialmente era resulta 



AE X k 



E = Ei\l + —) . (2.57) 

El tiempo en el cual se realiza cada ciclo esta determinado por como las particulas difunden 
en cada medio, es decir, de los coeficientes de difusion D Q y D 1 , a traves de la expresion 
^cicio ~ {^/c){Do/vq + Di/vi) (Protheroe 1999). El "tiempo de aceleracion" es el tiempo 
requerido para que las particulas alcancen una energia E a traves de un mecanismo de 
aceleracion dado; para el mecanismo de Fermi de tipo I resulta: 



1 dE 
E~dt 



(AE) Lcicl ° (C-l)v \v • Vl 



— + — . (2.58) 



Asi, el tiempo de aceleracion depende de la geometria de B. A continuacion damos expre- 

— * 

siones sencillas para t ac correspondientes a configuraciones particulares de B. 

■ Choques paralelos (B \\ v ch ): suponiendo que D — D ± — D\\ y que B = Bi se tiene 
que el tiempo de aceleracion es 

6 q e JJ v ch 

■ Choques transversales (B _L v ch ): suponiendo que D = D 1 = D± y que B\ ~ AB se 
tiene que el tiempo de aceleracion es 

~ l^rk^- (2 - 60) 

Por otro lado, debido a que en los choques relativistas el calculo de D es aiin mas complicado 
que en los NR, se puede considerar que en los primeros t ac ~ E/(0.1 q e B c) (como en el 
caso de la supernova del Cangrejo - Crab-), sin especificar la geometria de B. 

El numero de ciclos que realiza la particula depende basicamente del tamaho L ac del 
acelerador. Las particulas pueden cruzar el choque sucesivas veces antes de que r g (oc E) 
crezca lo suficiente, es decir, r g ~ L ac , como para escapar del acelerador (Hillas 1984). Al 
cabo de k c ciclos el espectro de las particulas aceleradas e inyectadas en el medio chocado 
sera Q(E) = K E~ P , donde 

p= 1 < E ' E ^ (261) 
P \n(l + AE/E) [Z - bi) 

y [Q] = erg -1 s _1 . En la ecuacion anterior la dependencia con k c esta imphcita en E/Ei. 
Esta inyeccion Q(E) de particulas por unidad de tiempo en una determinada region de 
la fuente dara lugar, al cabo de un tiempo t, a una distribucion N(E, t) de particulas no 
termicas. 
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2.6.1. Poblacion de particulas no termicas 

Dada una distribution Q(E) de particulas aceleradas e inyectadas en una region, nos 
interesa conocer cual sera la distribution de particulas N(E, t) en un tiempo t. 

A lo largo de esta tesis consideraremos que tanto la region donde se aceleran las particu- 
las, el acelerador, como aquella en la cual radian, el emisor, son homogeneag^). Es decir, 
las magnitudes fisicas no dependeran de las coordenadas espaciales. Bajo esta hipotesis, 
la ecuacion que describe la evolution temporal de N(E, t) es la siguiente (Ginzburg & 
Syrovatskii, 1964): 



dN(E,t) d fdE 
dt ~ dE V ~d7 



N(E,t) ) -^2^. + Q(E), (2.62) 



r; „: / *esc 



donde [N(E,t)] = erg -1 . Esta ecuacion tiene en cuenta la inyeccion de particulas Q(E) 
que en un tiempo t nos da la distribution N(E,t). Las particulas desde que se inyectaron 
en el tiempo to = hast a el tiempo t sufren perdidas radiativas que estan contempladas 
en el termino dE/dt\ Ta ,d = E Ta &. Ademas, puede haber perdidas de energia en el sistema 
debido a que las particulas pueden escapar del acelerador. El escape de las particulas 
esta considerado en el termino N(E,t)/t esc , donde 

= (t~ + ■ (2-63) 

Esta escala de tiempo tiene en cuenta el escape por convection (t conv ) y por difusion (t dif ) . 
Mientras que el primero contempla el arrastre de las particulas por el movimiento del 
medio chocado (que se mueve a una velocidad ~ t> c h/4), el segundo considera las perdidas 
de energia por difusion de las particulas. 

La solution mas general de la ecuacion f )2.62p considerando que t esc = t conv , es decir, 
^conv <S tdiff ? es (Khangulyan et al. 2007) 

1 f EeS ( t(E E')\ 

N(E, t) = -r- Q(t, E') exp - \' ' dE?, (2.64) 

| E I J E V ''esc / 

donde E e g es la energia maxima de las particulas que en el tiempo t pueden enfriarse hasta 
una energia E, es decir, 

m = J, W{ y r(E,E^j E w . (2.65) 

Dada una energia E', r es el tiempo en el cual una particula de energia E' se enfria hasta 
obtener una energia E. Notemos que r < t. La ecuacion (I2.64[) toma una expresion muy 
sencilla en dos casos limites que veremos a continuation, si p > 1: 



6 Notamos que el campo magnetico debe tener una componente inhomogenea para que se el mecanismo 
de Fermi pueda desarrollarse. 
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■ Si t(E) <C t csc entonces t(E,E') <C t esc y entonces el termino exp(— r/t csc ) ~ 1, con 
lo cual N(E,t) ~ (-Q(E) + Q(E eS ))(E cS - E)/\E\ ~ Q(E) t(E) y las parti'culas se 
acumulan a medida que t crece. Este proceso de acumulacion de particulas de energia 
E continua hasta que E e g coincide con la energia maxima de la distribucion Q(E) o 
el tiempo de escape se hace menor que el tiempo de enfriamiento. 

■ Si t csc t(E): N(E) ~ Q(E) t csc y el sistema llega a un estado estacionario. 



Para conocer la distribucion N(E) necesitamos conocer las perdidas radiativas que 
sufren las particulas. Para esto, en el proximo capitulo repasaremos los procesos radiativos 
que hemos estudiado durante esta tesis. 



38 Ondas de choque y aceleracion de parti'culas 



Capitulo 3 
Procesos radiativos 



Todo lo que conocemos de la mayorfa de las fuentes astrofisicas es gracias a que ellas 
radian y a que es posible detectar parte de los fotones emitidos con instrumentos adecua- 
dos. Por esta razon, conocer los procesos a traves de los cuales se genera la emision es 
indispensable para modelizar correctamente las fuentes que nos interesa estudiar. 

Como mencionamos en el Capitulo [U, la emision de rayos gamma no es posible mediante 
procesos en equilibrio termodinamico. La produccion de estos fotones tan energeticos ocurre 
a traves de procesos no termicos, esto es, por la interaccion de partfculas relativistas con 
campos de materia, de fotones y magneticos. Todas estas interacciones pueden describirse 
como dispersiones de las partfculas relativistas por la interaccion con otras partfculas (en 
ocaciones virtuales) y que como resultado producen fotones. Por esto, para estudiar la 
produccion de radiacion (gamma, en particular), es conveniente definir algunos conceptos 
basicos que son comunes a todos los procesos radiativos no termicos y que usaremos a lo 
largo de toda la tesis. 



3.1. Conceptos basicos 

Dada una poblacion de partfculas de tipo i contenidas en un volumen V, si todas ellas 
tienen energfas E{ diferentes tales que Ef 1111 < E( < _E 4 max , la densidad total de estas 
particulas sera 

f E ^ dn- 

rk= I j^dE t , (3.1) 



mm 



donde [n*] = cm 3 . La densidad de energfa de estas partfculas es 



* dn ■ 

u t = I ^E.dE,, (3.2) 



siendo [ui] = erg cm 3 . Si la distribucion fuese monoenergetica, es decir, Ei constante, 
entonces Uj = rii E^ Si la velocidad de las partfculas es v iy la luminosidad de estas sera 



L{ — Ui ViS, 



(3.3) 
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donde la superficie S = dV y [Li] = erg s _1 . Si las particulas son relativistas, con una 
energfa cinetica (1^ — 1) to, c 2 , la ecuacion (13.31) se escribe como (I2.43p . 

Supongamos que inyectamos una partfcula relativista, con energfa E i: en un medio en 
el cual existe una densidad de particulas que llamaremos bianco y cuya energfa es E h . 
En el SR de la partfcula relativista, las particulas del medio se mueven a una velocidad 
~ c y el numero de interacciones que ocurren por unidad de tiempo es 

N int = a{Ei)n h c, (3.4) 

donde a(Ei) es una medida de la superficie efectiva ([a] = cm 2 ) de interaction del proceso 
de dispersion y conocida como seccion eficaz0. Mientras mayor sea a, mayor es la proba- 
bilidad de dispersion. A la energfa de los fotones creados en la interaction la llamaremos 
E p h y depende de la energfa initial E h , de la energfa de la partfcula relativista, E iy y del 
angulo de interaction. Sin embargo, E p ^ no esta completamente determinada ya que de- 
pende de las caracterfsticas de la colision. Por lo tanto, el estado final, caracterizado por 
Eph y por el angulo de dispersion Q, tendra asociada una fruition de probabilidad. Esta 
information esta contenida en la seccion eficaz diferencial da(E p ^, E i ,Q)/dE p i l dQ, de tal 
manera que la seccion eficaz total es 




dV(ff ph ,E t ,^ F 



dE ph dtt 

Debido a la naturaleza del proceso, la distribution energetica de los fotones creados o 
dispersados dependera fuertemente del angulo de interaction. Sin embargo, ya que en 
la mayorfa de las situaciones astroffsicas este angulo no se conoce, se suelen tomar las 
magnitudes isotropizadas, es decir, integradas en todos los posibles angulos de interaction 
(o de dispersion). En esta tesis no consideraremos las dependencias angulares (ni espaciales: 
emisor homogeneo) de las magnitudes. De la ecuacion (13.4)) . si x = ct = (unb) -1 , entonces 
se tiene que ocurre solo una interaction en el tiempo t. Se define el camino libre medio 
como la distancia A que las particulas recorren entre dos interacciones sucesivas: 

A EE — lr . (3.6) 

El tiempo de enfriamiento t ra( j es una medida de la eficiencia de un proceso de 
interaction de una partfcula relativista con un bianco de densidad n h . Si en este proceso 
la partfcula relativista pierde una cantidad de energia Ei — Ef, donde Ei y Ef son las 
energfas initial y final, respectivamente, y definimos la inelasticidad / del proceso como 
/ = (Ei — Ef)/Ei, entonces 

^rad = , , p > • (3.7) 



1 Debido a que los valores de a son muy chicos, la unidad de medida que se suele usar es el Barn: lb 



lO" 24 cm 2 . 
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Las perdidas de energia que mencionamos en el capitulo anterior, en la ecuacion (I2.62p . se 
definen ahora como 



dE, 



dt 



Ei 

= -. (3.i 

rad ^ r ad 



Supongamos que inyectamos una distribution Qi(Ei) de particulas relativistas en el 
medio. Al cabo de un tiempo > t r ad/esc, la distribution de estas sera Ni(Ei) ~ Qi(Ei) t ra( j/. 
El numero de estas particulas acumuladas en el volumen V sera3 



esse • 



rrel 



N?= / NiiEjdEi. (3.9) 

Para conocer el espectro de fotones producidos en la interaction de las N™ 1 particulas 
relativistas y cuya distribution en energias es Ni(Ei), necesitamos conocer la emisividad 
Jph de la fuente. Si las particulas son relativistas podemos suponer que se mueven a una 
velocidad ~ c y entonces 

Jph(£ P h) ~ ^ / ' n h da[ 7» h > Ei} N^Ei) dE u (3.10) 
4irJ Ei=Eph dE ph 

donde [Jph] = erg -1 s _1 . Si conocemos J p h luego la luminosidad especifica en el caso de un 
emisor homogeneo e isotropico es 

L Eph ~E ph J ph (3.11) 

A continuation haremos una somera descripcion de los cuatro procesos radiativos que 
hemos considerado en esta tesis: radiation sincrotron, dispersion Compton inversa (IC), 
Bremsstrahlung relativista e interacciones proton-proton (pp). 



3.2. Procesos radiativos 

Existen varias maneras de producir fotones. Una de ellas es a traves de la aceleracion 
de particulas cargadas a traves de la fuerza ejercida sobre ellas por algun campo. Otra 
manera es a traves del decaimiento de particulas. 

3.2.1. Radiacion sincrotron 

Las particulas cargadas sienten la fuerza de Lorentz que ejerce sobre ellas el campo 
electromagnetico. Esta fuerza hace que las particulas describan un movimiento helicoidal 

2 A las magnitudes intensivas las llamamos con letras minusculas y a las extensivas con mayusculas. 
Por ejemplo, las distribuciones de energia de las particulas seran m(Ei) si son por unidad de volumen y 
Ni{Ei) en caso contrario. De esta manera, bajo la suposicion de que m(Ei) no depende de las coordenadas 
espaciales, Ni(Ei) = rn(Ei) V. 
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F(x) 




x 



0.29 1 2 3 4 5 6 



Figura 3.1: Grafico de la funcion F(x), donde x = E ph /E c . Como puede apreciarse, el 
maximo corresponde a x = 0.3, esto es, E p ^ = 0.3 E c . Gracias a esta forma tipo funcion 
5 que muestra P sin (E ph , E e ) es posible adoptar la aproximacion de que los electrones de 
energia E e emiten fotones de energia ~ 0.3E c (E e ). 



alrededor de las lmeas de campo magnetico B. Las particulas relativistas que son aceleradas 
por la fuerza de Lorentz producen la llamada radiacion sincrotron. Como veremos luego, 
este proceso radiativo es mucho mas eflciente para leptones que para hadrones, por lo cual 
en esta tesis solo consideraremos la radiacion sincrotron producida por los primeros. Otra 
caracteristica de este proceso es que los fotones son emitidos en una direction preferential, 
aquella correspondiente al movimiento de la particula. 

La radiacion es emitida en un cono cuyo angulo de apertura llamaremos 9 ~ l/r e , donde 
T e = E e /(m e c 2 ) es el factor de Lorentz de un electron con energia E e . La distribution en 
energia de la potencia sincrotron por particula es 

P sm (E e , E ph ) = l^B ± ^ ^ K 5/3 (() dC, (3.12) 
h m e c z E c J Eph /E c 

donde B± es la componente de B perpendicular a la direction del movimiento de la particu- 
la. Si asumimos que B es isotropico, entonces Bj_ = a/2/3S. La energia caracteristica de 
los fotones producidos es E c (E e ) = 5.1 x 10~ 8 B E 2 erg. La funcion de Bessel de segunda 
especie y de orden 5/3, K 5 / 3 , tiene una forma tal que la integral de ella por el cociente 
E e /E c puede aproximarse de la siguiente manera 

§^ r iMCK - 1-85 (§^) 1/3 -p i^) ■ (3 - i3) 

La funcion P sin (E e , E ph ) tiene un maximo muy pronunciado en la energia E ph ~ 0.3^, 
como se muestra en la Figura 13.11 

La perdida total de energia radiada por un electron de energia E e se obtiene integrando 
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la distribution (13. 12 j) en todas las energias radiadas E ph : 



dE e 



P sin (E e ) = - / P siri (Ee, E ph ) dE, 



(3.14) 



dt 



ph- 



sm 



Luego, el tiempo de enfriamiento resulta 



2tt (m e c 2 ) 2 
3ca T E e B 2 



4.1 x 10 2 



B 2 E e 



s . 



(3.15) 



Siendo t sin oc m 4 , donde m es la masa de la particula que esta radiando, el tiempo de 
enfriamiento por radiacion sincrotron sera ~ 10 12 veces mas corto para electrones que para 
protones, ya que m e ~ lCT 3 m p . 

Si ahora consideramos que inyectamos una poblacion Q e (E e ) oc E~ p de electrones 
relativistas, las perdidas que sufriran por radiacion sincrotron modificaran Q e (E e ) dando 
lugar a la distribution de electrones 



Dada la distribution N e (E e ), para obtener la potencia radiada por todos estos electrones 
lo que debemos hacer es integrar el siguiente producto: 

gmax 



donde E™ m y E™ ax son las energfas minima y maxima, respectivamente, de la distribu- 
tion de electrones relativistas. Sin embargo, debido a que cada electron radiara la mayor 
parte de su energfa E e en fotones de energfa E p ^ ~ 0.3E c (E e ), podemos calcular la inte- 
gral (13.171) considerando una aproximacion tipo 5. De esta manera, utilizando una aproxi- 
macion £(-E pn — 0.3E C ) podemos hallar una solution analftica aproximada de P S i n (£ , p h). 

Para obtener la luminosidad a la energfa E ph calculamos E ph L ph = E ph P S i n (E ph ) . Es 
interesante notar que el espectro de fotones emitidos tiene una distribution en energfa tipo 
ley de potencias de fndice a (E^Ls^ oc E~^), si el espectro de electrones tambien es una 
ley de potencias N e oc E~ p ' . Las relaciones entre ambos indices es a = (p' — l)/2. 

3.2.2. Radiacion Compton inversa 

Consideremos un gas de fotones de energfa -E p h m cuya densidad de energia es M p h m y 
un electron relativista que atraviesa dicho gas. Si en el SR del laboratorio los fotones son 
menos energeticos que el electron, entonces los primeros seran dispersados por el segundo. 
Como producto de esta interaction los fotones ganan energfa, a diferencia de lo que ocure 
en la interacciones Compton (directas), y de aqui que se las llame dispersiones Compton 
inversas. 

En el SR del electron, si se acerca un foton poco energetico, este sufrira un cambio 
mayor en su energia y momento que uno con mas energfa. Esto se refleja en la section 




(3.16) 




(3.17) 
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eficaz de la interaction, <ji C . Si integramos en todos los posibles angulos de incidencia, die 
depende de la energfa de los fotones semilla y del electron a traves de la expresion (Vila & 
Aharonian 2009) 



vie 



3o"x 



y 



, , ln(l + 2y) 
V J 



1 
2 



y 2(1 + 2^ 



(3.18) 



donde y = E p ^ m E e / '{m 2 e c A ) y or ~ 0.67 b es la section eficaz de Thomson. Notamos que 
para j/ < 1, crjc ~ cr T y se dice que la interaction ocurre en el regimen de Thomson 
(Th), mientras que si y 1 entonces crjc ~ (3ctt/8) ln(Ay)/y y se dice que interaction 
se desarrolla en el regimen de Klein- Nishina (KN). Como se muestra en la Figura l3~2l 
en el regimen de KN aic cae abruptamente. Sin embargo, las perdidas de energia por 
interaction son mayores en el regimen KN que en el Th, siendo en el primero catastroficas, 



E 



P h 



E e . En el regimen Th, la energia maxima que pueden alcanzar los fotones 

2 



esto es 

es £ p m h ax ~ (4/3) r e 2 E ph __ 

Dado un electron relativista inmerso en un gas de fotones, el numero de interacciones 
por unidad de tiempo puede calcularse a traves de la ecuacion ( 13 .4p . considerando que 
nb = n phm , con lo cual obtenemos que el numero de interacciones por unidad de tiempo 
es iVi n t = n p hm cic c. El numero de interacciones es igual al numero de fotones dispersados. 
Si ahora queremos conocer la distribution energetica de los fotones dispersados por un 
electron, debemos tener en cuenta que la energia -Eph de estos fotones no esta fija sino que 
es una distribution y esto se muestra en la section eficaz diferencial (Blumenthal & Gould, 
1970) 

daic f „ p , 3cr T 

"l-frph, tie) 



dR 



(3.19) 



con 



fix* 



2 (4e ph r e 



2x ic In x ic + x ic +l 2x( c + 1 ^ ' + 4 ' e ' ph^ea . ic) , ' 



Pu 



ATI 



1, x 



1 J- 5 -^lC I ) 



(3.20) 



donde e p h = E ph /m e c 2 . La funcion P ic es igual a 1 cuando I /ATI — x ic < 1 y -Pic = en 
cualquier otro caso. Siendo e p h m = E phra /m e c 2 , el parametro adimensional x ic se define de 
la siguiente manera 

= 7- ~ 7t-tt - (3.21) 



4e phm r2(l - e ph /r e )' 

La interaction de un electron con un campo de fotones dispersa a estos dando un 
espectro de fotones dispersados con energias E™™ < E ph < E™ x . Para calcular la potencia 
radiada por un electron por interacciones IC debemos integrar sobre todas las energias 
P p hm a traves de la expresion 



Pic(E e 



dE 
~dt 



ic 



phm 



E nun 
phm 



Ph 



E mm 
ph 



daic(E ph ,E e ) dn phm (£ , phr 



dE, 



P h 



dE, 



phm 



E p h dEphm di? p h- 

(3.22) 
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Figura 3.2: Section eficaz de la interaction IC en los regfmenes Th y KN para fotones 
semilla con energfa E phm = 1 eV. 



Para el caso en el cual los fotones semilla siguen una distribution (quasi) monoenergetica con 
energia -Ephm, las perdidas de un electron pueden obtenerse a traves de una aproximacion 
^(-Ephm ~~ -Ephm)- El tiempo de enfriamiento tanto en el regimen de Th como en el de KN 
queda determinado por la expresion (Bosch-Ramon & Khangulyan, 2009) 

f = 6.1xl0 12 E phm (1 + 8.3?/) (l + l-3y 2 ) 

10 Wphm ln(l + 0.2y)(l + 0.5y + 1.3i/ 2 ) ' l ' ' 

donde ahora y = E phm E e / (mlc 4 ) y w p hm es la densidad de energia de los fotones ambienta- 
les. Considerando los hmites a bajas energias de (I3.23p . es posible hallar t IC en el regimen 
de Th, result ando 

tj c h ~ s- (3-24) 

- C/ e ""phm 

Si ahora inyectamos una distribution Q e (E e ) cx E~ p de electrones relativistas en el gas 
de fotones, las perdidas por IC en el regimen Th modifican el espectro de la misma manera 
que las perdidas por radiation sincrotron, esto es, 

N e (E e ) ~ Q e (E e ) tj£ cx E~ p ~ x . (3.25) 

Dada la distribution N e (E e ), para obtener la distribution de fotones dispersados calculamos 
la emisividad de los mismos a traves de 

pmax run ax , 7 _ <_ , , 

47r Jcmin JE™™ Cl£/ph Cliiphm 

e phm 
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Auto Compton 

Si los fotones semilla son creados externamente a la fuente (es decir, fuera del volumen 
V) entonces el proceso se dice Compton externo (EC, por External Compton). Por otro lado, 
en el caso particular de que los fotones semilla sean producidos por radiacion sincrotron 
de la misma poblacion de electrones N e (E e ) que interactiian por IC, entonces el proceso se 
llama auto Compton sincrotron (SSC, por Synchrotron Self Compton). Bajo esta situation, 
la distribution de fotones semilla es una ley de potencias y para calcular el espectro de 
los fotones producidos debemos usar la ecuacion (I3.26p . considerando que dn phm / dE phra = 
dn S i n / dE sin es la densidad de fotones producidos por radiacion sincrotron. 



3.2.3. Bremsstrahlung relativista 

Un electron relativista inmerso en un campo de materia (no relativista) sera acelerado 
por el el campo coulombiano producido por los nucleos de los atomos que forman este 
ultimo. La materia puede estar constituida por nucleos desnudos (es decir, ionizada) o por 
atomos (nucleos apantallados por los electrones). En ambos casos, la interaction puede 
esquematizarse de la forma 



e + (e,N) — y e + (e,N) +7, 



(3.27) 



aunque la section eficaz no es la misma. Esta viene dada por la expresion (Bosch-Ramon 
2006) 

da Brem (E e , E ph ) Aa ei r 2 Z 2 

-(p{E e , E ph ), 



dR 



P h 



(3.28) 



P h 



donde a e f ~ 1/137 es la constante de estructura fina y r e = q 2 /m e c 2 es el radio clasico del 
electron. Para el caso de un nucleo desnudo, la funcion (f> esta dada por 



<t>{E e , E ph ) 



E, 



P h 



E e 



E p h 
~E~, 
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2E e (E e — E ] 



ph) 



m e c 2 Ep h 



(3.29) 



mientras que para el caso en que el nucleo esta completamente apantallado por todos los 
electrones 



(f>(E e , Eph) 



P h 



E„ 



E, 



P h 



E e 



In 
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^V3 
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E, 



ph 



E„ 



(3.30) 



El electron entrega casi toda su energia a los fotones en estas interacciones, es decir, 
Eph ~ E e , y por esto se dice que las perdidas que sufre son catastroficas. Sin embargo, 
es posible hallar una expresion continua para el tiempo de enfriamiento. Para calcularlo, 
hacemos el mismo analisis que en el caso de las interacciones IC, sin embargo, debido a que 
hay dos parametrizaciones de crerem, para el caso en el que la materia esta apantallada o 
aquel donde el medio esta completamente ionizado, obtendremos dos expresiones diferentes 
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para el tiempo de enfriamiento. Siendo n m la densidad del campo de materia, en el primer 
caso tenemos 

1.4 x 10 16 

Wn ~ 7 7— \ 7T" S ' ( 3 - 31 ) 

n m Z? f In ' 183 * 1 



m e V U e 1/3 7 is 
mientras que en el segundo 

ion I- 4 x 10 16 

^Brem ~ 7 7 ~ \ T S - ( 3 - 32 ) 

^^ 2 (ln(^)+0.36j 

Si ahora inyectamos una distribucion Q e (E e ) oc de electrones relativistas, las perdi- 
das por las interacciones con la materia casi no modifican Q e (E e ) ya que la dependencia 
con E e de ^Brem es despreciable. Por esto, 

N e (E e ) ~ Q e (£ e ) t B rem oc E~* '. (3.33) 

La emisividad producida por una distribucion de electrones N e (E e ) esta dada por la 
ecuacion (13. lOj) . donde n D es la densidad de materia y a = 0Brem, con lo cual hallamos 

4tt y Eph dE ph 



3.2 .4. Interacciones proton-proton 

Ademas de procesos asociados a cargas aceleradas, los fotones tambien pueden producir- 
se por decaimientos de partfculas. Por ejemplo, los piones neutros 7r° decaen en dos rayos 
gamma. Una manera de producir 7r° es a traves de colisiones inelasticas proton-proton, 
en las cuales un proton relativista interactua con un proton no relativista. El canal de 
interacciones pp mas importante para generar 7r° es el siguiente (Romero 2010): 

p + p — >p + A + + n 7r° + h±(n + + (3.35) 

donde tc + y ir~ son piones cargados positiva y negativamente, respectivamente, y no y n± 
son las multiplicidades (numeros enteros positivos). En este proceso, el proton relativista 
pierde ~ 50% de su energia, con lo cual la inelasticidad de la interaccion es f pp ~ 0.5. 
El decaimiento A + — > p + 7r° entrega al 7r° aproximadamente un 17 % de la energia E p 
del proton relativista. El 33 % restante de E p es entregada al resto de los piones creados. 
Un 11 % va a los piones neutros y un 22 % a los cargados. Sin embargo, debido a que las 
multiplicidades h y n± son grandes, la energia de cada una de estas particulas es poca. 
Es por esto que al 7r° con energia ~ 0.17£ , p se lo llama "pion lider". 

Para que la interaccion pp ocurra, la energia del proton relativista debe ser mayor que 
un valor umbral E u ~ m p c 2 + 2m 7r c 2 (l + m 7r /4m p ) ~ 1.22 GeV. La parametrizacion mas 
reciente de la seccion eficaz total de la interaccion pp es la dada por Kelner y colaboradores 
(2006): 

2 

mb, (3.36) 



a pp {E p ) = (34.3 + 1.88 L + 0.25 L 2 ) 



E, 



48 



Procesos radiativos 



CT inel • mb 

: — i i 1 1 1 iii| — i i 1 1 1 1 ii| — i i 1 1 1 ni| — i i 1 1 1 1 1 1 1 — i i 1 1 1 1 1 ■ | — i i 1 1 1 iii| — i i 1 1 1 1 11 




E p , TeV 



Figura 3.3: Section eficaz de la interaction pp (Kelner et al. 2006). 



donde L = ln(E p /l TeV) y su forma se muestra en la Figura [3731 

Como mostramos en la ecuacion (13.7p el tiempo de enfriamiento depende de la section 
eficaz de la interaction, con lo cual para diferentes parametrizaciones de esta tendremos 
diferentes valores del tiempo de enfriamiento por pp. Sin embargo, debido a que todas las 
parametrizaciones halladas son parecidas y que la dependencia con E p es chica, usaremos 
una expresion sencilla pero util para los fines de esta tesis. Considerando que a pp ~ 30 mb 
y que la inelasticidad de la interaction es f pp ~ 0.5, el tiempo de enfriamiento para un 
medio de densidad n m resulta 

2 x 10 15 

tpp ~ s. (3.37) 

n m 

Como ocurre con el proceso Bremsstrahlung relativista, las perdidas por pp no modifican 
la forma del espectro de inyeccion. Si inyectamos protones relativistas con una distribution 
Q P (E P ) oc E~ p , al cabo de un tiempo > t pp tendremos una distribution estacionaria de 
protones tal que 

N p (E p ) ~ Q p (E p ) t pp oc E-p. (3.38) 

Estos protones relativistas, al interactuar con el medio con densidad n m , produciran ir° 
cuya emisividad esta dada por: 

JAE^o) = ± n m / P da ^ E f^ E P> Np{Ep) dEp . (3.39) 

Para E p < 0.1 TeV, los piones creados tendran una energfa promedio E n o ~ K n E p , donde 
n n ~ 0.17 (Gaisser 1990), y por lo tanto es posible considerar la aproximacion 8(E 7r o — 
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K n Ep U ) para calcular J n o. Siendo Ep in = E p — m p c 2 la energfa cinetica de los protones, se 
llega al siguiente resultado: 



J^E^o) ~ n m iVp m p c H a pp I m p c H . (3.40) 

Sin embargo, la aproximacion 5 no es muy buena si E p > 0.1 TeV ya que en este caso la 
distribucion de piones se ensancha en energfa y la parte de mas baja energfa se hace mas 
"blanda" para E p mayores, ya que se crean muchos piones de poca energfa (es decir, la 
multiplicidad crece). En este rango de energfa de los protones, las formulas adecuadas para 
calcular J n o son las dadas por Kelner y colaboradores (2006). De esta manera, usando la 
aproximacion 5 (modificada) para E p < 0.1 TeV y las parametrizaciones de Kelner para 
E p > 0.1 TeV, la multiplicidad no puede considerarse un parametro libre y ajustarlo de tal 
manera que J n o sea una funcion continua. Finalmente, calculada J n o, podemos obtener la 
emisividad de los rayos gamma a traves de 

J ph (E ph ) = 2 f En ° .^° (jM 2 dE w o , (3.41) 

donde E™ n = E p ^ + m^c 4 / '(4_E p h) . El factor 2 en la expresion (I3.4ip tiene en cuenta que 
por cada 7r° se producen dos rayos gamma. 



Creacion de pares electron-positron 

De acuerdo a (13.351) . ademas de 7r° tambien se crean 7r en las interacciones pp. Estos 
ultimos luego decaen en muones que a su vez producen neutrinos y pares electron-positron 
(e 1 * 1 ). La mayorfa de los pares producidos tienen energfas E e2 3> m e c 2 y por esto nos interesa 
conocer la distribucion N e2 (E e2 ) de estas partfculas ya que son una poblacion adicional de 
leptones relativistas que se enfriaran de la misma manera que los electrones primarios 
acelerados, por ejemplo, en frentes de ondas de choque. 

Si N p (E p ) es una ley de potencias, entonces la distribucion de pares e ± creados tambien 

sera una ley de potencias N e2 (E e2 ) = K' e2 E e2 P2 , donde tanto p' 2 como K' e2 dependen de 
los parametros de N P (E P ). Calculando N e2 (E e2 ) traves de las formulas dadas por Kelner y 
colaboradores (2006) y considerando K' e2 = 1 podemos hallar p' 2 de hacer un ejuste grafico 
del espectro calculado. Luego, conociendo el fndice del espectro de los pares podemos hallar 
K' e2 a traves de la igualdad u e , 2 = f e ±u p , donde f e ± depende de p 2 . Las cantidades u e2 y u p 
son las densidades de energfa de los pares e ± y de los protones relativistas, respectivamente. 



3.3. Absorcion 

Anteriormente hemos visto diferentes procesos no termicos que producen fotones. Sin 
embargo, estos pueden a su vez ser absorbidos mediante interacciones con otras partfcu- 
las. Estas interacciones pueden ocurrir con campos de materia, con campos magneticos 
(Sturrock 1971) o con campos de fotones (Coppi & Blandford 1990). 
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Supongamos que un foton de energfa E ph se propaga en un medio de densidad n h y 
contenido en una region de tamaiio L m . Si la emisividad de los fotones es J ph , luego de 
atravesar la region de tamano L m , la emisividad sera reducida a J p h de manera tal que 

J P h(E ph ) = J° h (E ph )e~ T . (3.42) 

La profundidad optica del medio, r, se define de la siguiente manera: 

/•Lm POO 

r(E ph )= / / n h (E h ,r)a ahs (E ph ,E h )dE h dr, (3.43) 

donde <7 a b s es la seccion eficaz del proceso y E™ m es la energia minima de las particulas del 
medio para que ocurra la interaccion entre una de ellas y un foton. Si r 1 entonces no 
hay absorcion y el medio se dice opticamente delgado. Por otro lado, si r > 1 la absorcion 
puede ser significativa y el medio se dice opticamente grueso. 



3.3.1. Absorcion por creacion de pares electron-positron 

En esta tesis solo consideraremos un mecanismo de absorcion; aquel en el cual dos 
fotones se aniquilan produciendo un par e de la siguiente manera 

7 + 7 — >e + + e~. (3.44) 

Si consideramos que el foton de energfa E ph producido en la fuente se propaga en un 
gas isotropico y (quasi)monoenergetico de fotones de energia E® h , para que la interaccion 
ocurra, debido a que se crean dos particulas con masa m e , la energia umbral de los fotones 
debe ser tal que E^E^ > 2m 2 e c l . La seccion eficaz de esta interaccion foton-foton es 
(Coppi & Blandford 1990) 

( x _ i)V3 / l 3 \ 
a 77 (£ ph , E ) ~ a T v 77 5/2 —-^ + - ln(x 77 ) H(x^ - 1) (3.45) 

donde x 77 = E ph E^ h / m 2 e c A y i?(x 77 — 1) es la funcion de Heaviside. Como se muestra en 
la FiguraE3l cr 77 tiene un maximo cuando E ph E® h ~ 3.7?7igC 4 . 

Si la densidad del gas de fotones en el cual se propaga el rayo gamma es n phm , la 
profundidad optica sera: 

r 77 (i?ph, E ph ) = riph cr 77 (E phm ,E° ph )L m . (3.46) 
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Figura 3.4: Section eficaz isotropica (cr 77 ) del proceso de creation de pares e ± por aniqui- 
lacion de fotones. El grafico corresponde a E® h = 1 eV. 
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Parte II 

Estudios de fuentes a diferentes 

escalas 



Capitulo 4 

Objetos estelares jovenes 



4.1. Introduccion 

Las estrellas se clasifican de acuerdo a su masa en de gran masa (M* > 8M Q ) 
y de baja masa (M* < 8M ). Numerosas de estas ultimas estan cerca y son faciles de 
detectar. Por esto se conoce bastante de su formacion y evolution. Sin embargo, no ocurre 
lo mismo con las estrellas de gran masa, las cuales se encuentran embebidas en grandes 
condensaciones de gas y polvo con lo cual la extincion de la luz que emiten es significativa 
y poco llega de ella a nuestros detectores. Si ademas tenemos en cuenta que el tiempo 
de vida de estas estrellas es muy corto (~ 1 — 10 Myr), la observation y detection de las 
estrellas tempranas (tipos espectrales O y B) en cada estado evolutivo es extremadamente 
dificil. Es por esto que el estudio de la formacion de las estrellas de gran masa es uno de 
los grandes topicos de la astrofisica actual. 

La formacion estelar comienza cuando una nube de gas en el espacio se vuelve inestable 
y colapsa bajo la action de su propia gravedad. Durante el colapso, la nube se fragmenta. 
Esto es todo lo que puede decirse respecto de la etapa initial de la formacion de una 
estrella, ya que se sabe muy poco de como ocurre la fragmentation de la nube. Luego, la 
nube queda dividida en varias partes, cada una de ellas con una distribution de densidad 
inhomogenea la cual induce el proceso de acrecion y asi se forman nucleos mas densos. 
Dependiendo de la masa de estos nucleos y de la tasa de acrecion de materia, se formaran 
estrellas de diferentes masas. 

Para los nucleos poco densos que luego daran lugar a las estrellas de baja masa, la 
secuencia de eventos hasta llegar a la formacion de la estrella pareciera ser clara y la 
resumimos de la siguiente manera (Shu et al. 1987): 

1. Un estado initial en el cual la region central de los nucleos densos se contrae, inten- 
sificando asi su campo gravitational y formando una protoestrella en cada nucleo. 

2. Un estado de acrecion caracterizado por la formacion de un disco alrededor de cada 
protoestrella y a traves del cual la misma acreta materia del medio circundante. 

3. La fase en la cual se producen los flujos bipolares {outflows o jets dependiendo del 
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Figura 4.1: Region de formation estelar masiva conocida como Messier 20. Esta imagen ha 
sido tomada en el IR por el satelite Spitzer. 



angulo de colimacion de los mismos), por los cuales la protoestrella deposita materia 
con momento angular y energia cinetica en sus alrededores. Estos flujos bipolares, 
al chocar violentamente con el medio circundante, producen los llamados objetos 
Herbig-Haro (HH). 

4. Finalmente, la etapa en la cual la protoestrella se situa en la secuencia principal 
de edad cero (ZAMS, por Zero Age Main Sequence). Es en este momento cuando 
comienza la combustion eficiente de hidrogeno en el nucleo de la estrella. 

Por otro lado, el mecanismo de formation de las estrellas de gran masa no es tan claro 
hoy en dia. Si bien no hay dudas de que las estrellas tempranas se formarian en nubes 
moleculares gigantes, el proceso de fragmentation de dicha nube para formar grumos auto- 
gravitantes de gran masa es aun desconocido, como asi tambien los subsecuentes procesos 
dinamicos que darfan origen a estas estrellas. Se han planteado dos modelos para explicar 
la formation de las estrellas de gran masa: uno es aquel en el cual el mecanismo de forma- 
tion es similar al que opera en las estrellas de baja masa (Shu et al. 1987), mientras que el 
otro es un modelo tipo jerarquico donde la coalescencia de estrellas menos masivas daria 
origen a estrellas con mas masa (Bonnell et al. 1998). Observaciones recientes de regiones 
de formation estelar gigantes y masivas (ver la Figura 14.11) han detectado la presencia de 
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discos de acrecion y de flujos bipolares emanando de protoestrellas de gran masa, como 
describiremos a continuation. 

Cumulos: Observaciones en radio de regiones de formation estelar masivas de la Ga- 
laxia y de galaxias cercanas muestran que las estrellas de gran masa se forman en grupos. 
Las fuentes IRA 150 mas luminosas asociadas a regiones compactas de gas ionizado (y por 
ende a objetos estelares jovenes), muestran una morfologia compleja cuando se las observa 
en frecuencias radio. Garay y colaboradores (1993) sugieren que la presencia de estrellas 
embebidas en la region (un cumulo de estrellas O y B) excitarfan el gas de la misma y 
esto producin'a la compleja estructura observada. Algunas de las regiones mas estudiadas, 
que contienen estrellas masivas, son W3 Orion- Trapezium y NGC 2024. La densidad tipica 
de estos cumulos de estrellas jovenes es de ~ 10 4 estrellas por pc~ 3 y el tama no es de 
~ 0.2 - 0.4 pc. 

Discos de acrecion: Los discos detectados en regiones de formation de estrellas de 
gran masa tienen un diametro entre 0.1 y 1 pc, y sus masas pueden ir desde las 10 a las 
2000 M . Una de las evidencias observacionales mas importantes de la existencia de discos 
circumestelares en regiones de formation estelar masiva es la detection de la emision de 
gas y polvo en el infrarrojo (IR). 

Jets y outflows: La presencia de flujos bipolares en la formation de estrellas de 
gran masa fue determinada recientemente mediante observaciones en radio. Sobre una 
muestra de aproximadamente 120 regiones de formation de estrellas tempranas, en el 90 % 
de los casos se observo gas moviendose a velocidades altas (> 100 km s" 1 ). Luego, si las 
velocidades elevadas del gas estan asociadas a la presencia de outflows, puede concluirse 
que estos son comunes tambien en la formation de las estrellas de gran masa. Algunas de 
las propiedades de los outflows detectados en regiones de formation de estrellas tempranas 
son las siguientes: 

■ Tienen una masa promedio de 130 M , aunque en algunos caso puede alcanzar valores 
de hasta ~ 4800 M . 

■ La tasa de perdida de masa abarca desde ~ 3 x 1O _5 M yr -1 hasta ~ 3 x 1O~ 2 M 
yr -1 . 

■ La energia cinetica toma valores desde ~ 10 46 hasta ~ 6 x 10 48 erg. 

Estas cantidades son ~ 100 veces mayores que las correspondientes a los outflows presentes 
durante la formation de estrellas de baja masa. Luego puede concluirse que los outflows de 
alta masa inyectan mayor cantidad de energia al medio circundante que los de baja masa. 
En la Figura 14721 se muestra el outflow detectado en la fuente AFGL 2591. 

Los resultados observacionales descriptos anteriormente sugieren que las estrellas de 
gran masa se forman a traves de un mecanismo similar al que opera para formar las 
estrellas de baja masa, pero con una tasa de acrecion de materia a traves de los discos y 
de eyeccion a traves de los outflows ~ 100 veces mayor a la observada en estrellas de baja 
masa. Sin embargo, no se han detectado discos de acrecion en estrellas de > 50 M , 



Se conocen como fuentes IRAS a aquellas detectadas por el satelite IR que lleva dicho nombre. 
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Gemini Obsefvatory'Colin A spin 



Figura 4.2: Outflow masivo asociado a la fuente AFGL 2591. La imagen ha sido tomada 
en el IR por el telescopio Gemini Norte (Zinnecker & Yorke, 2007). 



4.2 La fuente IRAS 16547-4247 
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las cuales al ser tan luminosas (L* > 10 5 L Q ) podrian evaporar los discos a traves de 
un mecanismo conocido como "fotoerosion" . Zinnecker & Yorke (2007) sugieren que el 
mecanismo para formar a las estrellas mas masivas (M* > 50 M Q ) no seria un simple 
escaleo del proceso de formacion de las de baja masa. Los fuertes vientos de las estrellas 
mas masivas (junto a la fuerte radiacion) disiparian el disco y favorecerfan las interacciones 
entre protoestrellas cercanas, siendo entonces la coalescencia de estrellas o protoestrellas 
el mecanismo de formacion. 

Por todo lo expuesto anteriormente, es claro que todavia hay mucho para decir respecto 
de como se forman las estrellas de gran masa. Todos los estudios y observaciones de regiones 
de formacion estelar masivas son en general en las bandas de radio, IR y rayos X, mientras 
que a mas altas energias solo se han observado algunos cumulos y asociaciones de estrellas 
O y B. Hasta el momento, no se han detectado estrellas en formacion en el rango de los 
rayos gamma, aunque posiblemente se han detectado rayos gamma de estrellas de gran 
masa evolucionadas, en general Wolf-Rayet (WR) o estrellas luminosas variables azules 
(LBV, por Luminous Blue Variable). Las regiones en las que se ha detectado emision 
significativa en rayos gamma son Westerlund 2 (Aharonian et al. 2007), Cygnus OB2 y 
Carina (Tavani et al. 2009a). Las particulas relativistas que emiten esta radiacion gamma 
detectada podrian acelerarse en las ondas de choque producidas por las explosiones de 
supernovas dentro de la region o bien en choques de vientos de estrellas de gran masa. Sin 
embargo, estos cumulos o asociaciones de estrellas tempranas tambien pueden ser regiones 
de formacion estelar, aunque nada se habia dicho de la posible emision a altas energias de 
las protoestrellas masivas. Los rayos gamma, de ser producidos en regiones de formacion 
estelar, aportarfan information adicional a la que se obtiene de otras frecuencias sobre 
como es el mecanismo de formacion de las estrellas de gran masa, motivando esto el primer 
estudio realizado en esta tesis. 

Sabemos que para emitir rayos gamma es necesaria la presencia de particulas no termi- 
cas. Una de las evidencias de la existencia de electrones relativistas es la detection en 
frecuencias radio de emision no termica. Se han detectado en radio muchos YSOs de gran 
masa, pero solo algunos con emision no termica. Para nuestro estudio hemos elegido una 
fuente muy particular detectada en radio (no termico) por Garay y colaboradores (2003). 
Estos autores han detectado en frecuencias radio un sistema triple compuesto por una 
protoestrella de gran masa y dos lobulos. Este sistema se asocia con la fuente infrarroja 
IRAS 16547-4247. 

4.2. La fuente IRAS 16547-4247 

La fuente IRAS 16547-4247 es una region de formacion estelar de gran masa ubicada 
a una distancia d mn ~ 2.9 kpc y tiene un diametro angular de ~ 27" (Garay et al. 2003), 
que a la distancia d nm corresponde a un tamaho lineal D mn ~ 0.38 pc (~ 1.1 x 10 18 cm). 
La luminosidad bolometrica es ~ 6.2 x 10 4 L Q ~ 2.4 x 10 38 erg s _1 donde L Q = 
4x 10 33 erg s -1 es la luminosidad bolometrica del Sol. De la densidad columnar de hidrogeno 
estimada, resulta que la masa de la nube es M nm ~ 9 x 10 2 M© ~ 10 36 gr, siendo M ~ 
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2 x 10 33 gr una masa solar, y asi la densidad numerica promedio de particulas resulta 
n nm ~ 5.2 x 10 5 cm -3 . Los parametros conocidos y estimados de esta fuente estan listados 
en la Tabla I4TT1 



Parametro 


Valor 




Distancia 


4m — 


-- 2.9 kpc 


Tamaiio 




= 0.38 pc 


Masa 


M nm 


= 9 x 10 2 M 


Densidad 




= 5.2 x 10 5 cm" 3 


Luminosidad 


L* = 


6.2 x 10 4 L© 



Tabla 4.1: Parametros de la fuente IRAS 16547-4247. 

Mediante observaciones llevadas a cabo con el interferometro ATCA [Australia Telesco- 
pe Compact Array), Garay y colaboradores (2003) detectaron un sistema triple embebido 
en la nube molecular asociada a la fuente IRAS 16547-4247. Este sistema esta formado 
por una protoestrella de gran masa y dos lobulos alineados con esta y ubicados simetrica- 
mente a una distancia ~ 0.14 pc (~ 10 ) de la protoestrella central, como se muestra en 
la Figura 14.31 Con observaciones posteriores realizadas con el radiotelescopio VLA ( Very 
Large Array) por Rodriguez y colaboradores (2005) a las frecuencias v = 8.46 y 14.9 GHz, 
se calculo el mdice espectral a del ffujo observado S u (oc v~ a ) de cada componente del 
sistema triple. Para la fuente central se estimo a = —0.33 ± 0.05 mientras que para los 
lobulos norte y sur, a = 0.17 ± 0.39 y 0.59 ± 0.15, respectivamente. Estos valores son muy 
similares a los calculados por Garay y colaboradores (2003) con los datos de ATCA. 

Si la emision continua detectada en radio es producida unicamente a traves de inter- 
acciones libre-libre, entonces el espectro observado debe corresponder a un valor a < 0.1, 
independientemente de las caracterfsticas del emisor. Por otro lado, valores a ^> 0.1 son 
debidos a radiacion sincrotron opticamente delgada (Rodriguez et al. 1993). Con este cri- 
terio, la emision radio de la fuente central seria termica, asociada a la radiacion libre-libre 
de la base del jet y que se conoce como jet termico, mientras que la emision proveniente 
de los lobulos, en particular del sur, seria no termica. Sin embargo, para confirmar que la 
emision de los lobulos es no termica, deben hacerse estudios polarimetricos. 

La deteccion de emision no termica proveniente de los lobulos es una evidencia de la 
presencia de electrones relativistas en tal lugar (donde esta radiacion es producida). Estas 
particulas serian aceleradas hasta energias relativistas en los choques producidos cuando 
el jet es frenado por el medio en el cual se esta propagando, es decir, la nube molecular. 
En esta tesis modelamos a la fuente como se muestra en la Figura 14.41 

Ademas de radiacion sincrotron, los electrones relativistas podrian producir fotones a 
mas altas energias por otros mecanismos radiativos como el Bremsstrahlung relativista y la 
dispersion Compton inversa. Por otro lado, debido a que el mecanismo de Fermi puede ace- 
lerar tambien protones, las interacciones pp entre protones acelerados y protones frios de la 
nube son tambien posibles. Como consecuencia de esto se producen, ademas de los fotones, 
pares e^ secundarios los cuales se enfrfan de la misma manera que los electrones primarios. 
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Figura 4.3: Imageries tomadas con VLA de la fuente IRAS 16547-4247 en v = 8.46 (izquier- 
da) y 14.9 GHz (derecha). A estas frecuencias se ve claramente el sistema triple compuesto 
por la fuente central (central source) y los lobulos norte (N) y sur (S) (Rodriguez et al. 
2005). 



De esta menera, podriamos tener tres poblaciones diferentes de particulas no termicas: 
electrones y protones primarios y pares secundarios e . A continuation estudiamos como 
se producen estas poblaciones de particulas relativistas. 



4.3. Poblacion de particulas relativistas 

La emision detectada en frecuencias radio proveniente de los lobulos de la fuente 
IRAS 16547-4247 es no termica (de acuerdo al criterio determinado por Rodriguez et al. 
1993) y consideramos que dicha radiation es producida por electrones relativistas acelera- 
dos en los choques teminales de los jets del YSO o por pares e ± producidos por protones 
relativistas acelerados en los mismos choques. Mientras los jets se propagan (libremente) 
por la nube molecular van chocando el material de esta produciendo asi un bow shock en 
el medio. Cuando los jets son frenados significativamente por este medio externo chocado, 
se forma un choque reverso en los jets, como se comento en la section 12.5.21 A continua- 
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Figura 4.4: Esquema del modelo propuesto para la fuente IRAS 16547-4247 (Araudo et al. 
2007). 



cion haremos una description mas detellada que la expuesta en el Capitulo El estudiando 
especificamente el caso de los YSOs, para luego estimar la poblacion de particulas rela- 
tivistas en los lobulos de la fuente IRAS 16547-4247. En particular, nos concentramos en 
el estudio del lobulo sur, ya que la emision radio proveniente de este tiene un indice es- 
pectral puramente no termico, a ~ 0.6, mientras que la emision del lobulo norte podria 
estar contaminada de radiation termica, ya que el a (~ 0.2) es mas cercano al valor critico 
a c = 0.1. 

4.3.1. Choques terminales 

Para determinar las caracterfsticas de los choques terminales, un parametro importante 
es el contraste entre la densidad del jet y de la nube, x y so = n i/ n nm- Ademas de x yso es 
necesario conocer la velocidad del jet, v- v para poder estimar las velocidades del bow shock 
(t>bs) y del choque reverso (t> r ), de la siguiente manera (Bosch- Ramon et al. 2010): 

Vhs ~7i — ^17^7 y v T ~v i --v b „ (4.1) 

(1 + Xyso ) 4 

donde se ha supuesto que el bow shock es adiabatico. Debido a que la velocidad de los jets 
de la fuente IRAS 16547-4247 no se conoce, adoptamos para esta magnitud un valor que 
sea concordante con aquellos medidos en otros YSOs de gran masa. Marti y colaboradores 
(1995) han determinado para el jet termico de la fuente HH 80-81 una velocidad que ronda 
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los 600 y 1400 km s -1 . Luego, para la fuente que estamos estudiando en este capftulo 
consideramos Vj = 1000 km s _1 . 

Si bien no conocemos el valor de Xyso, porque no conocemos nj en la localization de los 
lobulos (Zj ~ 4 x 10 17 cm), podemos determinar t>bs a traves del tiempo de vida del jet, 
Tj. Esta escala temporal se define como aquella en la cual el jet recorre la distancia Z- } a la 
que se detectan los lobulos. Esto es, 73 ~ Zj/t>b s y en el caso de la fuente IRAS 16547- 
4247, Tj ~ 10 11 s (Garay et al. 2007) con lo cual, si Zj ~ 4.3 x 10 17 cm obtenemos 
v bs ~ 5 x 10 6 cm s _1 . Con este valor de t>b s y las ecuaciones (14. ip hallamos x y so ~ 2.5 x 10~ 3 
y luego rij(Zj) ~ 1.25 x 10 3 cm -3 . Siendo t>b s «C v- v resulta que v T ~ v r 

La region entre ambos choques puede tener un estructura muy compleja ya que estos 
choques pueden ser radiativos y en la superficie de discontinuidad entre ambos medios 
chocados (conocida como "superficie de trabajo") se pueden desarrollar inestabilidades de 
RT que mezclan material de ambos medios (Blondin et al. 1989). Para estimar si los choques 
son radiativos o no, comparamos el tiempo de enfriamiento por emision termica, t te r, con 
el tiempo caracterfstico t carac en el cual el medio chocado recorre una distancia ~ -Riob- 
Este ultimo resulta t 

carac ~ 4 -Riob/^ch ■ Para estimar tt er consideramos la ecuacion (I2.55P y 
obtenemos (Bosch- Ramon et al. 2010): 

t rad ~ 2 x 10- 25 ^— ~ 7 x 10- 18 ^ s, (4.2) 
n A(T) n 

donde hemos considerado la funcion de enfriamiento A(T) = 7x 10~ 19 T -06 (Bosch- Ramon 
et al. 2010) y las condiciones de salto (12.271) y ( 12. 28ft . Luego, comparando t ra( j con t C arac 
tenemos que los choques son radiativos si t ra d < t C arac, es decir, si 

A r> \ 1/4.2 

4 itiob Uq \ 1 a ^ g x/4.2 _i 



Vch < [tVW^) - L6 >< l0 ^o cms-. (4.3) 

Considerando v ch = t> bs y n = n nm para el bow shock y t> ch = v r y n = rij(Zj) para 
el choque reverso, obtenemos que mientras que el primero resulta radiativo, el ultimo es 
adiabatico. Luego la densidad del material chocado de la nube aumenta mas de un factor 
4, pudiendo llegar hasta valores mucho mas grandes que el valor adiabatico. Con esto, 
si el factor de mezcla de ambos medios chocados es grande, podemos considerar que las 
particulas aceleradas en el choque reverso pueden interactuar, ademas de con el material 
del jet chocado, con material chocado de la nube, cuya densidad es mucho mas alta y asi las 
perdidas por Bremsstrahlung relativista e interacciones pp mucho mas eficientes. 

En resumen, tenemos que los lobulos son regiones en donde hay choques fuertes que 
pueden acelerar particulas hasta energfas relativistas. Estas particulas, ademas de pro- 
ducir la radiation sincrotron observada, pueden tambien radiar por otros mecanismos no 
termicos, como veremos a continuation. 

4.3.2. Aceleracion de particulas y perdidas radiativas 

Como vimos en el Capftulo [2, la eficiencia para acelerar particulas mediante el mecanis- 
mo de Fermi de tipo I es oc w 2 h . Por esto, nos concentramos en la aceleracion de particulas 
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en el choque reverso, ya que v x ^> v^ s . 

El tiempo de aceleracion depende, ademas de v 2 , del campo magnetico en la region 
de aceleracion, que coincide con el lobulo sur, y denotaremos B\ oh . Para estimar £> lob 
consideramos que la densidad de energia magnetica, Bf oh /(87i), es igual a la densidad 
de energia u nt de las particulas no termicas. Sabemos que hay leptones relativistas en el 
lobulo que producen la emision sincrotron detectada en frecuencias radio. Estos leptones 
pueden ser electrones primarios (ei) acelerados en el choque reverso o bien pares e ± (e 2 ) 
producidos por interacciones inelasticas pp de protones (p) acelerados en el mismo choque 
con el material de la nube molecular. El caso mas general que podemos plantear es 

B 2 

-J2k = Unt = u + u + u (4.4) 

87T 

donde u ei y u p son las densidades de energia de los electrones y protones primarios, res- 
pectivamente, y u e2 corresponde a los pares e ± . En cada caso {i = ei,p,e 2 ). 

Ui = I Ein(Ei)dEi, (4.5) 



siendo n(Ej) = k' i E i Pi exp(— E'j/ i^ max ) la distribucion de energia de las particulas i por 
unidad de volumen, es decir, [n(.E'j)] = cm -3 erg -1 y Ef 1 ™ y i^ max las energias minima 
y maxima de la distribucion n(Ei). Como no tenemos ningun indicio observacional de la 
presencia de protones relativistas en el lobulo sur, la unica manera de estimar u v es a traves 
de un parametro fenomenologico a, imponiendo que u v = au ei . En esta tesis consideramos 
3 casos: a = 0, 1 y 100. El valor a = 100 se toma porque el espectro de los rayos cosmicos 
galacticos pareciera indicar que la fraccion de protones a electrones es ~ 100. Por otro lado 
u e2 = fe±u p , donde f e ± < 1 (Kelner et al. 2006). 

El espectro observado en radio del lobulo sur puede ajustarse con una ley de potencias 
(oc v~ a ) con indice a = 0.59. Luego, el indice espectral de los leptones relativistas que 
producen esta radiacion sincrotron resulta p' e = 2a + 1 = 2.18. En los casos a = 1 y a = 
100 consideramos ademas una distribucion n p (E p ) de protones relativistas cuya poblacion 
de pares e ± producidos en las interacciones pp puede estimarse como se explico en el 
Capitulo |3j En el caso con a = 1, u e2 = f e ± u ei con lo cual la emision de los primarios es 
tambien mayor que la de los secundarios y ajustamos el espectro observado con la emision 
de los electrones acelerados en el choque reverso. Esto es, j/ = 2.18 y fijando p' = j/ 
resulta p' e2 = 2.13. Por otro lado, en el caso con a = 100, resulta u e2 ~ 50u ei y asi el 
espectro observado se ajusta con la emision sincrotron producida por los pares. De esta 
manera, si p' e2 = 2.18, entonces p' = 2.27 y por lo tanto p' ei = 2.27 tambien. 

Fijados los valores de p[ para cada caso y las relaciones entre las densidades de energia Ui, 
podemos determinar las constantes k[. Para esto, nos valemos del flujo sincrotron observado 
en radio, cuya expresion es la siguiente (Ginzburg & Syrovatskii 1964): 

^ = 1.35xlO' 22 a sin ( P :) k >\y r 26Xl ° ) , (4.6) 

d nm V V J 
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a 


-Biob 
[G] 


V 

[erg Pe i _ cm~ 3 ] 


k' p 

[erg p p _1 cm -3 ] 


K 

[erg Pe 2 _1 cm -3 ] 





2.0 x 10" 3 


2.6 x 10" 9 






1 


2.5 x 10-3 


2.0 x lO" 9 


6.7 x lO" 9 


6.5 x 10~ 10 


100 


3.0 x 10~ 3 


2.7 x 10~ n 


1.7 x 10~ 8 


1.5 x 10~ 9 



Tabla 4.2: Campo magnetico y constantes de normalizacion (k^) para los diferentes casos 
considerados en esta tesis. 



donde a sin (p-) es una funcion complicada tabulada para diferentes valores de p[. Para 
p[ = 2.18 resulta a sin ~ 0.1, interpolando linealmente entre a s i n (2) y a S i n (2.5). El volumen 
del lobulo sur es Vj Q b ~ 5.6 x 10 48 cm 3 y d nm es la distancia a la fuente. Considerando que 
S(v = 8.46 GHz) = 2.8 mJy (Garay et al. 2003) y la ecuacion (I4.4p de equiparticion de 
la energia es posible obtener B\ oh y k\ {% = e\ 6 i = e 2 ) para cada valor del parametro 
a. Luego, teniendo en cuenta las relaciones entre las densidades de energia (u p = au ei y 
u e2 ~ f e ±u p)i hallamos k' y la constante leptonica que falta. En la Tabla 14721 se listan los 
valores de -Bi b, k' ei , k' p y k' e2 calculados para los 3 valores de a considerados. 

Conociendo el valor de _Bi b podemos estimar el tiempo de aceleracion tanto para proto- 
nes como para electrones primarios. Aunque no tenemos information sobre la geometria de 
-Blob en los lobulos, aqui adoptamos el valor mas conservative para la eficiencia acelerativa 
(-Biob || #ch) y es por esto que consideramos la formula (I2.59j) para calcular el tiempo de 
aceleracion, resultando 

t ac = 1.4 x 10 10 ~ 1-4 x 10- 6 §^ s, (4.7) 

-Dlob V r Mob 

si consideramos que v r ~ uj ~ 1000 km s _1 . 

Como vimos en el Capftulo [2], la aceleracion continua hasta que las particulas o se 
escapan del acelerador o bien se enfrian por perdidas radiativas. En el primer caso, el 
tiempo de escape t esc esta dado por la ecuacion (12.631) . siendo el tiempo de convection 

tco„v ~~4xl0 8 s, (4.8) 

Vbs/4: 

y el tiempo de difusion en el regimen de Bohm resulta 

t dif ~ 4.8 X 10-20 ^lob ^lob _ 5>g x 1Q 6^ob g _ (4 g) 

Ademas de acelerarse y/o escaparse, las particulas primarias y tambien los pares secun- 
darios pueden perder energia a traves de interacciones sucesivas con otras particulas o 
camp os. 
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Perdidas leptonicas 

Como vimos en el Capftulo [21 las interacciones de leptones relativistas con campos 
magneticos producen radiacion sincrotron, cuyas perdidas resultan (ver ecuacion (13.151) ): 

t sin ~4x 10 2 -2-^ s. (4.10) 

-"lob ^ei,e 2 

Debido a que -Bi OD cambia con a, para cada valor de este parametro tenemos un valor 
diferente de t 8 i n , al igual que de t ac . 

Los leptones tambien pueden interactuar con los fotones IR producidos por la proto- 
estrella masiva. La densidad de energia de estos fotones es w p h = L*/(47r B^ m c), donde 
L* es la luminosidad de la protoestrella y R nm = D nra /2 es el radio de la nube. Resulta 
M p h ~ 1.84 x 10~ 9 erg cm -3 y los fotones correspondientes a esta densidad tienen una 
energia E p ^ ~ 3KbT uui ~ 6.6 x 1CT 3 eV al estar reprocesados por el polvo. Con fotones se- 
milla de tan baja energia, las interacciones IC ocurriran en el regimen de KN para energfas 
tan altas como E ei>e2 > {m e c 2 ) 2 / ~ 3.8 x 10 13 eV. Tanto en el regimen de Th como en 
KN la fomula valida para el tiempo de enfriamiento es la ecuacion (13.231) . que en el caso 
del lobulo sur nos queda 

tlo ~3.3xl0"£ <I±Ml> (l±iVL_ s, (4.11) 
C p ln(l + 0.2?/) (1 + 0.5 ?/ + 1.3 ?/ 2 ) ' V ; 

donde y = E^E eu J(B.l x 10 5 eV) 2 . 

Por otro lado, debido a que la nube tiene una densidad alta, n nm ~ 5 x 10 5 cm -3 , las 
perdidas por Bremsstrahlung relativista seran importantes, con un tiempo caracteristico 
(ver ecuacion (I3.32p ) 

2.8 x 10 10 fA , 

^Brem ~ n B. (4.12) 

In (f^L) +0.36 

Para obtener la energia maxima que pueden alcanzar los electrones acelerados en el 
lobulo sur, comparamos las ganancias por aceleracion con las perdidas radiativas. Para 
esto hemos graficado los tiempos defmidos antes como se muestra en la Figura l4T5l (izquier- 
da). Como vemos en el grafico de la izquierda, los electrones alcanzan energfas maximas 
~ 1 TeV, para el caso a = 1, como consecuencia del balance entre la aceleracion y las perdi- 
das por radiacion sincrotron. Como puede apreciarse tambien en la Figura 1431 (izquierda) . 
para E e < 3 x 10 10 eV, las perdidas por sincrotron dejan de dominar para hacerse mas 
importantes aquellas por Bremsstrahlung relativista. A esta energia de quiebre, que 11a- 
mamos E q , el espectro de electrones primarios sufre un cambio en su dependencia con E e , 
ya que las perdidas dominantes cambian. El valor de E q puede determinarse igualando 
£ S in = ^Brem, con lo cual obtenemos E q ~ 1.8 x 10 5 /Bf oh eV. En la Tabla I4~41 listamos los 
valores de _E™ ax y E q para los tres valores de a considerados. 

Dado que el tiempo de vida de la fuente r vi d a (~ ij ~ 10 11 s) es mucho mayor que los 
tiempos de enfriamiento dominantes (sincrotron y Bremsstrahlung relativista), el espectro 
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log(E 5i [eV]) log(E p [eV]) 

Figura 4.5: Tiempos de aceleracion y de enfriamiento para los electrones primarios (izquier- 
da) y para los protones (derecha). Estos graficos corresponden al caso con a — 1, pero al 
resto de los casos les corresponden figuras similares, ya que la unica diferencia esta en el 
valor de -Biob, que afecta a £ ac , t S i n 

y tat, y este no cambia demasiado, como se muestra en 

la TablalU 



de energfa de los electrones relativistas inyectados en el lobulo se encuentra en el estado 
estacionario. Es decir, toda la energfa de estas particulas relativistas se radfa, basicamente, 
por emision sincrotron y Bremsstrahlung relativista dentro de la fuente (el lobulo sur), 
ya que los tiempos radiativos son similares al tiempo de escape (t esc ~ 4 x 10 8 s). Para 
mostrar esto hemos calculado el espectro de particulas resultante de la inyeccion q ei (E ei ), 
a diferentes tiempos de inyeccion r iny . Considerando r iny = 10 7 , 10 8 , 10 9 y 10 10 s obtuvimos 
los resultados que se muestran en la Figura 14.61 Finalmente, el espectro de electrones 
relativistas en el lobulo sur resulta 

„ (p \ _ / Qe(E e )t Brem si E e < E q 

ne{Ee) -\ qe (E e )t sm si E e >E q . ^ 

Perdidas hadronicas 

En el caso de los protones, estos pierden energfa principalmente por interacciones con 
los protones frfos de la nube molecular, cuya densidad es n nm y el tiempo de enfriamiento 
por pp resulta: 

tpp ~ 4.3 x 10 9 s. (4.14) 

Como se muestra en la Figura l4T5l (derecha) . las perdidas por difusion son importantes solo a 
energfas muy altas, cercanas a E™ x . Esta ultima queda determinada entonces comparando 
el tiempo de aceleracion con el de difusion, resultando E™ ax ~ 1.3 x 10 18 Bf oh eV, como se 
muestra en la Tabla [4741 
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Figura 4.6: Distribuciones de energfa de los electrones primarios. En cada panel mostramos 
las distribuciones para diferentes tiempos de inyeccion de partfculas (T in j). Las diferentes 
curvas corresponden a los siguientes valores de T m y. 10 7 s (lmea de rayas largas, verde), 
10 8 s ( lmea de rayas y puntos, azul), 10 9 s (lmea punteada, rojo) and 10 10 s (lmea solida, 
negro). Notar que para edades > 10 9 s el estado estacionario es alcanzado (Araudo et al. 
2007). 



a 


E q 
[eV] 


V 

, ei,q 

[erg Pe i cm -3 ] 


V 

, 2,q 
[erg Pe 2 cm -3 ] 





9.4 x 10 9 


3.7 x 10- 11 




1 


3.1 x 10 10 


9.8 x 10" 11 


3.2 x 10" 11 


100 


2.2 x 10 10 


9.5 x 10- 13 


5.4 x 10" 11 



Tabla 4.3: Energfa de quiebre del espectro de electrones para los diferentes casos discutidos 
en esta tesis. Las constantes de normalization k[ [i = ei, e 2 ) de las distribuciones de energfa 
de los leptones para E eue2 > E q estan tambien listadas. 



Considerando que el tiempo mas relevante es t pp , cuya dependencia con E p es despre- 
ciable (ver la Figura I4TB1) . el espectro de los protones relativistas resulta estacionario y con 
una dependencia con E p similar a la de la inyeccion, es decir 

n P {E p ) ~ q p (E p ) t pp cx E~ Vv . (4.15) 

Las interacciones pp dan lugar a la produccion de pares e ± , los cuales radian luego 
por los mismos procesos leptonicos descriptos en la seccion anterior. Como en el caso de 
los electrones primarios, hemos calculado la evolucion temporal del espectro de los pares 
obteniendo que para r vi d a ~ 10 11 s el espectro se encuentra en el estado estacionario. En la 
Figura 14.71 se muestran los resultados hallados. 
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a 


[eV] 


j£ max 

[eV] 


rp max 

[eV] 





3.1 x 10 12 






1 


3.1 x 10 12 


4.7 x 10 13 


1.1 x 10 12 


100 


2.8 x 10 12 


5.7 x 10 13 


1.8 x 10 12 



Tabla 4.4: Energias maximas obtenidas para los electrones y protones acelerados en el 
lobulo sur. Ademas, las energias maximas de los pares secundarios se muestran en la ultima 
columna. 



-12.5 
-13 



secondary pairs 
a=1 




secondary pairs 
a=100 




Figura 4.7: Distribuciones de energia de los pares secundarios. En cada panel mostramos 
las distribuciones para diferentes tiempos de inyeccion de particulas (T iny ). Las diferentes 
curvas corresponden a los siguientes valores de r iny : 10 7 s (lmea de rayas largas, verde), 10 8 
s ( linea de rayas y puntos, azul), 10 9 s (linea punteada, rojo) y 10 10 s (lmea solida, negro). 
Notar que para edades > 10 9 s el estado estacionario es alcanzado. La energia minima de 
los pares inyectados en el lobulo sur es ~ 5 x 10 7 eV, como se aprecia en las figuras 
(Araudo et al. 2007). 
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4.4. Distribuciones espectrales de energfa 

En las secciones previas vimos que en los choques terminales (en particular el choque 
reverso, en el caso de la fuente que estamos estudiando) se pueden acelerar partfculas hasta 
energfas > 10 13 eV y luego estas se enfrfan eficientemente por diversos procesos radiativos 
no termicos. En esta seccion calculamos las SEDs producidas por estos procesos radiati- 
vos, teniendo en cuenta las distribuciones rii(Ei) estimadas anteriormente y los campos 
magneticos, de materia y de fotones presentes en el lobulo sur de la fuente IRAS 16547- 



Asumiendo que el emisor es homogeneo y esferico con un radio i?i Q b y volumen Viob, 
consideraremos en todos los casos que la luminosidad especifica es 



donde j p h(E ph , Ei) es la emisividad de cada proceso radiative 

4.4.1. Interacciones leptonicas 

Considerando los valores de mostrados en la Tabla 14.21 calculamos la luminosidad 
especifica emitida en el lobulo sur por radiacion sincrotron. Por otro lado, con los valores 
de M p h y n nm calculamos las emisividades de los procesos IC y Bremsstrahlung relativista 
usando las formulas descriptas en el Capftulo 3. En la Figura 14.81 se muestra la SED 
correspondiente al caso puramente leptonico, a = 0. 

Como se muestra en la figura antes mencionada, la radiacion a energfas altas esta domi- 
nada por la emision producida por el mecanismo Bremsstrahlung relativista, con un pico 
de ~ 10 32 erg s" 1 a E pil ~ 1 MeV. A energfas E ph > 1 GeV, la fuente presenta lumino- 
sidades de ~ 10 31 erg s _1 , con un cut-off a E ph ~ 10 GeV. En rayos X, en el rango 1-10 
keV, las luminosidades esperadas son de ~ 10 30 erg s _1 , presentando un ablandamiento del 
espectro debido al cut-off exponencial exp(— E ei jEf^) de la distribution de los electrones 
primarios. En las energfas mas bajas, el espectro calculado (sincrotron) ajusta muy bien 
los puntos observados a las frecuencias v = 8.4 y 14.9 GHz. 

4.4.2. Interacciones hadronicas 

En los casos con a = 1 y a = 100 la emision de los protones es significativa, ya que 
la densidad n nm de la nube molecular es alta. Calculando la emisividad pp a traves de la 
formula f !3.39|) obtenemos las luminosidades espeefficas que se muestran en la Figura H~9| 
ademas de las leptonicas (sincrotron, IC y Bremsstrahlung relativista). En el caso con 
a = 1 la luminosidad emitida por pp es similar a la emitida por Bremsstrahlung relativista 
(~ 10 32 erg s _1 ) pero la primera se extiende hasta energfas mas altas, ~ 1 TeV. Los leptones 
primarios dominan la emision por sincrotron y por Bremsstrahlung relativista en el caso 
a = 1, mientras que en el caso con a = 100 la contribucion mas importante a la SED 
es producida por los pares secundarios. En ningun caso la emision por IC es significativa. 



4247. 




(4.16) 
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Figura 4.8: Distribution espectral de energia para el caso puramente leptonico, a = 0. La 
lmea negra representa la suma de las tres contribuciones (sincrotron, IC y Bremsstrahlung 
relativista) a la SED total. En las frecuencias v = 8.4 y 14.9 GHz (-Eph ~ 3.5 x 10~ 5 y 
6.2 x 10~ 5 eV, respectivamente) se representa con cuadrados rojos la luminosidad observada 
con VLA a estas frecuencias. En el rango de los rayos X se muestra la sensibilidad del 
satelite XMM, mostrando que nuestra prediction teorica esta de acuerdo con el hecho de 
que este instrumento no haya detectado a la fuente (Araudo et al. 2007). 



Basicamente la emision para > 1 GeV es debida al decaimiento de los 7r°, mientras 
que la emision en rayos X y en rayos gamma blandos es producida por Bremsstrahlung 
relativista. Para los casos a = 1 y a = 100, las componentes espectrales debidas a la 
radiation sincrotron y Bremsstrahlung relativista de los pares secundarios tienen la misma 
forma que la emision de los primarios (para los mismos casos) pero diferentes a bajas 
energias debido a que la energfa minima de cada distribution de particulas (electrones 
primarios y pares secundarios) es distinta. En el caso de los electrones primarios la energia 
minima considerada es ~ 2m e c 2 mientras que los pares tienen una energia minima 
determinada por la energia umbral para la creation de los piones cargados: E n ~ 1.22 GeV. 

4.5. Discusion 

Las fuentes de rayos gamma del piano galactico son usualmente asociadas con regiones 
de formation estelar (Romero et al. 1999). Sin embargo las fuentes detectadas por egret| 
son mucho mas luminosas que los YSOs de gran masa, de acuerdo a los calculos realizados 
con nuestro modelo. Las fuentes EGRET no identificadas tendrian como contrapartidas 
a piilsares, remanentes de supernovas, estrellas tempranas (Romero 2001) y MQs (Bosch- 
Ramon et al. 2005). La emision producida por YSOs de gran masa esta por debajo de 

2 E1 instrumento EGRET (a bordo del satelite Compton) funciono en los anos 1991-2000 y detectaba 
fotones con energias 30 MeV < E p h < 300 GeV. Con los datos obtenidos se confecciono el primer catalogo 
de fuentes de rayos gamma que permitio realizar estudios poblacionales. 
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Figura 4.9: Distribuciones espectrales de energia para los casos con a = 1 (izquierda) y 
100 (derecha), en los cuales consideramos tanto una poblacion leptonica como hadronica 
de partfculas primarias y la poblacion de pares producida. Con (1) y (2) indicamos la 
contribution de los leptones primarios y secundarios, respectivamente. Como menciona- 
mos en la Figura [OJ la lmea negra indica la SED total (suma de todas las contribuciones, 
leptonicas y hadronica) y los cuadrados rojos indican la luminosidad observada con VLA. 
La sensibilidad de XMM es tambien graficada mostrando nuevamente que los lmiites ob- 
servacionales en rayos X tampoco se violan en los casos con a = 1 y 100 (Araudo et al. 
2007). 



la sensibilidad de EGRET, pero estos podrian ser detectados por Fermi. En el rango de 
las VHE, los futuros arreglos de telescopios Cherenkov como HESS II, MAGIC II y CTA^I 
podrian detectar fuentes como IRAS 16547-4247. Ya que estos instrumentos deberian llegar 
a medir hasta el cut-off de altas energias, la medicion de este valor daria una valiosa 
information acerca de la eflciencia de aceleracion de parriculas en los choques terminales 
de los jets de YSOs de gran masa. 

En rayos X, la fuente deberia ser detectable por Chandra y XMM-Newton, probable- 
mente como una fuente puntual a traves de una observation profunda en la region. En 
los datos de archivo de XMM-Newton, no se ve una emision significativa por encima de la 
radiation de fondo en la region donde la fuente IRAS esta localizada. Esto nos proporciona 
una cota maxima de ~ 10 31 erg s _1 (a la distancia ~ 2.9 kpc) en el rango de 1 a 10 keV 
a las luminosidades que obtenemos con nuestros modelos teoricos, como se muestra en las 
Figuras 14.81 y 14.91 Sin embargo, una observation mas profunda (mas de 30 ks) deberia 
detectar radiation por encima de la emision de fondo de esta region. Si las temperaturas 
del material chocado son suficientemente altas, la componente termica deberia ser tambien 
detectable en rayos X. En la banda ultravioleta (UV), esta emision termica podria ser aun 



3 Los instrumentos HESS II y MAGIC II serian los que reemplazarfan a HESS y MAGIC, respectiva- 
mente. Por otro lado, los instrumentos que formaran parte del arreglo de telescopios Cherenkov CTA (por 
Cherenkov Telescope Array) ya comenzaron a contruirse, aunque el sitio donde estaran montados aiin no 
esta dcfmido. 
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mayor que la producida por las partfculas no termicas, pero sen'a dificil de detectaiQ 

Un hecho relevante respecto de la factibilidad del escenario discutido en este capftulo 
es que las luminosidades no termicas son bastante mas bajas que la luminosidad cinetica 
del jet de la fuente IRAS 16547-4247. La perdida de masa del jet ha sido estimada en 
~ 1O _5 M yr _1 y con una velocidad t>j ~ 1000 km s _1 con lo que la luminosidad cinetica 
del jet resulta ~ 10 36 erg s _1 , mas de tres ordenes de magnitud mayor que la luminosidad no 
termica predicha por nuestro modelo. Esto es, nuestro modelo predice niveles de emision 
que respetan los lfmites energeticos de la fuente. Por lo tanto, todavfa serfa posible un 
incremento de la eficiencia radiativa, incrementando asf las posibilidades de detecteccion 
de este tipo de objetos. Ademas de la fuente IRAS 16547-4247, otras protoestrellas de gran 
masa son tambien potenciales candidatos a fuentes de rayos gamma. Entre ellas podemos 
mencionar HH 80-81 (Bosch-Ramon et al. 2010) y W3(OH) (Araudo et al. 2008a). 



4 Para una discusion mas detallada sobre la emision termica, como asf tambien un analisis mas preciso 
del impacto del medio en la evolution de las partfculas y en la radiation producida, ver Bosch-Ramon et 
al. (2010). 
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Capitulo 5 
Microcuasares 



5.1. Introduccion 

Las estrellas en general no se encuentran aisladas, sino que la mayoria de ellas forman 
sistemas binarios o cumulos. En los primeros, si ambas componentes tienen masas diferen- 
tes, la estrella mas temprana evolucionara mas rapido que la mas tardfa y asi las estrellas 
estaran en diferentes estados evolutivos. Si la masa de la estrella mas evolucionada que 
queda luego de la explosion de supernova es 3M > M* > 1.4M , entonces el objeto re- 
manente al final de la evolucion sera una estrella de neutrones, mientras que si M+ > 3M Q 
sera un agujero negro. 

Consideremos un sistema binario formado por una estrella no degenerada, la cual puede 
estar en diferentes estados de su evolucion, y el objeto colapsado, el cual puede ser un 
agujero negro o una estrella de neutrones. Cuando la estrella companera se convierte en 
gigante roja, las capas externas de la atmosfera de la misma llenan de material el lobulo de 
Roche. La materia aqui dentro se escapa a traves del punto lagrangiano LI siendo acretada 
por el objeto compacto. La materia cae a la superficie del objeto a traves de un disco de 
acrecion debido a la conservation del momento angular. La materia en el disco rota a una 
velocidad que disminuye con la distancia r& al objeto compacto. Anillos a diferentes r& rotan 
a velocidades distintas lo que produce disipacion de energia y calentamiento por roce de 
anillos contiguos. La materia del disco se puede calentar hasta temperaturas suficientemente 
altas como para radiar termicamente en rayos X. 

Los sistemas binarios que emiten fuertemente en rayos X debido a la acrecion, son 
llamados binarias de rayos X (XRB, por X-Ray Binaries). Si la estrella primaria es de 
gran masa (tipo espectral O 6 B, M+ > 8M ) entonces el sistema se dice de gran masa 
(HMXB, por High Mass X-ray Binary), mientras que si es vieja el sistema binario es una 
binaria de rayos X de baja masa (LMXB, por Low Mass X-ray Binary). Una subclase de 
las XRB son los MQs, los cuales presentan jets extendidos en la banda de radio (Mirabel & 
Rodriguez 1999). Estos jets son flujos colimados de materia que se mueve a una velocidad 
de co nj unto relativista. La materia que forma los jets puede ser materia arrancada del 
disco de acrecion por las lfneas de campo magnetico (Blandford & Payne 1981) o bien 
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Figura 5.1: Esquema de un MQ. Se muestra una estrella evolucionada transfiriendo masa al 
objeto compacto a traves del disco de acrecion. Jets bipolares emergen del objeto compacto. 
(Credito: ESA, NASA y Felix Mirabel.) 



pares generados en la ergosfera del agujero negro (Blandford & Znajek 1977) y lanzados 
en forma de chorros por procesos MHDs. El campo magnetico juega un rol importante en 
la formation y colimacion de los jets. 

Se ha detectado emision no termica proveniente de los jets de los MQs. Esta emision 
abarca desde las frecuencias radio (Ribo 2005) hasta los rayos X (Corbel et al. 2002), 
aunque emision de mas alta energfa (rayos gamma) tambien puede producirse (Kaufmann 
Bernado et al. 2002, Bosch-Ramon et al. 2006). Los telescopios MAGIC y AGILE han 
detectado emision transitoria (flares) en HE y VHE asociada a los HMMQs Cygnus X-3 
(Tavani et al. 2009b; Abdo et al. 2009a) y Cygnus X-l (Albert et al. 2007; Sabatini et al. 
2010), respectivamente. Eventos de emision transitoria en rayos gamma tambien han sido 
observados en las HMXBs LS 5039 y LS 1+61 303 por los telescopios HESS (Aharonian 
et al. 2005) y MAGIC (Albert et al. 2006), respectivamente, aunque la naturaleza de las 
fuentes no ha sido aun confirmada en estos casos. (Una discusion sobre el tema puede 
encontrarse en Romero et al. 2007.) Por otro lado, el satelite Fermi ha detectado emision 
esporadica proveniente de fuentes no identificadas del piano galactico (ATeQ 1394, Abdo 
et al. 2009). Esta emision altamente variable detectada en rayos gamma y generada en 



1 La sigla ATel (por Astronomical Telegrams) se utiliza para designar los reportes sobre des- 
cubrimientos observacionales. La pagina de internet en la cual pueden hallarse estos reportes es 
http: / /www. astronomerstelegram.org/. 
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HMMQs podria tener un origen similar. Por ejemplo, la interaction entre el jet y el viento 
de la estrella companera podria producir flares en el rango de las HE y VHE. 

La perdida de masa de las estrellas de gran masa forma un viento estelar que se propaga 
a velocidades supersonicas. Observaciones en rayos X de las lmeas de emision producidas 
en el viento indican que este no tiene una estructura uniforme sino porosa, con clumps (o 
grumos de materia) que se forman a una distancia del orden de un radio estelar de 
la superficie de la estrella (e.g. Puis et al. 2006, Owocki & Cohen 2006). Sin embargo, las 
caracteristicas de estos clumps no se conocen ya que la resolution espacial de los telescopios 
aun no es suficiente. Por esto las propiedades de los clumps, como el tama no, la densidad 
y el factor de llenado, no son bien conocidas y se estudian a traves de metodos indirectos 
como el analisis de lmeas espectroscopicas (Moffat 2008). 

En este capftulo proponemos un modelo para explicar los flares en rayos gamma de- 
tectados en algunas binarias de gran masa, basado en la interaction de los jets del HMM Q 
con las inhomogeneidades del viento de la estrella companera. Los clumps pueden eventual- 
mente penetrar en el jet y producir emision no termica transitoria, generada al convertir 
parte de la energia cinetica del jet en energia de particulas relativistas, las cuales emitiran 
radiation sincrotron, IC y sufriran colisiones pp. 



5.2. Escenario 

Para estudiar la interaction entre un clump del viento estelar con uno de los jets de 
un HMMQ adoptamos un escenario con caracteristicas similares a las del sistema binario 
Cygnus X-l. Fijamos la separation a mq entre el objeto compacto y la estrella masiva en 
3 x 10 12 cm (~ 0.2 UA). Para la luminosidad y la temperatura de la estrella adoptamos 
los siguientes valores: = 10 39 erg s _1 y = 3 x 10 4 K, respectivamente. Por otro 
lado suponemos que la perdida de masa de la estrella es M* = 3 x 1O~ 6 M yr 1 con una 
velocidad terminal del viento v v ~ 2.5 x 10 8 cm s _1 . Un esquema grafico de este escenario 
se muestra en la Figura I5T2"1 

Modelo para el clump 

Debido a las incertezas en la determination de las caracteristicas de los clumps, su- 
ponemos que estos son esfericos y consideramos dos valores para el radio: R c = 10 10 
y 10 11 cm (~ 3 x 10~ 3 — 3 x 10~ 2 a mq ). A su vez adoptamos una distribution unifor- 
me de densidad en el clump, n c = 10 12 cm -3 , que corresponde a un factor de llenado 
f c = M*/47r m p v c n c ~ 0.005, donde v c es la velocidad de los clumps y consideramos 
que estos se mueven a la velocidad del viento, esto es, v c = v v (Owocki et al. 2009 y refe- 
rencias en ese articulo). La temperatura de los clumps, T c , es fijada en 10 4 K (Krticka & 
Kubat 2001), siendo moderadamente menor que la temperatura superficial de la estrella 
primaria. 
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Figura 5.2: Esquema de un MQ. El viento de la estrella companera presenta clumps, algunos 
de los cuales pueden llegar a penetrar en el jet (Romero et al. 2008). 
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Tabla 5.1: Parametros adoptados y estimados para el escenario en el cual desarrollamos el 

estudio presentado en este capitulo. 

Parametro [unidades] Clump Jet 

Radio [cm] R c = 10 10 - 10 11 R i = 1.5 x 10 11 

Velocidad [cm s" 1 ] v c = 2.5 x 10 8 v ] = 10 10 

Densidad [cm -3 ] n c = 10 12 n- s = 4.7 x 10 7 

Sistema binario 



Tamano del sistema [cm] 




= 3 x 10 12 


Luminosidad de la estrella [erg s _1 ] 


L* = 


10 39 


Temperatura de la estrella [K] 


T it = 


3 x 10 4 


Tasa de perdida de masa [M yr -1 ] 


M* = 


= 3 x icr 6 


Velocidad del viento [cm s _1 ] 


v v = 


2.5 x 10 8 



Modelo para el jet 

Para los jets, adoptamos un modelo hidrodinamico, es decir, consideramos un jet do- 
minado dinamicamente por protones frios con una velocidad de conjunto moderadamente 
relativista: Vj = 0.3 c, que corresponde a un factor de Lorentz Tj = 1.05. Observaciones 
en radio muestran que los jets de los MQs son extremadamente colimados (Miller- Jones 
et al. 2006). En esta tesis hemos supuesto que la relation entre el radio y la altura es 
-Rj(zj) = 0.1 Zj, lo que corresponde a un angulo de apertura (f> = 6°. La velocidad de ex- 
pansion del jet resulta v exp = 0.1 Vj, que para un jet hidrodinamico en expansion libre 
implica un numero de Mach en la base de ~ Vj/v exp ~ 10, con lo cual el jet es supersonico. 
La luminosidad cinetica del jet es fijada en Lj = 3 x 10 36 erg s -1 , similar a la estimada 
para Cygnus X-l (e.g. Gallo et al. 2005, Russell et al. 2007). Usando la ecuacion ( I2.43[) 
podemos estimar la densidad de particulas nj del jet en el sistema de referenda del labora- 
torio. A la altura de la interaccion jet-clump, que hemos fijado en z int = a mq /2, obtenemos 
rij = 4.7 x 10 7 cm -3 . De esta manera, el cociente entre la densidad del clump y la del jet 
resulta Xmq = 2.1 x 10 4 . Finalmente mencionamos que despreciamos la curvatura del jet 
producida por la interaccion con el viento estelar. En los jets de HMMQs con luminosidades 
cineticas > 10 36 erg s _1 , la geometria de estos no deberia ser modificada considerablemen- 
te por el viento (Perucho & Bosch- Ramon 2008), aunque este efecto puede ser importante 
en sistemas con jets tipo HH (protoestelares) interactuando con el viento de una estrella 
(Raga et al. 2009). 

Los valores de los parametros del jet y del clump que suponemos o estimamos en este 
capitulo estan listados en la Tabla I5TT1 

5.2.1. Interaccion jet-clump 

Los clumps se forman en la region de aceleracion del viento de la estrella, aproxima- 
damente a una distancia ~ R± de la superficie de la misma (Puis et al. 2006). Algunos 
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de estos clumps pueden llegar hasta el jet y penetrar en el, debido al gran contraste de 
densidades Xmq (= n c /n^). 

Habiendo especificado las caracteristicas que adoptamos para los clumps y los jets de 
nuestro modelo, nos concentramos ahora en la interaction de uno de estos clumps con uno 
de los jets del HMMQ. De acuerdo a los parametros considerados, las escalas de tiempo 
de los procesos dinamicos descriptos en la Section 12.5.11 toman los valores que exponemos 
a continuation y que resumimos en la Tabla I5T21 

El tiempo de penetration del clump en el jet es: 

R c ' '" 



mientras que el tiempo de cruce a una altura z [nt = a mq /2 = 1.5 x 10 12 cm es 

fc ~ 1.2 x 10 3 ( ^- ^ s. (5.2) 

J V1.5 x 10 12 cm/ v ; 

El bow shock en el jet alcanza el estado estacionario en un tiempo muy corto. Para los 
parametros de este sistema y considerando que fji ~ Uj/4 y que Z = 0.2 R c (van Dicke & 
Gordon, 1959) tenemos que 

t h8 ~ 0.8 [ , ^ c 1 (—7^ A 1 s. (5.3) 

Vl0 10 cm J VlO^cms" 1 / v ; 

El tiempo caracateristico de este tipo de interacciones esta dado, como mencionamos 
en el Capitulo [21 por aquel en el cual el choque que se propaga en el clump lo recorre 
completamente. Esta escala temporal esta caracterizada de la siguiente manera en un 
HMMQ: 



t rr ~ 5 x 10 3 



R c \ ( n c \V2 / z int 



10 10 cm/ V10 12 cm- 3 / \1.5xl0 12 cm 



/rj-iy 1/2 / Lj Y 1/2 f v i y 1/2 

X V 0.06 J \3 x lO^ergs^ 1 ) V10 10 cms- 1 / ^ ' ' 

La materia del jet acelera al clump aplicandole una fuerza ~ M c g, donde M c es la masa 
del clump y la aceleracion g ~ 3 x 10 4 y 3 x 10 3 cm s -2 , para R c = 10 10 y 10 11 cm, 
respectivamente. Esta fuerza puede acelerar al clump hasta la velocidad del jet en un 
tiempo 



1.3 x 10 b 



Rc \ ( \ / ^int 



10 10 cm/ V10 12 cm- 3 / V 1 - 5 x 10 12 cm / 
0.06 J V3 x lO^ergs" 1 ) S ' 
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Tabla 5.2: Valores obtenidos para las escalas de tiempo descriptas en la seccion 15.2. H 
usando los valores de los parametros listados en la Tabla 15.11 



Escala de tiempo [s] 


R c = 10 10 cm 


R c = 10 11 cm 


tc 


80 


800 




1.2 x 10 3 


1.2 x 10 3 


tcc 


5 x 10 3 


5 x 10 4 


tbs 


0.8 


8 




7 x 10 10 


7 x 10 11 



Sin embargo, antes de que el jet acelere al clump y este comience a ser arrastrado por el 
flujo, el clump puede ser destruido por las inestabilidades de Rayleigh- Taylor y Kelvin- 
Helmholtz, las cuales crecen en un tiempo ^rt/kh ~ t cc <C t g hasta longitudes de escala 
~ R c (sin considerar el efecto estabilizador que puede llegar a tener el campo magnetico). 

A modo de resumen y de acuerdo a las escalas de tiempo estimadas en los parrafos 
previos, el clump puede penetrar completamente en el jet si t c < t cc . En nuestros calculos 
radiativos no consideraremos los detalles de la penetration del clump en el jet, sino que 
para un tiempo > t c el primero se encuentra completamente dentro del ultimo (esto es, 
consideramos un sistema con simetria cilmdrica como se muestra en la Figura 15. 3p . El bow 
shock se forma rapidamente en un tiempo tbs mucho menor que t c y que t cc . Por ultimo, 
notamos que el clump podria no escapar del jet ya que para los parametros considerados 
y el Zi n t que hemos fijado, tj > £rt/kh- Sin embargo, solo podemos cuantificar el tiempo 
de las inestabilidades de una manera muy somera y no podemos afirmar que el clump se 
destruira necesariamente antes de poder escapar del jet. Simulaciones numericas muestran 
que las escalas de tiempo de las inestabilidades pueden ser varias veces el tiempo de cruce 
del choque en el clump, es decir, ^rt/kh > t cc (Klein et al. 1994). 

Respecto de las propiedades de los choques y dadas las caracteristicas espeefficas de 
este escenario (HMMQ), el choque en el clump es fuerte, radiativo y lento, mientras que 
el bow shock es tambien fuerte pero adiabatico y rapido. Por estas razones, el material 
chocado y calentado del clump radia una fraction significativa de la energfa que el choque 
le ha transferido (pero debido a que Xmq 1, esta energia es baja comparada con la que 
transporta el jet). 



5.2.2. Emision termica del clump 

Para estimar la densidad n c \(x) y la temperatura T cl (x) del material del clump chocado 
a una distancia x del choque usamos las relaciones (I2.33P y (12.321) considerando que A(T) = 
7 x i o~ 19 T~ ' 6 erg cm -3 s _1 y que la temperatura y la densidad en la zona adiabatica de la 
region chocada del clump toman los valores 8.5 x 10 6 K y 4x 10 13 cm -3 , respectivamente. En 
la Figura l5~4l se muestran los graficos de n c i(x) y T c \(x). Luego, el tiempo de enfriamiento 
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Figura 5.3: Esquema de la interaction jet-clump (Araudo et al. 2009). 



por radiation termica (de continuo y de lmeas) resulta 



^ter = 3xl0 2 cU / s~10s. (5.6) 



n cl (x) 



Siendo t tei < t cc podemos decir que el choque en el clump es radiative Este tiempo de 
enfriamiento corresponde a una distancia x ia< ^ ~ t^ cv v cc / 4 ~ 2 x 10 cm, que es menor que 
2R C . La velocidad v cc de propagation del choque en el clump es calculada a traves de la 
ecuacion ( I2.47P y resulta v cc ~ 7 x 10 7 cm s _1 . A una distancia > x ra d la temperatura del 
material chocado del clump es muy baja y la densidad crece hasta valores ~ 10 14 cm -3 . 

Aunque en este capitulo estamos interesados en la emision de rayos gamma produci- 
da por la interaction jet-clump, hemos estimado por completitud la radiation libre-libre 
generada por el material chocado y calentado del clump. Considerando n cl (x) y T cl (x) 
estimamos la luminosidad por emision libre-libre integrando la emisividad a lo largo del 
clump chocado (Lang 1999). Las luminosidades bolometricas obtenidas son L ter ~ 5 x 10 30 
y 5 x 10 32 erg s _1 para R c = 10 10 y 10 11 cm, respectivamente, con un maximo alrededor 
de los rayos X blandos (~ 1 keV). La luminosidad especi'fica se muestra juntamente con la 
emision no termica en la Figura 15.81 

Contrariamente al choque en el clump, el bow shock es adiabatico y rapido. Por esta 
razon es un lugar propicio para la aceleracion de particulas hasta energfas relativistas, 
como veremos a continuation. 
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Figura 5.4: Variation de la temperatura (izquierda) y la densidad (derecha) en la region 
chocada del clump en funcion de la distancia x al choque. 



5.3. Poblacion de particulas relativistas 

En lo que sigue consideramos que las particulas relativistas son aceleradas en el bow 
shock unicamente. Los electrones y protones acelerados en este choque son inyectados 
siguiendo una distribution 

Q ei , P = K ei , p E-* p exp ■ (5.7) 

Determinamos la constante de normalization K ei p suponiendo que el 25% de la lumino- 
sidad del jet inyectada en el bow shock, L hs ~ (a c /ay) Lj, donde a c = nR^ es la section 
efectivaH del clump, se convierte en potencia de inyeccion de las particulas relativistas. Asi, 
fijando que L ep = 0.25 L hs hallamos 

K ei _ = 0.25 ( —) — ^ — . (5.8) 



5.3.1. Aceleracion de particulas y perdidas radiativas 

Para calcular el tiempo de aceleracion t ac , necesitamos conocer el valor del campo 
magnetico en la region del bow shock. Consideramos dos valores para esta magnitud. En 
primer lugar, estimamos el B^ s resultante de imponer que la densidad de energfa magnetica 
ub = B^ s /(8n) es el 10 % de la densidad de energfa del material chocado del jet, Uji, cuya 
expresion es la siguiente: 

«ji = 2 Pcin = g n i m v v h ( 5 - 9 ) 



2 Despreciamos la region donde el bow shock se hace muy oblicuo y nos focalizamos en aquella en la cual 
el choque es mas fuerte, es decir, en el frente del clump. 
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donde P cin = pv 2 es la presion cinetica de un medio con densidad p que se mueve con 
velocidad v. Fijando 

|^ = 0.1 Mjl (5.10) 

obtenemos _£>bs ~ 150 G, lo que nos da un tiempo de aceleracion para choques perpendi- 
culares tf c ~ 10~ 2 E eitP s (ver la ecuacion f )2.60p ). Por otro lado, hemos adoptado un valor 
mucho mas bajo, B^ s ~ 1 G, para chequear el impacto de considerar un campo magneti- 
co mucho menos intenso que 150 G. Con B\, s ~ 1 G, el tiempo de aceleracion resulta 
ti ~ 0.2£ eiiP s. 

Ademas de acelerarse, las partfculas pueden escapar de la region chocada del jet ya sea 
por perdidas difusivas o convectivas. El tiempo de convection por los costados del clump 
puede estimarse de la siguiente manera: 



R c . I R, 



v hs /A \W 10 cm 



s. (5.11) 



Por otro lado, el tiempo de difusion desde el bow shock hasta el clump considerando regimen 
de Bohm resulta 

Z 2 ( R c V fB hs \ 1 , . 

Notamos entonces que las partfculas mas energeticas pueden difundir hasta el clump antes 
de ser arrastradas por el material chocado del jet. Tanto en el clump como en la region del 
bow shock, los electrones y protones relativistas pierden energfa a traves de los diferentes 
procesos radiativos que hemos descripto en el Capitulo [3] y cuyas formulas aplicaremos a 
continuation. 



Perdidas leptonicas 

En el caso con B^ s = 150 G, la radiation sincrotron es el mecanismo mas eficiente de 
enfriamiento de los electrones en la region chocada del jet, con una escala de tiempo 

1.8 x 10~ 2 

£ sin ~ s, (5.13) 

mientras que si B\, s = 1 G, t S i n ~4x 10 2 /E ei s. Para este ultimo valor de Bf, s el proceso 
radiativo dominante resulta ser la dispersion IC. 

A la altura del jet a la cual estamos considerando que ocurre la interaction jet-clump, 
z [nt = 1.5 x 10 12 cm, la densidad de energia de los fotones provenientes de la estrella 
companera del HMMQ es u pn * ~ 2.4 x 10 -2 erg cm -3 y la energia promedio de estos fotones 
es £ ph „ ~ 3K B T* ~ 10 eV. Para y = E ph *E e /(m e c 2 ) 2 > 1, esto es, E e > 2.5 x 10 10 eV, la 
interaction IC ocurre en el regimen de KN. El tiempo de enfriamiento tanto en el regimen 
de Th como en el de KN esta dado por la expresion (ver (13.231) ): 

, n , (l + 8.3y) (l + l-3y 2 ) 

tlC ~ °' 4 In(l + 0.2j,) (l + 0.5y + 1.3y=») 8 " (5 ' 14) 
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Figura 5.5: Tiempos de aceleracion y de perdidas radiativas (sincrotron e IC) para elec- 
trones en la region del bow-shock para los casos B hs = 1 G (izquierda) y B hs = 150 G 
(derecha). Los tiempos de convection y difusion se muestran para R c = 10 10 (linea fina) y 
10 11 cm (linea gruesa). 



Por otro lado, las perdidas por Bremsstrahlung relativista no son relevantes en ninguno 
de los dos casos (ni con B^ s = 150 G ni con £?b s = 1 G), ya que la densidad del jet, rij, es 
muy baja a la altura z int , con lo cual la densidad en la region del bow shock tambien resulta 
pequena (n^ s = 4n- 3 ~ 2 x 10 8 cm -3 ). Con este valor de nt, s el tiempo de enfriamiento por 
Bremsstrahlung relativista resulta muy alto (ver (I3.32p ) 

7 x 10 7 . 

t B rem = 7 n S. (5.15) 

In -^2 + 0.36 

Teniendo en cuenta la ganancia por aceleracion, las perdidas radiativas y los tiempos de 
escape calculamos la energia maxima que pueden alcanzar los electrones acelerados en el 
bow shock. En el caso con B^ s = 1 G, la energia maxima esta determinada por las perdidas 
difusivas, mientras que en el caso con B^ s = 150 G es la radiation sincrotron el mecanismo 
de enfriamiento dominante a altas energfas, como se muestra en la Figura [5751 Como puede 
observarse en estos graficos, la energia de quiebre E q del espectro se obtiene igualando 
^conv = tdif , en el caso con B\, s = 1 G mientras que en el caso B^ s = 150 G el quiebre ocurre 
cuando t conv = t sin . Los valores de las energfas maximas y de quiebre obtenidos en cada 
caso se listan en la Tabla [5731 

Teniendo en cuenta las perdidas radiativas que sufren los electrones relativistas en la 
region del bow shock, calculamos la distribution espectral de energia de estos, N ei (E ei ), 
a traves de la ecuacion (12.621) y teniendo en cuenta la constante de normalization del 
espectro de inyeccion dada por la ecuacion (15. 8p . Consideramos un tiempo de duration de 
la interaction r vi d a > t cc y obtenemos un espectro en estado estacionario ya que los tiempos 
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Tabla 5.3: Energfas maximas alcanzadas por las partfculas aceleradas en el bow shock. Los 
valores listados han sido calculados teniendo en cuenta los diferentes valores de R c y B^ s 
considerados en este capftulo. 



R c 


cm] 


10 10 


10 10 


10 11 


10 11 


Bbs 


[G] 


1 


150 


1 


150 


r^max 
h *l 


[eV] 


1.5 x 10 11 


8 x 10 11 


1.5 x 10 12 


8 x 10 11 


rp max 


[eV] 


6 x 10 11 


9 x 10 13 


6 x 10 12 


9 x 10 14 




eV] 


3.1 x 10 10 


2.8 x 10 9 


3.1 x 10 11 


2.8 x 10 8 




Figura 5.6: Distribution de energfa de los electrones acelerados en el bow shock para los 
casos .Bbs = 1 (izquierda) y 150 G (derecha). 



de perdidas radiativas (sincrotron e IC) y de escape son menores que t cc . Los resultados 
obtenidos se muestran en la Figura l5~6l 

Finalmente, los electrones mas energeticos (E ei > 0.3-E™ ax ) se podrian difundir hasta 
el clump. Esto ocurre si tdn < t conv y tdif < £rad> donde la distancia que deben difundirse 
es el grosor del bow shock. Si los electrones llegan al clump, probablemente radien alii toda 
su energfa. 



Perdidas hadronicas 

De la misma manera que el Bremsstrahlung relativista es un canal de enfriamiento muy 
lento para los electrones en la region chocada del jet, el pp, con un ritmo de enfriamiento 
similar al Bremsstrahlung relativista, es tambien un proceso poco eficiente para los proto- 
nes. En la region del bow shock t pp resulta ~ 10 7 s, ya que la densidad del jet chocado es 
nbs ~ 2 x 10 8 cm -3 . La energfa maxima de los protones relativistas acelerados en el bow 
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shock queda determinada por el tiempo de difusion, dando 

E r =» >< "> u (lofe) (£) eV - < 5 - 16 > 

Los protones con energia -E p > 0.025 -E™ ax pueden difundir hasta el clump, ya que 
para estas energfas tdif < ^conv en la direction perpendicular al choque. Una vez allf, estos 
protones pueden interactuar con el material chocado del clump y emitir rayos gamma en 
un tiempo t pp ~ 500 s. Para que los protones puedan estar confinados en el clump, el campo 
magnetico de este debe ser B c > 7.5 x 10 _3 -£>b s , es decir, B c > 7.5 x 10~ 3 6 1.125 G, para 
el caso con B\> s = 1 y 150 G, respectivamente. Luego, los protones menos energeticos de 
aquellos que lleguen al clump podran ser confinados para valores de B c razonables. Pero 
este lfmite de confinamiento es muy laxo, y en realidad el campo B c necesario para confinar 
a los protones podrfa ser mucho mayor. 

Aquf consideraremos el caso mas conservative para calcular la distribucion de energia de 
los protones en el clump. Es decir, teniendo en cuenta que estos permanecen en el clump solo 
el tiempo que tardan en cruzarlo a una velocidad cercana a c, esto es, t C ruce ~ R c / C ~ 0.3 
6 3 s si R c = 10 10 6 10 11 cm, respectivamente. Siendo t cruce < t pp , la distribucion de energia 
de los protones en el clump resulta 

N p (E p ) = ^Q p (E p ), (5.17) 

donde 0.025 £™ ax < E p < E™ x . 

5.4. Distribuciones espectrales de energia 

5.4.1. Emision asociada al bow shock 

Como mostramos en la seccion de perdidas leptonicas, los mecanismos mas eficientes 
de perdidas radiativas son el sincrotron y la dispersion IC. Sin embargo, debido a que la 
densidad de energia de la emision sincrotron es menor que la correspondiente al campo 
magnetico (u&) o a campos de radiacion externos, la emision por SSC no sera importante y 
por lo tanto no la tenemos en cuenta en nuestros calculos. Solo consideramos los procesos 
sincrotron e IC externo. 

Calculamos -Eph£ S m(_E p h) y E p hLic(E p h) para los diferentes valores de B\, s y R c consi- 
derados en este capitulo y los resultados se muestran en la Figura 15.71 Como puede verse 
en los graficos, la componente debida a la radiacion sincrotron es mas luminosa que la co- 
rrespondiente al IC en los casos con B^ s = 150 G, alcanzando luminosidades bolometricas 
-^sin ~ 10 33 y 2 x 10 35 erg s _1 para R c = 10 10 y 10 11 cm, respectivamente. Por el contrario, 
para B^ s = 1 G, el proceso radiativo dominante es la interaccion IC, con luminosidades 
bolometricas Lie ~ 2 x 10 32 y 10 35 erg s _1 para R c = 10 10 y 10 11 cm, respectivamente. Las 
energfas maximas alcanzadas por los fotones emitidos son E™ x ~ 1 TeV. Notamos que 
tanto el espectro debido a la radiacion sincrotron como a la dispersion IC en el caso con 
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Figura 5.7: Radiation sincrotron (lmea segmentada azul) y IC (lmea roja) producida en la 
region del bow-shock, para R c = 10 10 (izquierda) y R c = 10 10 cm (derecha), con = 1 
(lmea fina) y 150 G (lmea gruesa). Los efectos de la absorcion gamma-gamma son mostrados 
(lmea gruesa), juntamente con el espectro no absorbido (lmea punteada fina). 



-Bbs = 150 G se quiebra de manera clara por efecto de las perdidas por convection de las 
particulas de la region del bow shock. Ademas, el efecto de las perdidas IC en el regimen 
de KN endurecen el espectro de los electrones en el caso con R c = 10 11 cm y B^s = 1 G, lo 
cual se ve reflejado en el espectro de los fotones producidos por los mecanismos sincrotron 
e IC. Sin embargo, notamos que los electrones relativistas aunque son arrastrados de la 
region del bow shock antes de que emitan significativamente en frecuencias radio, tambien 
pueden emitir en esta banda de energias en otras regiones del je^H- 

En los calculos de las SEDs hemos tenido en cuenta la absorcion (gamma-gamma) por 
creation de pares en el campo de los fotones producidos por la estrella. Debido a que 
no nos focalizamos en la geometria del sistema HMMQ/observador, hemos supuesto que el 
campo de fotones semilla es isotropico (tambien hemos despreciado los efectos angulares en 
el calculo de la interaccion IC). Como puede verse en la Figura 157ft la absorcion gamma- 
gamma reduce los niveles de emision varios ordenes de magnitud a energias de cientos de 
GeV. Solo en algunos casos, con geometrfas especfficas en la interaccion gamma-gamma, 
la atenuacion puede ser despreciable (Khangulyan et al. 2008). 

5.4.2. Emision asociada al clump 

Las particulas mas energeticas aceleradas en el bow shock pueden difundir hasta el 
clump y radiar allf. En el caso de los protones, aquellos con E p > 0.025£'™ ax pueden llegar 

3 E1 calculo de esta emision escapa a los intereses de nuestro estudio. Por esta razon, no hemos tenido 
en cuenta los efectos de la autoabsorcion sincrotron en el espectro en frecuencias radio y no haremos 
predicciones en este rango de energias. 
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al clump. Por otro lado, en el caso B hs = 1 G, los electrones con E ei > 0.3£'™ ax son los que 
pueden difundir una distancia Z y perder toda su energfa en el clump. 

A la emision pp de los protones que difunden hasta el clump la calculamos usan- 
do las ecuaciones (I3.39|) y (13.411) . y considerando que N p (E p ) esta determinada por la 
formula (I5.17p . La emision por pp alcanza una luminosidad ~ 10 32 erg s _1 a una energfa 
E ph ~ 50 GeV, como se muestra en la Figura |5T8| junto con las contribuciones leptonicas 
de la region del bow shock. Ademas de rayos gamma, en las interacciones pp tambien se 
producen pares e ± y neutrinos de muy altas energfas, siendo la luminosidad de estos ulti- 
mos ~ L pp (Aharonian et al. 2006; Reynoso & Romero 2009). Los pares radian casi 
toda su energfa en el clump por los procesos sincrotron, IC y Bremsstrahlung relativista, 
pero la contribution de los electrones primarios tanto en la region del bow shock como en 
el clump es mayor que la de los pares (Bosch- Ramon et al 2005; Orellana et al. 2007). 

Para estimar la radiation producida en el clump por los electrones primarios acelerados 
en el bow shock y que difunden hasta allf, suponemos dos valores para el campo magnetico 
del clump: B c = 1 y 100 G, siendo el ultimo similar al valor de equiparticion entre la 
densidad de energfa magnetica y termica. A las luminosidades especfficas de cada proceso 
(sincrotron, IC y Bremsstrahlung relativista) las calculamos con las formulas dadas en el 
CapftuloEJ Como se muestra en la Figura la componente sincrotron domina a la debida 
a las interacciones IC en el caso con B c = 100 G, alcanzando una luminosidad ~ 10 34 erg s _1 
(R c = 10 11 cm) a Ep^ > 1 MeV. Por otro lado, para B c = 1 G la componente IC absorbida 
(por gamma-gamma) alcanza una luminosidad similar a la sincrotronica y es ~ 10 33 erg s" 1 . 
La emision por Bremsstrahlung relativista es despreciable en ambos casos. Las SEDs para 
el caso R c = 10 10 cm, y para ambos valores de B c , son morfologicamente similares a 
las correspondientes mostradas en la Figura I5.9[ pero las luminosidades alcanzadas son 
aproximadamente dos ordenes de magnitud menores, ya que estas son oc a c . 

Notamos finalmente que la emision no termica del clump es similar, aunque levemente 
menos intensa, a la producida en la region del bow shock, siendo el espectro mas duro. Esto 
es debido a que las partfculas inyectadas en el clump tienen energfas > 0.3i?™ ax . Por clari- 
dad, no graficamos la emision leptonica del clump jutamente con las demas componentes 
(leptonicas de la region del bow shock y hadronica del clump). 



5.5. Interacciones simultaneas 

Hasta aqui solo hemos considerado la interaccion de un clump con el jet a la altura 
Z'mt = ctmq/2. Sin embargo, muchos clumps pueden estar simultaneamente interactuando 
con el jet a diferentes Zj (Owocki et al. 2009). 

La altura minima 2^™°, a la cual los clumps pueden penetrar completamente dentro 
del jet sin ser destruidos en el proceso, es aquella para la cual t c < t cc . Esto determina 
un z ; ™J n ~ 4 x 10 11 cm. De esta manera, el modelo presentado en este capftulo para la 
interaccion jet- clump es valido solo para z- mt > 4 x 10 11 cm (y para R c < Ry). Considerando 
que el jet presenta una geometrfa conica, calculamos el numero de clumps N c que pueden 
estar simultaneamente dentro del jet. Para esto integramos desde z = z^l n hasta z = a mq 



90 



Microcuasares 





Figura 5.8: SEDs de la emision sincrotron, IC, pp y termica para diferentes valores de B^ s 
y R c ; las curvas de la componente absorbida y no absorbida (tineas finas) de la radiacion 
IC y pp son graficadas. 



5.6 Discusion 



91 



35 
33 
31 
29 

" 27 
m 35 
•i 33 
I 31 

29 

27 

25 

-5 -3 -1 1 3 5 7 9 11 13 15 

log (e [eVJ ) 

Figura 5.9: Emision no termica del clump, para el caso con R c = 10 11 cm y para los dos 
valores de B c asumidos. Las curvas correspondientes a la emision por IC y Bremsstrahlung 
relativista se muestran absorbidas (lmeas gruesas) y no absorbidas (lmeas finas). 



y asumimos que el factor de llenado de clumps en el jet es el mismo que en el viento, 
/ = 0.005. Asf hallamos N c ~ 350 y 0.5 para R c = 10 10 y 10 11 cm, respectivamente. Como 
consecuencia de estos resultados, los flares producidos por la interaccion de un clump con 
el jet seran un fenomeno esporadico para iV c bajo (R c = 10 11 cm), o sera una modulation 
de la emision continua de la fuente, para N c alto (R c = 10 10 cm) (ver Owocki et al. 2009). 
En el ultimo caso, la SED resultante sera morfologicamente similar a las mostradas en 
las secciones previas, pero multiplicada por iV c , como se muestra en la Figura 15.101 Sin 
embargo, notamos que por un lado el jet puede verse dinamicamente afectado si muchos 
clumps estan simultaneamente dentro de el. Por el otro, a diferentes alturas de interaccion, 
la luminosidad no termica disponible para radiar es distinta (disminuyendo con z$) y por lo 
tanto la SED resultante no es un simple escaleo con N c como el mostrado en la Figura IBTTOl 

5.6. Discusion 

En este capitulo exploramos los procesos fisicos mas relevantes y la naturaleza de la 
radiacion producida por la interaccion de un clump del viento de la estrella primaria con el 
jet de un HMMQ. La interaccion jet-clump produce dos choques: uno en el jet y otro en el 
clump. El primero alcanza rapidamente el estado estacionario formando un bow shock en 
el jet mientras que el segundo se propaga a traves del clump empujado por la presion del 
medio chocado, que se equilibra con la presion del material del jet chocado. El bow shock 
es adiabatico y rapido, y particulas cargadas pueden acelerarse hasta energfas muy altas 
alii mediante el mecanismo de Fermi. En la region del bow shock los electrones relativistas 
se enfrian eficientemente por radiacion sincrotron e IC. Por otro lado, el choque en el clump 
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Figura 5.10: Distribuciones espectrales de energia de la emision sincrotron, IC, pp y termica 
para el caso con clumps de R c = 10 10 cm y 5b s = 1 (derecha) y 150 G (derecha). Hemos 
tenido en cuenta la contribution de todos los clumps que simultaneamente se encuentran 
dentro del jet. 



es lento y radiativo, no siendo eficiente para acelerar particulas. Sin embargo, la emision 
termica del material chocado del clump podrfa ser significativa, como asi tambien podrian 
serlo las interacciones pp entre los protones relativistas acelerados en el bow shock que 
difunden hast a el clump y el material chocado de este. Si los electrones acelerados en el 
bow shock llegan hasta el clump, estos pueden a su vez radiar en el clump eficientemente a 
traves de los procesos sincrotron e IC. Las SEDs de las componentes radiativas mencionadas 
anteriormente han sido calculadas en el contexto de un HMMQ con parametros similares 
a los del sistema Cygnus X-l y los resultados se muestran en las Figuras I5T81 y 15.91 

En rayos X, la emision es producida por radiation sincrotron (en el bow shock y en el 
clump) y termica (en el clump). En los casos con B^ s ~ 1 G, la emision termica alcanza 
luminosidades L ter ~ 10 32 erg s -1 siendo asi mayor que la sincrotron emitida en la region 
del bow-shock, pero no mayor que la emitida en el clump si R c = 10 11 cm. La emision 
sincrotron de la region del bow-shock es dominante en rayos X para i?bs = 150 G, alcanzando 
L sin ~ 10 35 erg s -1 . En una fuente como Cygnus X-l, estos niveles de emision en rayos X 
serian superados por la radiation del disco de acrecion. Sin embargo, en el caso de fuentes 
poco luminosas en rayos X, como LS 5039 y LS I +61 303 (Bosch-Ramon et al. 2007, Paredes 
et al. 2007), los rayos X producidos via el proceso sincrotron durante la interaction jet- 
clump deberfan ser detectables, y aiin la componente termica debena ser detectable bajo 
determinadas condiciones (clumps grandes con densidades relativamente bajas). 

Las dispersiones IC en la region del bow-shock y en el clump producen rayos gamma 
hasta VHE, dominando la SED en los casos con campos magneticos relativamente bajos 
(5bs = 1 G). En nuestros calculos, la luminosidad mas alta alcanzada es Ljc ~ 10 35 erg s _1 
para R c = 10 11 cm, aunque la absorcion gamma-gamma puede reducir sustancialmente la 
emision por encima de los 100 GeV (Romero et al 2010). Las interacciones pp en el clump 



5.6 Discusion 



93 



pueden tambien producir rayos gamma a energias tan altas como ~ 10 14 eV (B hs = 150 G). 
La luminosidad maxima obtenida por pp es sin embargo modesta, L pp ~ 10 32 erg s _1 para 
R c = 10 11 cm, aunque clumps mas densos y/o mas grandes, y jets mas poderosos produ- 
cirian cantidades detectables de fotones fuera del rango (0.1-10 TeV) donde la absorcion 
gamma-gamma es importante. Recordamos que una geometria especifica del sistema bina- 
rio/observador mas un emisor de altas energias lejos del objeto compacto puede dar una 
atenuacion de los rayos gamma mucho menor (e.g. Khangulyan et al. 2008). 

Como consecuencia de las caracteristicas de la interaccion, la emision esperada es tran- 
sitoria (tipo flare). La duracion de esta emision esta relacionada con la permanencia del 
clump dentro del jet, lo cual depende fuertemente de las inestabilidades de RT y KH, las 
cuales pueden destruir el clump. Debido a que el clump puede ser acelerado dentro del jet, 
el tiempo de vida del primero es de algunas veces el tiempo caracteristico t cc . Si el clump 
no ha sido destruido, este puede eventualmente salir del jet luego de haber sido chocado y 
calentado. Dadas las escalas de tiempo dinamicas de la interaccion, el evento tendra una 
duracion de entre unos pocos minutos y algunas horas. 

Los flares producidos en las interacciones jet-clump pueden tener asociados componen- 
tes espectrales a bajas (radiation sincrotron y emision termica) y altas energias (interaccio- 
nes IC y pp), las cuales no tienen que estar correlacionadas con la actividad de acrecion del 
disco. El nivel total de emision, la importancia relativa de las diferentes componentes de 
la SED y la duracion de los flares pueden proveer de information sobre la potencia del jet, 
como asi tambien del tamano y la densidad de los clumps y el valor del campo magnetico 
en la region donde ocurre la interaccion (Romero et al. 2007). Por lo tanto, ademas de las 
propiedades del jet mismo, las caracteristicas de los clumps pueden ser testeadas a traves 
de observaciones en HE y VHE (y probablemente tambien en frecuencias radio) de los 
flares producidos en HMMQs, abriendo una nueva ventana del espectro electromagnetico 
para estudiar los vientos de las estrellas de gran masa. 

Dependiendo del factor de llenado / del viento (o de la densidad de los clumps) y de 
R c , el numero de clumps -N c - que simultaneamente pueden estar dentro del jet puede ser 
bajo (< 1) o alto, esto es, un clump eventualmente o muchos simultaneamente. Luego, los 
flares producidos por las interacciones de clumps con los jets en HMMQs pueden ser un 
fenomeno esporadico (N c pequeno) o pueden aparecer como una modulation estacionaria 
en el espectro (flickering) (N c alto). Sin embargo, notamos que el jet puede verse sustan- 
cialmente afectado si muchos clumps estan simultaneamente dentro de el. Asumiendo que 
la ruptura del jet tiene lugar para Oj < N c x o c , para los parametros del viento y del jet 
adoptados en este capitulo, el jet podria ser destruido si R c < 10 10 cm con el valor de / 
adoptado. Sin embargo, calculos mas detallados de la dmamica de la interaccion jet-clump 
son requeridos para clarificar este hecho. En el proximo capitulo, haremos un estudio un 
poco mas minucioso de la interaccion simultanea de muchos obtaculos con jets, pero en el 
contexto de los AGNs. 
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Capitulo 6 

Niicleos de galaxias activas 



6.1. Introduccion 

Los niicleos de las galaxias albergan agujeros negros supermasivos (SMBHs, por Super 
Massive Black Holes) los cuales abarcan un amplio rango de masas: 10 6 < M sm ^ < 
1O 1O M . El proceso de formacion de estos SMBHs es a traves de la captura de material 
del medio circundante, ya sean nubes de gas, estrellas o cumulos de estrellas (Rees 1984). 
Dependiendo del momento angular del SMBH y de la materia, la acrecion sera esferica 
o no. Estos SMBHs pueden estar activos o no, dependiendo de la tasa de acrecion de 
materia. Podemos decir que la actividad es una etapa en la vida de las galaxias que depende 
fuertemente de la cantidad de materia que haya en las cercanias del SMBH. Si la acrecion 
es suficiente como para que se forme un disco y consecuentemente los jets, entonces se dice 
que tenemos una galaxia activa. 

Actualmente, el modelo mas aceptado para describir a los AGN es el que se conoce 
como modelo estandar. Este considera que las galaxias activas son un sistema compuesto 
por un SMBH rotante como objeto central, circundado por un disco de acrecion el cual 
a su vez esta rodeado por un toro de gas y polvo. Ademas el sistema se compone de 
dos jets paralelos al eje de rotacion del SMBH y que se propagan en sentidos opuestos. 
En la Figura 16.11 se muestran las principales componentes de un AGN. De acuerdo al 
modelo estandar, la actividad de los AGN es un fenomeno intrinsecamente anisotropico. 
Consecuentemente, la fenomenologfa observada dependera del angulo de inclination entre 
la lmea de la visual y algiin eje de simetria de la fuente, como por ejemplo el eje de rotacion 
del SMBH. Cuando este angulo es cercano a 90° estamos observando o bien una galaxia 
de lmeas delgadas (NLRG, por Narrow Line Region Galaxy) o bien una galaxia Seyfert 
del tipo II; mientras que si dicho angulo es practicamente nulo, tenemos un blazar. En 
cualquier caso intermedio el objeto observado sera un cuasar radio-silencioso, una galaxia 
de lmeas anchas (BLRG, por Broad Line Region Galaxy) o una galaxia Seyfert del tipo I. 
Por otro lado, las radiogalaxias son AGN que presentan una fuerte emision en frecuencias 
radio. Los jets de estas fuentes no estan alineados con la lmea de la visual, es decir, no 
son blazares, y se clasifican de acuerdo a la luminoisidad de los mismos. Las radiogalaxias 
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Figura 6.1: Modelo estandar de AGN. Se muestran las componentes mas importantes de 
acuerdo a este modelo: el SMBH, el disco de acrecion, el toro, los jets y las regiones donde 
se emiten las lmeas anchas y delgadas (Urry & Padovani 1995). 



6.1 Introduction 
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Figura 6.2: Imagen de la radiogalaxia Centaurus A. La imagen muestra esta galaxia activa 
en radio y en el optico. La emision a bajas frecuencias corresponde a los jets y a los lobulos, 
y se piensa que es producida por el proceso sincrotron. La emision optica corresponde a 
la radiacion termica del gas y del polvo que se encuentran tapando la region nuclear de la 
galaxia. (Credito: NASA.) 



Faranoff-Riley I (FR I) son menos luminosas que las FR II (Fanaroff & Riley, 1974) por 
lo cual los jets de las primeras pueden propagarse una distancia mayor que los jets de las 
FR II antes de ser frenados por el medio externo. 

Una de las caracterfsticas principales de los AGN es que la emision continua abarca 
casi todo el espectro electromagnetico, desde frecuencias radio hasta los rayos gamma. 
Esta radiacion proviene basicamente del disco (termica) y de los jets (no termica). En la 
Figura IBT21 se muestra una imagen compuesta de la galaxia activa Centaurus A (Cen A) en 
radio y en optico. 

Ademas del espectro continuo, los AGN emiten tineas en optico y en UV. El mecanismo 
mas aceptado para explicar la emision de estas lineas es que el medio cerca del SMBH no 
es homogeneo sino que tiene una estructura porosa en forma de nubes o estrellas. La region 
donde se forman las lmeas de emision mas anchas (BLR, por Broad Line Region) esta for- 
mada por nubes de material del disco o del medio circundante. Este material esta confmado 
por el medio externo caliente, cuya temperatura es ~ 10 8 K, (Krolik et al. 1981) o por 
campos magneticos (Rees 1987). La materia que forma estas nubes puede ser ionizada 
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por los fotones emitidos en el disco de acrecion produciendo tineas de emision. Luego es- 
tas tineas son ensanchadas debido a que las nubes se mueven en el pozo de potencial del 
SMBH con una velocidad v n > 1000 km s -1 . Otro modelo para la production de las tineas 
anchas detectadas en algunos AGN es aquel en el cual la BLR esta compuesta por estrellas 
evolucionadas (gigantes rojas) cuyas cromosferas son fotoionizadas (Penston 1988). 

La BLR rodea al agujero negro y por lo tanto la interaccion de algunas de las nubes que 
la componen con la parte mas interna de los jets es factible. En este capftulo estudiamos 
la interaccion de nubes de la BLR con la base de los jets, realizando un tratamiento 
similar al desarrollado en el capitulo anterior para los HMMQ. Los choques producidos 
por la penetracion de estas nubes en los jets pueden acelerar particulas hasta velocidades 
relativistas, las cuales luego pueden enfriarse por diferentes procesos no termicos produciedo 
niveles detectables de emision en rayos gamma. La detection de esta radiation nos proveerfa 
information valiosa sobre las condiciones ambientales en las cercamas de la base de los jets 
como asf tambien de las propiedades de la BLR. 

6.2. Escenario 

Bajo ciertas relaciones entre la presion cinetica del jet y la densidad y el tamano de 
las nubes, la penetracion de una nube en el jet es factible. Los detalles del proceso de 
penetracion en si mismo son complejos, y en esta tesis no tratamos esto en detalle aunque 
supondremos que la penetracion ocurre si se satisfacen ciertas condiciones. Un esquema de 
la interaccion se muestra en la Figura 16.31 

A diferencia del tratamiento hecho en el capitulo anterior de la interaccion jet-clump, 
donde consideramos un escenario similar a la fuente Cygnus X-l, en este capitulo no plan- 
teamos un escenario definido, asociado a un AGN especffico, sino que hacemos un tra- 
tamiento mas general de la interaccion. Dejamos fijos a lo largo de todo el capitulo solo 
aquellos parametros estandares que, en principio, no varfan sustancialmente de una galaxia 
a otra. 

Modelo para las nubes 

Supondremos en este estudio nubes esfericas, cuyo radio fijamos en R n = 10 13 cm 
(Risaliti 2009), y con una densidad uniforme tipica n n = 10 10 cm" 3 . La velocidad de las 
nubes esta fijada en v n = 10 9 cm s -1 (Peterson 2006). 

Modelo para los jets 

Suponemos que los jets son hidrodinamicos y relativistas, con un factor de Lorentz 
Tj = 10, lo cual implica una velocidad Vj ~ c. A su vez fijamos el angulo de apertura 
cf) ~ 6°, es decir, la relation entre el radio y la altura del jet resulta i?j = tan(0) Zj ~ 0.1 z- y 
La densidad de los jets a una altura z- } y en el SR del laboratorio puede estimarse a traves 
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Tabla 6.1: Valores asumidos en este trabajo para las nubes de la BLR y para los jets. 



Description Valor 



Tama no de las nubes 


R n 


= 10 13 cm 


Densidad de las nubes 


n n 


= 10 10 cm" 3 


Velocidad de las nubes 


v n -■ 


= 10 9 cm s -1 


Temperatura de las nubes 




= 2 x 10 4 K 


Factor de Lorentz de los jets 




= 10 


Angulo de semi-apertura de los jets 




i 6° 



de la relation (12.431) . con la cual obtenemos 

"< = m-i^ " 8 x 104 (ic4^) pf 1 ) 1 (io^)" 2 cm " 3 - (61) 

donde Lj es la luminosidad cinetica de los jets J a i — irRf- En la Tabla [67T1 se listan los 
valores de los parametros de las nubes y de los jets que quedaran fijos a lo largo de todo 
el capitulo. 



6.2.1. Interaccion jet-nuhe 

El modelo de interaccion de una nube de la BLR con uno de los jets del AGN es similar 
al descripto en el capitulo anterior pero en este escenario el jet es relativista. 

Una de las condiciones que deben satisfacerse para que la nube pueda entrar entera en 
el jet es que la presion cinetica de este no la destruya en el proceso de penetration. Esto 
significa que el tiempo de penetration, 

tn ~ = 2 X 10 - r S, (6.2) 

v n \10 13 cmy V10 9 cms- l ) V ; 

debe ser menor que el tiempo de vida de la nube dentro del jet. Para estimar esta es- 
cala temporal debemos conocer la velocidad del choque en la nube, v cn , que se obtie- 
ne, como vimos en el Capitulo El igualando las presiones cineticas del jet y de la nube: 
(Tj — 1) rij m p c 2 = n n m p v 2 n , valido mientras v cn <C c . Asi obtenemos 



1/2 



cm 



^ cn i/2 " 3X 10 ll0 10 cm- 3 J ll0 16 cnJ V lO^ergs- 1 J s' (6 ' 3) 

donde Xagn = n n /n^(zy). Luego, el tiempo de cruce del choque a traves de toda la nube 
result a 

«- ~ ^ = 7 x 10* C ^ f (^) ( ^A^) ~"\ . (6.4) 



v cn ' V 10l3cm 7 Vl0 10 cm- 3 7 V10 16 cm/ VlO^ergs- 1 



100 
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Figura 6.3: Esquema, no a escala, del modelo unificado de los AGN a las escalas espaciales 
de la BLR. La interaction entre una nube de la BLR y uno de los jets se muestra en la 
parte superior de la figura. El bow shock en el jet y el choque en la nube tambien son 
esquematizados. 
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Para un tiempo tan corto como ~ t cn , la nube se comporta como un obstaculo efectivo 
para el material del jet. Fijando t n ~ t cn nos permite obtener valores mmimos para Xagn y 
2j, por debajo de los cuales la nube no llega a penetrar efectivamente en el jet. 

Debido a la interaccion con el material del jet, inestabilidades hidrodinamicas afectan 
a la nube. El jet ejerce una fuerza en la nube a traves de la superficie de discontinuidad. 
La aceleracion aplicada a la nube puede estimarse a traves de la ecuacion ( I2.49P y en el 
caso que estamos estudiando resulta 

g = \r = \f' (6 - 5) 

Dada la aceleracion g, las inestabilidades de RT se desarrollaran en la nube con una escala 
de tiempo 




for - \ - = \ A T 2 ■ (6-6) 



Para perturbaciones con longitudes de escala I ~ i? n , que son aquellas asociadas a la 
fragmentation significativa de la nube, el tiempo de crecimiento de la inestabilidad resulta 
t RT ~ t CQ . Por otro lado, las inestabilidades de KH tambien pueden crecer suficientemente 
como para destruir la nube. Dada la alta velocidad relativa entre el jet chocado y el material 
de la nube, t> re i ~ Vj, se obtiene 

tKH ~= ^ • (6.7) 

c 

Para I ~ R n , obtenemos nuevamente que £kh ~ ^cn- Notamos que, dada la ecuacion (I6.5p . el 
tiempo necesario para acelerar la nube hasta la velocidad del choque v cn es ~ t cn , mientras 
que para acelerar la nube hasta Uj el tiempo requerido es 3> t cn . Por esto, antes de que la 
nube comience a moverse conjuntamente con el jet, probablemente sera fragmentada. 

Finalmente, hay dos escalas de tiempo adicionales que tambien son relevantes en nuestro 
estudio: el tiempo de formacion del bow shock, tbs, y el tiempo requerido para que la nube 
cruce el jet, t- y Considerando las ecuaciones (I2.45P y f )2.46p . y que en el caso de un bow shock 
relativista este se separa del obstaculo una distancia Z ~ 0.3 R n (obtenido considerando 
que en un plasma relativista las particulas se escapan de la region chocada del jet a la 
velocidad del sonido c/y/3), estas escalas de tiempo resultan 

t ^i = 2xW e (^—) ( /" Y 1 s, (6.8) 
J v n V10 16 cm/ VlOQcms" 1 / V ; 

y 

t hs ~ - = 10 2 (-^-) s. (6.9) 
c \10 id cm/ 

A modo de resumen del estudio anterior sobre la dinamica de la interaccion jet-nube, 
graficamos en la Figura l6\4l t cv (para diferentes valores de Lj), tj, t n y tbs en funcion de z- } . 
Como se muestra en la figura, para algunos valores de z- } y Lj la nube podn'a ser destruida 
por el jet antes de entrar completamente, es decir, t cn < t n . Esto nos provee de una 
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Figura 6.4: Los tiempos de cruce de la nube a traves del jet (lmea punteada azul), de 
penetracion de la nube (lmea punteada-rayada roja), de formacion de bow-shock (lmea 
rayada violeta) y el tiempo en el cual el choque en la nube la recorre entera (lmeas verdes) 
son graficados. Todos ellos han sido calculados usando los valores dados en la Tabla 16. 11 
El tiempo t crL es graficado para Lj = 10 44 , 10 46 y 10 48 erg s _1 (Araudo et al. 2010). 



condicion para determinar la altura del jet a la cual la nube puede penetrar entera dentro 
de el. Notamos tambien que, en general, t hs es mucho mas corto que las demas escalas de 
tiempo. 

Altura de la interaccion 

La nube puede penetrar completamente en el jet si el tiempo de vida de la nube despues 
del impacto con el jet es mas largo que el tiempo de penetracion. Ademas, la condicion de 
que la presion lateral del jet sea < n n m p se satisface automaticamente. Esto determina 
la altura de interaccion minima, z™™, que no permite la destruccion de la nube antes de 
la penetracion completa. Por otro lado, la interaccion no puede ocurrir mas abajo que la 
region de formacion del jet, que se da en z ~ 100 R g ~ 1.5 x 10 15 (M smbh /10 8 M ) cm 
(Junor et al. 1999), donde R g es el radio gravitacional del agujero negro. Finalmente, 
para que la interaccion entre nubes de la BLR y los jets pueda occurrir, a una altura que 
llamaremos z int , el tamano de la BLR debe ser i?bir > Zint > z o- 

El tiempo de vida de la nube depende del tiempo de fragmentacion (dado por el cre- 
cimiento de las inestabilidades), que esta fuertemente relacionado con t cn . El valor de z- m t 
puede luego ser estimado fijando t n < t cn , ya que queremos que la nube entre completa- 
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mente en el jet antes de ser sustancialmente distorcionada por el impacto con este. Una 
vez chocada, la nube puede sufrir expansion lateral y calentamiento por conduction lo cual 
haria que las inestabilidades crezcan mas rapidamente. Sin embargo no hemos tenido en 
cuenta estos dos ultimos efectos en nuestro estudio y determinamos z; nt imponiendo que 
t cn = 2 t n , con lo cual 



W 5 X 10l U 10^1] llb^O ( lO^ergs-0 ^ ((U0) 

Notamos que la potencia disponible en el bow shock es L bs ~ (<7 n /oj)Lj oc 2j~ 2 , donde 
cx n = 7r i? 2 , por lo cual la interaction mas luminosa tendra lugar en Zj ~ z- mt . Ahora debemos 
verificar si la election de z int que hemos hecho cumple con los requisitos geometricos z < 

z int < -Rblr- 

El tamano de la BLR puede estimarse a traves de las relaciones empfricas obtenidas 
para galaxias FR II, las cuales poseen una BLR bien determinada por las observaciones. 
Los ajustes que se obtienen a partir de las observaciones de numerosas galaxias son en 
general del tipo Rur oc L^ r , donde Lbi r es la luminosidad de la BLR y aur ~ 0.5 — 0.7 
(Peterson et al. 2005; Bentz et al. 2006). En este capitulo usamos las siguientes relaciones: 



R b]x ~ 6 x 10 16 ( AA LhlT : ) cm, (6.11) 
1 lO^ergs- 1 ' 



0.7 



0.55 

,16 / ^blr 



y 

R hlT ~ 2.5 x 10 ie ( ^ ) cm, (6.12) 
\10 44 erg s _i y 

obtenidas por Kaspi y colaboradores (2005, 2007). En galaxias FR I, en las cuales la 
detection de la BLR es muy imprecisa o ni siquiera se llega a detectar, las relaciones (16. lip 
y (I6.12p deben tomarse con cautela. 

En la Figura 16.51 mostramos como varian z int y i^bir con Lj , asumiendo que Lbi r es un 
10 % de la luminosidad del disco, L d , y esta ultima es tomada igual a L y Como puede verse 
en la figura, para valores razonables de los parametros, la condition z int < i? blr se satisface 
en un amplio rango de valores de Ly En la misma figura tambien mostramos la relation 
entre z y M sm bh, de la cual se deduce que para M sm bh > 10 9 M el jet podria no estar 
(completamente) formado a las escalas de la BLR para valores bajos de L y 



6.3. Poblacion de particulas relativistas 

En el bow shock y en el choque en la nube se pueden acelerar particulas a traves de 
un mecanismo de aceleracion difusiva como el descripto en el Capitulo 121 (Section 12.61) . Sin 
embargo, debido a que el bow shock es mas fuerte que el choque en la nube (t>b s ^ 
la aceleracion de particulas sera mas eficiente en el primero. Ademas, la luminosidad del 

1 /2 

choque en la nube es menor que la luminosidad del bow shock por un factor ~ l/(2x a gn), con 
lo cual la energfa disponible para entregar a las particulas aceleradas es menor tambien. Por 
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Figura 6.5: La altura de interaction, z- m t (lmea solida roja), y el tamano de la BLR, i?bir 
(lmea rayada verde), se grafican para diferentes valores de L s (eje horizontal inferior). 
Derivamos Ru r (Lj) fijando que Lbir = 0.1 L s y graficamos Rut usando las ecuaciones (16. lip 
y (16.121) . En la misma figura, la altura de la base del jet, z$ (lmea punteada azul), es 
graficada en funcion de M smbh (eje horizontal superior). 
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estas razones, solo consideraremos la aceleracion de electrones y protones en el bow shock 
y en esta seccion describiremos someramente la inyeccion y evolution de estas particulas 
no termicas. 

La luminosidad inyectada por el jet en el bow shock localizado en z int es L^ s = (cr n / cr^)L y 
Una fraction rj nt de esta cantidad sera transferida a las particulas que se aceleran en este 
choque, con lo cual la luminosidad de estas es 



L e , p = Vnt L hs « 4 x 10 39 (H) (, 



Rn \ ( L\ \ erg 
" x 1 ' » — . (6.13) 



10 13 cm/ VlO^ergs- 1 / s 

La distribution de energfa de la poblacion de particulas relativistas inyectadas en la 
region del bow shock es una ley de potencias con un cut-off exponencial: 

Q e , P = K e , p E~™ exp-*-/*- . (6.14) 

A la constante K ep la determinamos a traves de la ecuacion (16.131) y considerando que 

pmax 
/ e,p 

L e , P — / Qe, P E e ^dE ep , (6.15) 

J Emin 

donde fijamos E™™ = 2m 6iP c 2 pero E™^ es desconocida aiin. Sin embargo, igualando las 
ecuaciones (16. 14|) y (16 . 1 5[) obtenemos 

(0.2 \ 0.2 

^ . _ 02 — nn ~ Vnt L hs — — 02, (6.16) 
&e,p ~ ^e,p J &e,p 

S1 ^e,p > ^e,p ■ 

6.3.1. Aceleracion de particulas y perdidas radiativas 

El campo magnetico -Bbs en el acelerador/emisor puede determinarse relacionando las 
densidades de energfa magnetica y no termica de la siguiente manera: ub = )]B«e,p, donde 
la densidad de energfa de las particulas relativistas es u e ^ p = L e>p / (<r n c). De manera tal que 
la radiation IC sea importante en rayos gamma, t]-q <C 1 es requerido. En este contexto, si 
?7b = 0.01, -Bbs puede se parametrizado de la siquiente manera: 

B hs ~ 10 (*L) ( J?B V ( "° \ G . (6.17) 
V0.01/ Violins" 1 / V10 10 cm- 3 / V ; 

Debido a que el bow shock es relativista y el tratamiento de estos choques es complejo, 
adoptamos la prescription E™ ~ 0.1 eB^c erg s _1 (de Jager et al. 1996) para la tasa de 
aceleracion, con la cual el tiempo de aceleracion esta dado por la relation 

-Bbs N 

■-...II / H I 



°- 7E ^\ wk< s - (6 - 18) 
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Las particulas estan afectadas por diferentes tipos de perdidas que compensan las ga- 
nancias por aceleracion. Las perdidas por escape de las particulas contemplan la convection 
del material chocado, 

Wv~^ = 10 3 (-^-) s, (6.19) 
c \10 id cm/ 

y la difusion, la cual para que las particulas puedan difundir una distancia Z = 0.3i? n el 
tiempo requerido es 

5xl ° 6 (lo^) 2 (^)^ s - < 620 > 

Perdidas leptonicas 

Los electrones en la region chocada del jet radian principalmente por los procesos 
Bremsstrahlung relativista, radiacion sincrotron y dispersiones IC. En el ultimo caso los 
fotones semilla pueden ser externos (EC), producidos en el disco o en la BLR, o generados 
en la misma fuente por radiacion sincrotron (SSC). Dada la baja densidad del jet en z int 
(rij ~ 3.2 x 10 5 cm~ 3 ) las perdidas por Bremsstrahlung relativista seran despreciables, con 
una escala de tiempo ts rem ~ 10^ s. Considerando el valor de dado en la ecuacion ( 16.171) . 
las perdidas por radiacion sincrotron resultan: 

* sin ~4.i(^r 2 -L s 

V0.01/ E e 

Para estimar la relevancia de las interacciones IC (EC o SSC), necesitamos conocer la 
densidad n p h de fotones semilla para asi luego calcular £ec y ^ssc a traves de las formulas 
dadas en el Capitulo [31 A continuation describimos como calculamos n p ^ de cada campo 
de fotones semilla considerado en este capitulo. 

■ SSC: Estimando la luminosidad de la radiacion sincrotron, L sin ~ f N e (E e ) P sin (E e , E ph 
donde V(E e , E ph ) esta dada por la formula ( 13.121) y considerando que la distribution 
iV e se encuentra en el estado estacionario, podemos hallar la densidad de energfa de 
los fotones sincrotron: M sin ~ L S i n /(a n c). Como la distribution de estos fotones es 
una ley de potencia, la densidad de fotones es n p ™ = d£7 S m/d-Eo m . Para calcular tssc 
desarrollamos un codigo numerico que describiremos someramente mas abajo. 

■ EC (disco): La densidad de energfa de los fotones emitidos por el disco pue- 
de estimarse de la siguiente manera: u d ~ L d / zf Qt c) . De acuerdo al modelo de 
disco delgado de Shakura & Sunyaev (1973), la temperatura T d del material acre- 
tado decrece con la distancia r d al SMBH de acuerdo a T d (r d ) = T in (r d /r in )~ 3 / 4 , 
donde r in = 3r g = 4.5 x lO 5 (M smbh /M ) cm y T in = T d (r in ). Debido a que la ma- 
yor intensidad de radiacion se produce en las partes mas internas del disco, esto 
es, en r d ~ r in , podemos considerar que T d ~ T in y luego L d ~ 47r cjsb T^rf n , donde 
csb = 5.67 x 10~ 5 gr K~ 4 s -2 . Despejando de aqui la temperatura, la energfa de los fo- 
tones emitidos por el disco resulta E$ ~ 2.4 x 10 2 (L d /10 44 ) 1/4 (M smbh /10 6 M Q )- 1 / 2 eV 



(6.21) 
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y los electrones con E e > (L d /lO 44 )- 1/4 (M smbh /lO 6 M ) 1 / 2 GeV interactuan con los 
fotones del disco en el regimen de KN. 

■ EC (BLR): La densidad de energfa de los fotones emitidos por la BLR es calculada 
mediante uu r ~ -^bir/( 7r -Rbir c ) ~ 0.08(L b i r /10 44 ) -0,4 , donde hemos considerado la 
relation (16. lip . Para nubes con temperaturas T n ~ 2 x 10 4 K, la energfa de los 
fotones resulta Eq 1 * ~ 5.3 eV y la interaction en el regimen de KN se dara para 
dispersiones con electrones cuya energfa sea E e > 48 GeV. 

En la Figura I6U1 (izquierda) se muestran las escalas de tiempo de las perdidas radiativas, 
junto con el tiempo de aceleracion y de convection para un bow shock ubicado en z int . Para 
tjb adoptamos un valor igual a 0.01 y las perdidas por SSC son graficadas para el estado 
estacionario de NE e (E e ). Notamos que tanto B\> s como L e ^ p y z int son constantes para 
diferentes valores de Lj y fijadas las fracciones 7/b y r/ nt . Solo Lbi r y se espera que varien 
con L y Asf, mientras los campos de fotones externos no sean relevantes, la energfa maxima 
que pueden alcanzar los electrones en z int no cambia para diferentes potencias del jet. 

Los procesos radiativos dominantes a las energfas mas altas de la poblacion de electrones 
son el sincrotron, el EC y el SSC. El Bremsstrahlung relativista es despreciable para todas 
las energfas y el escape por convection es relevante en la parte de mas bajas energfas. Esto 
produce un quiebre en el espectro de los electrones, cuando las perdidas por sincrotron, 
EC o SSC se igualan a las de convection. La transition del regimen de Th al de KN se ve 
claramente en las curvas de enfriamiento por EC, pero es menos notorio en el caso del SSC 
debido a la naturaleza no termica de los fotones semilla (sincrotron). La energfa maxima 
de los electrones es de ~ 40 TeV (para r/ B = 0.01) y esta determinada por la igualdad 

^ac fgiil- 

A primer orden, la evolution de la distribution de los electrones relativistas en la re- 
gion del bow shock puede ser calculada asumiendo condiciones homogeneas en el acelera- 
dor/emisor y los mecanismos de aceleracion y escape mencionados anteriormente. Debido 
a que el SSC es un canal de enfriamiento importante, el calculo de N e (E e ,t) debe rea- 
lizarse numericamente. Para esto desarrollamos un codigo en el cual para un tiempo de 
observation t, dividimos el perfodo [0-t] en intervalos At. En cada intervalo, N e se calcula 
teniendo en cuenta las perdidas por SSC con los fotones sincrotron radiados en el intervalo 
anterior. Asf se contimia hasta que se alcanza el estado estacionario. La duration de cada 
intervalo debe ser menor que la duration de los intervalos previos de modo de calcular 
apropiadamente el crecimiento de u S i n en el emisor. En la Figura 16.61 (derecha) se muestra 
el resultado obtenido considerando r vi d a > t cn . El estado estacionario se alcanza cuando 

t tconv 

Perdidas hadronicas 

Debido a la baja densidad del jet en z- mt , los protones no se enfrfan eficientemente en 
la region del bow shock por interacciones pp. Considerando que rij(zi n t) ~ 3.2 x 10 5 cm -3 , 
el tiempo de enfriamiento por pp resulta muy largo, t^l ~ 2 x 10 9 s. La energfa maxima de 
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Figura 6.6: Izquierda: Las ganancias por aceleracion (lmea punteada-rayada anaranjada), 
convection (lmea rayada azul) y los tiempos de enfriamiento radiativos son graficados. 
El tiempo SSC (lmea verde de rayas y puntos) graficado corresponde al calculado con la 
distribution de energfa estacionaria de los electrones. El tiempo de enfriamiento por EC 
con los fotones de la BLR (cuadrados turquesa) y del disco (lmea punteada marron) son 
graficados para fuentes debiles (BLR: 10 44 erg s -1 ; disco: 10 45 erg s _1 ) y brillantes (BLR: 
10 46 erg s _1 ; disco: 10 47 erg s _1 ). Las perdidas por radiation sincrotron (lmea roja) y 
Bremsstrahhmg relativista (lmea violeta) son tambien graficadas. Hemos fijado t/b = 0.01. 
Derecha: Evolution temporal de la distribution de los electrones relativistas en la region del 
bow shock, N e . Las lmeas de rayas indican el espectro a diferentes tiempos ti, mientras que 
la lmea llena indica el espectro una vez que el estado estacionario ha sido alcanzado. Como 
puede observarse, el quiebre en el espectro estacionario se produce en E e ~ 2 x 10 9 eV, 
como se deduce de la figura de la izquierda. 
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los protones queda entonces determinada por el tiempo que tardan estos en difundir una 
distancia Z = 0.3i? n desde el bow shock hasta la nube, lo cual nos da una energia maxima 



Los protones con energfas E p > 0.4_E™ ax , es decir, t conv > t dif , podran difundir hasta la 
nube antes de ser convectados por el material chocado del jet. 

En la nube, el campo magnetico necesario para confinar a estos protones tan energeticos 
es £>„ onf > 0.7(-E>bs/10 G). Si los protones estuviesen confinados en la nube, entonces ra- 
diarfan una fraction significativa de su energia por interacciones pp con el material chocado 
de la nube en un tiempo t pp ~ 2 x 10 15 / (4n n ) ~ 5 x 10 4 s. Sin embargo, como no conocemos 
el valor del campo magnetico de las nubes de la BLR, adoptamos que las particulas solo 
permanecen en la nube el tiempo que tardan en cruzarla a una velocidad ~ c, con lo cual 
el tiempo de cruce es t crU ce ~ Rn/c ~ 700 s < t pp . Luego, la distribution energetica de los 
protones relativistas en la nube sera 



6.4. Distribuciones espectrales de energia 

Una vez que las distribuciones de energia de las particulas relativistas (N e (E e ) y N P (E P )) 
han sido calculadas (en el estado estacionario), estamos en condiciones de calcular las SEDs 
de la radiacion no termica. 

6.4.1. Emision asociada al bow shock 

En el bow shock los electrones pierden energia basicamente por radiacion sincrotron, 
SSC y EC (con los fotones del disco y de la BLR), aunque este ultimo no es el proceso 
dominante si -Bbs ~ 10 G. La autoabsorcion sincrotron es tenida en cuenta, pero solo 
afecta a la parte menos energetica de los fotones sincrotron. En rayos gamma, la absorcion 
(gamma-gamma) por creation de pares e ± producida por la emision del disco de acrecion 
y de la BLR tambien debe ser considerada, pero la absorcion interna debida a la radiacion 
sincrotron es despreciable y no la tenemos en cuenta. Dadas las energfas tipicas de los 
fotones emitidos por el disco y la BLR, ~ 1 keV y ~ 10 eV, respectivamente, los rayos 
gamma con energfas por encima de 1 y 100 GeV pueden ser fuertemente afectados por la 
absorcion gamma- gamma. Por otro lado, en la mayoria de los casos los fotones con energfas 
< 1 GeV van a escapar del denso campo de fotones emitidos por el disco. 

En la Figura 16.71 se muestra la SED de los procesos sinctrotron y SSC en la region 
del bow shock. Debido a que este calculo fue hecho para un caso general, sin especificar 
las luminosidades ni del jet, ni de la BLR, ni del disco, la absorcion gamma-gamma no 




(6.22) 




(6.23) 



donde Q p = K p E p 2,2 exp(— Ep/E™*) y K p esta dada por la ecuacion f)6.16p . 
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Figura 6.7: Distribution espectral de energia calculada en el caso general considerando los 
procesos sincrotron y SSC. Puede observarse el efecto de la autoabsorcion sincrotron, a 
E ph < 10~ 3 eV. 



fue tenida en cuenta, pero si la autoabsorcion sinctrotron. Las luminosidades bolometricas 
alcanzadas por ambos procesos de emision son L sin ~ 4 x 10 38 erg s _1 y Lssc ~ 8 x 
10 38 erg s _1 . 

6.4.2. Emision asociada a la nube 

Aunque los protones no se enfrien eficientemente en la region del bow shock, si pueden 
hacerlo en la nube aquellos que con E p > OAE™^ llegan hasta alll Sin embargo, debido 
a que estos protones son muy energeticos y el campo magnetico necesario para confinarlos 
en la nube es muy alto, estos protones cruzaran la nube a la velocidad de la luz. Con lo 
cual, solo la pequeiia fraction ~ 0.17t cruce //f:pp ~ 1CT 3 de la energia promedio por proton 
sera radiada por interacciones pp en la nube y la luminosidad bolometrica emitida resulta 
L pp ~ 2 x 10 36 erg s -1 . Debido a que no hemos considerado confinamiento, este valor es un 
limite inferior para la emision por pp de la nube. 



6.5. Interacciones multiples 

El tamano de la BLR determina la region en la cual pueden ocurrir interacciones jet- 
nube. Debido a que hay muchas nubes en la BLR, es licit o pensar que algunas de ellas 
pueden interactuar simultaneamente con el jet a diferentes alturas z™[ n < Zj < Rhh J cada 
una de las cuales produce una cierta cantidad de radiation no termica. Esta luminosidad 
total puede ser mucho mayor que la producida por la interaction de una unica nube con 
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Figura 6.8: Esquema de la interaction simultanea de muchas nubes con el jet. 
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el jet, que es ~ L ep . El numero de nubes dentro de ambos jets, N^, puede ser calculado a 
traves del volumen de cada jet (Vj) y de cada nube (In), resultando 

donde el factor 2 es debido a la presencia de dos jets y / ~ 10~ 6 es el factor de llenado 
de nubes en la BLR (Dietrich et al. 1999). Este calculo de Af n es correcto si no tenemos en 
cuenta que las nubes pueden ser destruidas dentro del jet y los fragmentos eventualmente 
dilmdos. Por ejemplo, Klein y colaboradores (1994) estimaron que el tiempo de vida de una 
nube chocada es varias veces t cn , y Shin y colaboradores (2008) encontraron que aun un 
campo magnetico debil en la nube puede incrementar significativamente el tiempo de vida 
de estas. Finalmente, aun en el caso de que las nubes se fragmenten, bow shocks fuertes 
pueden formarse alrededor de cada fragmento antes de que estos sean acelerados hasta 
una velocidad ~ v- y Todas estas consideraciones hacen que el niimero real de nubes que 
simultaneamente estan interactuando con el jet sea dificil de estimar, pero este numero 
debe estar entre (t cn /tj) y N^. La presencia de muchas nubes dentro del jet, no solo en 
z int sino tambien a z- 3 mas altos, implica que la luminosidad no termica total disponible en 
la intersection entre el jet y la BLR es 

donde dA^ es el numero de nubes localizadas en un volumen del jet dVj = n (O.lz) 2 dzy 
En ambas ecuaciones (I6.24p y f)6.25p . L blr ha sido fijada en 0.1 Lj, aproximadamente como 
en las galaxias FR II, y R h \ T ha sido derivada usando la ecuacion (16. lip . 

En la Figura 16.91 mostramos estimaciones de la luminosidad en rayos gamma predicha 
para el caso de la interaction simultanea de muchas nubes con el jet. Para esto hemos 
seguido un procedimiento muy simple, asumiendo que toda la luminosidad no termica es 
radiada en forma de rayos gamma. Este puede ser el caso si los tiempos de escape y de 
enfriamiento por radiation sincrotron son mas largos que el tiempo de enfriamiento por IC 
(EC + SSC) para los electrones mas energeticos. Dado que se conoce muy poco de la BLR 
en las galaxias FR I, no consideramos este tipo de fuentes en la figura. 

6.6. Aplicaciones 

En las proximas dos subsecciones presentamos los resultados obtenidos de aplicar el 
modelo expuesto en las secciones 16.21 y 16.31 a dos fuentes caracteristicas: Cen A (FR I, 
unica interaction) y 3C 273 (FR II, multiples interacciones) . 

6.6.1. Galaxias FR I: Cen A 

Cen A es la galaxia activa mas cercana, ubicada a una distacia d ~ 3.7 Mpc (Israel 
1998). Esta fuente ha sido clasificada como una radio-galaxia FR I y como un objeto 
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Figura 6.9: Lmiites superiores para la luminosidad de rayos gamma producida por 
nubes dentro del jet como una funcion de Lj en fuentes FR II. Dos casos son graficados, 
uno asumiendo que las nubes cruzan el jet sin ser fragmentadas (lmeas solidas verdes), y 
otro en en el cual las nubes son destruidas en un tiempo ~ t cn (lmeas rayadas verdes). Las 
lmeas gruesas (solidas y rayadas) y delgadas (solidas y rayadas) corresponden a i?bir oc L^ T 
(Kaspi et al. 2007) y .Rbir 

oc Zv^j^ 5 (Kaspi et al. 2005), respectivamente. Ademas, los niveles 
de sensibilidad de Fermi en el rango 0.1-1 GeV (lmeas punteadas marrones) son graficadas 
para tres distancias diferentes d = 10, 100 y 1000 Mpc (Araudo et al. 2010). 
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Seyfert 2 en el optico. La masa del agujero negro es M smbh »s 6 x 10 7 M (Marconi et 
al. 2000). El angulo entre los jets y la lmea de la visual es grande, > 50° (Tingay et al. 
1998), de tal manera que la radiacion producida en el jet no esta significativamente corrida 
en frecuencia por efecto Doppler. Los jets de Cen A son frenados a una distancia ~ kpc, 
formando dos radio lobulos gigantes que se extienden ~ 10° en el cielo del hemisferio sur. En 
el optico, la zona del nucleo de Cen A esta oscurecida por una region densa de gas y polvo, 
formada probablemente en una colicion reciente con otra galaxia (Thomson 1992, Mirabel 
et al. 1999). A energias mas altas, los satelites Chandra y XMM-Newton detectaron emision 
continua en rayos X proveniente de la region nuclear, con una luminosidad ~ 5 x 10 41 erg s _1 
entre 2-7 keV (Evans et al. 2004). Estos rayos X pueden haber sido producidos en el disco 
de acrecion y en el jet mas interno, aunque su origen es aiin desconocido. En HE, Cen A 
ha sido detectada por encima de las 200 MeV por Fermi, con una luminosidad bolometrica 
de ~ 4 x 10 40 erg s _1 (Abdo et al. 2009b, y por encima de ~ 200 GeV por HESS, con una 
luminosidad bolometrica de ~ 3 x 10 39 erg s" 1 (Aharonian et al. 2009). En ambos casos, la 
emision en HE es asociada a la region nuclear. Cen A ha sido propuesta como una fuente 
de rayos cosmicos de alta energfa por Romero y colaboradores (1996). 

Si bien no se han detectado aun lfneas de emision que indiquen la presencia de una 
BLR en Cen A (Alexander et al. 1999), esto puede ser una consecuencia de que la region 
nuclear de esta fuente esta oscurecida en el optico por la "estela de polvo". Es por esto 
que puede haber nubes circundando el SMBH de Cen A (Wang et al. 1986, Risaliti et al. 
2002) pero, debido a que el disco de acrecion de este AGN es debil, no se espera un nivel 
alto de fotoionizacion de estas posibles nubes y por lo tanto el proceso sera ineficiente 
para producir lfneas. Bajo esta hipotesis podemos considerar que existe una poblacion de 
nubes oscuras en la region nuclear de Cen A y aplicar el modelo desarrollado a esta fuente, 
considerando interacciones EC solamente con los fotones emitidos por el material acretado. 

Adoptando Lj = 10 44 erg s -1 para la luminosidad del jet de Cen A y aquellos valores que 
han sido listados en la Tabla [6TT1 para el resto de los parametros, z- mt resulta m 5 x 10 15 cm. 
A esta altura del jet, la emision producida por la interaction con una nube es calculada su- 
poniendo que r] B = 0.01. La correspondiente SED se muestra en la Figura I6TTU1 (izquierda) . 
La parte de mas baja energfa del espectro sincrotron es autoabsorbida a energias menores 
que ~ 10~ 4 eV. En rayos gamma, la absorcion gamma-gamma se desprecia debido a que 
la densidad de fotones ambiente es muy baja (e.g. Rieger & Aharonian 2009, Araudo et 
al. 2009, 2010). En HE, la emision por SSC domina la SED, siendo la luminosidad de este 
proceso a energias mayores que 100 MeV ~ 2 x 10 39 erg s -1 , y por encima de 100 GeV 
aproximadamente 10 veces menor. Estos valores son un orden de magnitud menores que 
los correspondientes a los ffujos detectados por los telescopios Fermi y HESS. Notamos 
sin embargo que L 6iP oc R\, y para nubes apenas mas grandes L e , p puede aumentar hasta 
niveles detectables. La penetration de una nube de tamano R n > 10 13 cm en el jet de 
Cen A podrfa producir un flare de aproximadamente un dfa de duracion. 
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Figura 6.10: Distribuciones espectrales de energfa calculadas para las fuentes Cen A (iz- 
quierda) y 3C 273 (derecha). En ambos casos se han graficado las luminosidades especfficas 
producidas por los procesos sincrotron, EC, SSC y pp obtenidos de la interaction de una 
nube con uno de los jets. Tambien mostramos la emision detectada por Fermi y HESS y 
el espectro emitido considerando muchas nubes simultaneamente interactuando con ambos 
jets en el caso de 3C 273. 



Tabla 6.2: Parametetros adoptados para Cen A y 3C 273. 





Cen A 


3C 273 


Distancia [Mpc] 


3.7 


6.7 x 10 2 


Masa del SMBH [M ] 


6 x 10 7 


7 x 10 9 


Angulo de inclination del jet [°] 


> 50 


~ 15 


Luminosidad del jet [erg s _1 ] 


10 44 


4 x 10 47 


Luminosidad del disco [erg s _1 ] 


5 x 10 41 


2 x 10 46 


Densidad de energfa de los fotones del disco [eV] 


~ 5 x 10 3 


54 


Luminosidad de la BLR [erg s _1 ] 




4 x 10 45 
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6.6.2. Galaxias FR II: 3C 273 

A la distancia d = 6.7 x 10 2 Mpc, 3C 273 es el AGN activo en radio mas potente 
(Courvoisier 1998) con un SMBH de masa M smbh ~ 7 x 10 9 M Q (Paltani & Tiirler 2005). El 
angulo entre la direction del jet con la lmea de la visual es chico, ~ 6°, lo que implica que 
3C 273 es un blazar (Jolley et al. 2009). El espectro completo de esta fuente es variable 
(e.g. Pian et al. 1999) con perfodos que van desde anos (en radio) hasta unas pocas horas 
(en rayos gamma). En HE, 3C 273 fue el primer blazar detectedo en la banda MeV por 
el satelite COS-B, y luego por EGRET en la banda 0.1 - 10 GeV (Hartman et al. 1999). 
Esta fuente fue tambien detectada en energfas GeV por Fermi y AGILE, pero ami no ha 
sido detectada en TeV. Dada la luminosidad del jet de 3C 273, L } 4 x 10 47 erg s _1 
(Kataoka et al. 2002), z- m t resulta »3x 10 17 cm. La luminosidad de la BLR de esta fuente 
es L M r ~ 4 x 10 45 erg s" 1 (Cao & Jiang 1999), y su tamano Ruv ~ 7 x 10 17 cm (Ghissellini 
et al. 2010), lo que implica que las interacciones jet-nube pueden ocurrir por ser z int < Rur- 
La luminosidad del disco es alta, L d 2 x 10 46 erg s -1 , con energfas tipicas de los fotones 
de ^ 54 eV (Grandi & Palumbo 2004). 

La SED de la radiacion no termica generada por las interacciones jet-nube en 3C 273 
es mostrada en la Figura 16.101 (derecha). A la altura z- mt , los procesos radiativos mas 
importantes son el sincrotron y el SSC. Las luminosidades bolometricas generadas por estos 
procesos en una interaccion en z- mt son L sin ~ 6 x 10 38 erg s _1 y Lssc ~ 2 x 10 39 erg s _1 . 
Dados los fuertes campos de radiacion del disco y de la BLR, la emision por encima de 
~ 10 GeV es absorbida por interacciones foton-foton, y el maximo de la emision ocurre 
alrededor de 0.1-1 GeV. Siendo que el numero estimado de nubes en la BLR de 3C 273 es 
~ 10 8 (Dietrich et al. 1999), el factor de llenado que se deduce es / ~ 3 x 10~ 7 y el numero 
de nubes en ambos jets resulta ~ 2 x 10 3 y 5 x 10 5 , para los valores mmimo y maximo dados 
en la Section 16.51 En el caso mas optimista, la luminosidad SSC podria alcanzar valores 
~ 2 x 10 44 erg s _1 . Este valor esta por debajo de la luminosidad detectada por Fermi en 
HE, ~ 3 x 10 46 erg s _1 en el estado estacionario y ~ 1.7 x 10 47 erg s _1 en flares (Soldi et al. 
2009). Sin embargo, la emision detectada probablemente este muy amplificada por efecto 
Doppler y esto tapa la emision que no sufre tal amplificacion. Sin embargo, para AGN 
que no son blazares (galaxias FR II, como Cygnus A), cuya emision no sufre amplificacion 
Doppler, la radiacion en GeV producida por interacciones jet-imbe podrfan ser detectadas. 
En este caso, dado que muchas nubes pueden interactuar simultaneamente con el jet, la 
emision sera estacionaria. 

6.7. Discusion 

En este capitulo hemos estudiado la interaccion de nubes de la BLR con la base de 
los jets en los AGN. Considerando valores razonables para los parametros de las nubes y 
de los jets, estimamos las escalas de tiempo de los procesos dinamicos mas relevantes en 
el escenario analizado, concluyendo que las nubes de la BLR pueden entrar en el jet solo 
a Zj > Zjnt- Para alturas menores que este valor, el jet es muy compacto y las presiones 
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magnetica y cinetica de este destruiran a la nube antes de que esta entre completamente 
en el jet. Cuando la nube interactua significativamente con el jet, choques fuertes son 
generados con la subsecuente emision en rayos gamma de las partfculas aceleradas en tales 
choques. 

Valores del campo magnetico en la region chocada del jet (B hs ) menores que el co- 
rrespondiente a la equiparticion con la energfa de las partfculas no termicas permiten una 
emision significativa de rayos gamma. Para valores mas altos de B\, s (esto es, jets domina- 
dos por el flujo de Poynting), el tratamiento desarrollado aquf no es valido. En tal caso, 
z int puede aun ser definido si en vez de considerar la presion cinetica del jet, consideramos 
la magnetica. Si una nube entra al jet, la aceleracion de partfculas en el bow shock puede 
darse por reconecion magnetica, por ejemplo. El estudio de este caso requiere de un trata- 
miento completamente diferente al expuesto en este capftulo. En general, para valores de 
Bbs mayores que el de equiparticion la emision de rayos gamma por IC es suprimida en 
favor de la radiacion sincrotron al menos que la disipacion magnetica reduzca la intensidad 
del campo magnetico lo sufiente como para que el canal de perdidas IC sea el dominante. 

En fuentes cercanas, como por ejemplo Cen A, la interaction de nubes grandes con los 
jets se deberfa detectar como un evento esporadico, aunque el mimero de estas nubes mas 
grandes y el duty cyclt^ de los flares son diffciles de estimar. Dado que los campos de fotones 
externos son debiles en estas galaxias, fotones de VHE pueden escapar del emisor sin ser 
absorbidos por interacciones con fotones ambientales. Las interacciones jet-nube en galaxias 
FR I cercanas serfan detectables tanto en HE como en VHE como emision esporadica con 
escalas de tiempo de aproximadamente un dfa. La detection de esta emision nos proveerfa 
de information sobre las condiciones del medio de la region nuclear de los AGN (oscurecida 
en general en otras frecuencias por la absorcion del polvo), como asf tambien de la base de 
los jets. 

En fuentes tipo FR II, muchas nubes pueden interactuar simultaneamente con el jet. 
El mimero de estas -N^- depende fuertemente del tiempo de vida de las nubes dentro del 
jet, que es del orden de algunas veces t cn . Sin embargo, notamos que aun despues de la 
fragmentation de las nubes, los bow shocks asociados a cada fragmento pueden todavfa 
formarse y acelerar partfculas eficientemente si los fragmentos se mueven lentamente, a 
una velocidad <C v- y Debido a que las galaxias FR II tienen tasas de acrecion altas, la 
radiacion por encima de 1 GeV producida cerca de la base del jet puede ser sustancialmente 
atenuada por los densos campos de fotones externos (del disco y de la BLR) aunque los 
rayos gamma de energfas < 1 GeV no deberfan ser afectados significativamente. Ya que la 
emision producida en las interacciones jet-nube es isotropica, esta puede ser tapada por la 
emision producida en los jets de los blazares y amplificada por efecto Doppler. No obstante, 
los jets de las galaxias FR II cercanas no muestran un incremento muy importante del flujo 
y estas fuentes pueden emitir rayos gammas por interacciones jet-nube. En el contexto 
de la unification de los AGN (Urry & Padovani 1995), el mimero de AGN que no son 
blazares (intensos en radio) debe ser mucho mayor al de blazares para el mismo valor de 



1 Se define el duty cycle de una fuente como la fraction del tiempo de observation durante la cual se 
producen variaciones significativas de flujo. 
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L y Como se muestra en la Figura [6^9| fuentes cercanas e intensas podrian ser detectables 
por observaciones suficientemente profundas del satelite Fermi. Luego de unos pocos anos 
de exposition, sefiales significativas de estas fuentes podn'an ser encontradas, aportando 
esto una evidencia de que los jets estan dominados por la materia a las alturas de las 
interacciones con las nubes, ademas de informacion sobre las caracterfsticas de la BLR. 



Capitulo 7 
Ciimulos de galaxias 



7.1. Introduccion 

Las galaxias no estan aisladas en el Universo sino que se agrupan de a cientos en 
regiones de un tamano ~ 1 — 10 Mpc, donde forman parte de la estructura filamentosa 
del Universo. Para explicar la formation de estas estructuras a gran escala, el modelo 
actualmente aceptado contempla la existencia de materia oscura frfa (CDM, por Cold Dark 
Matter). Bajo este paradigma, las galaxias se forman por la fusion de estructuras menores 
que dan origen a sistemas de como asf tambien por la acrecion de material 

en forma de gas. Durante el proceso de acrecion y fusion, la materia se acumula en una red 
cosmica con sobredensidades moderadas en forma de filamentos, en cuyas intersecciones se 
ubican las galaxias y los ciimulos de galaxias. 

El proceso de acrecion contempla la caida de gas difuso proveniente del espacio circun- 
dante al cumulo. Este mecanismo es inestable y anisotropico, ya que ocurre a traves de 
flujos de material que se propagan a lo largo de los filamentos que convergen en el cumu- 
lo produciendo fuertes choques de acrecion al impactar con el medio intra cumulo (ICM, 
por Intra Cluster Medium). Por otro lado, la fusion de dos estructuras virializadas puede 
contemplar la colision entre dos galaxias, como se cree que ha ocurrido en Cen A, como 
asf tambien la fusion de dos ciimulos de galaxias. Tanto los procesos de acrecion como los de 
fusion son mucho mas intensos en los ciimulos de galaxias que en las galaxias individuales. 
Esto hace que se produzcan fuertes choques (de acrecion y de fusion), los cuales calientan 
el ICM y aceleran particulas cargadas hasta velocidades relativistas. Tanto la emision del 
gas caliente como de las particulas no termicas caracterizan a los ciimulos de galaxias en 
diferentes bandas del espectro electromagnetico. 

En rayos X, los ciimulos de galaxias son los objetos mas brillantes del cielo alcanzando 
una luminosidad Lx ~ 10 43 — 10 45 erg s -1 . Esta emision es termica, producida por el 
calentamiento del gas intra ciimulo que alcanza temperaturas ~ 10 s K. Sin embargo, la 
detection de estas estructuras en frecuencias radio como asf tambien en el extremo UV 
(0.07 — 0.4 keV) y en los rayos X duros (20 — 80 keV) indica la presencia de actividad 
no termica en estos sitemas. La emision en frecuencias radio detectada en los ciimulos de 
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Figura 7.1: Los cumulos de galaxias Abell 2163 (izquierda) , RX J1347. 5-1145 (medio) y 
ZwCl 2341.1 (derecha) observados en rayos X (amarillo-anaranjado-marron) y en frecuen- 
cias radio (contornos de lmeas llenas) son mostrados. Uno de los halos mas luminosos que 
se han detectado es aquel asociado al cumulo A2163. En el centro de la fuente RX J1347.5- 
1145, la mas luminosa en rayos X, se ha detectado un halo pequeno, mientras que rema- 
nentes dobles han sido detectados en los bordes del cumulo ZwCl 2341.1 (Ferrari 2010). 



galaxias es originada en la interaccion de electrones relativistas con el campo magnetico 
del cumulo (~ AiG). Por otro lado, la radiacion en el extremo UV debe ser producida o bien 
por la componente termica mas frfa o bien por interacciones IC con los fotones del fondo 
cosmico de radiacion (CMB, por Cosmic Microwave Background) con la misma poblacion 
de electrones que produce la emision sincrotron en radio. La radiacion sincrotron tambien 
podria explicar la emision de los rayos X duros, aunque esto es todavia tema de debate. 

Respecto de la emision en radio, aunque su naturaleza no termica esta bien establecida 
por las observaciones, esta puede ser debida a dos tipos de fuentes diferentes: halos o 
remanentes {relics). Los primeros estan localizados en los centros de los cumulos y la 
emision en radio de estos no esta polarizada y se superpone a la emision termica en rayos X. 
Los ultimos se encuentran en la periferia de los cumulos y son regiones de emision polarizada 
y con una estrucura mas irregular. El modelo actualmente aceptado para la formacion 
de estos halos y remanentes es que se producen por choques de acrecion y de fusion, 
respectivamente. En la Figura I7TT1 se muestran estos dos tipos de fuentes. 

Los remanentes han sido observados en muchos cumulos, aunque los mas poderosos 
y extendidos han sido detectados en fuentes con halos centrales, como el cumulo Coma 
(Giovannini et al. 1991) y las fuentes Abell 2163 (Feretti et al. 2001), Abell 2255 (Feretti 
et al. 1997), Abell 2256 (Rottgering et al. 1994) y Abell 2744 (Govoni et al. 2001). Solo 
muy pocos cumulos presentan dos remanentes (opuestos el uno al otro), siendo los ejemplos 
mas prominentes aquellos encontrados en los cumulos Abell 3667 (Rottgering et al. 1997) 
y Abell 3376 (Bagchi et al. 2006). 

Dado que hay una poblacion de particulas relativistas que emiten la radiacion sincrotron 
y quizas la emision de rayos X duros por IC, estas mismas particulas pueden producir 
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tambien emision de mas alta energia. La radiation gamma serfa generada por interacciones 
pp con el ICM si protones son acelerados (ver Volk et al. 1996), o por dispersiones IC entre 
los fotones del CMB y los electrones relativistas (Atoyan & Volk, 2000). Las simulaciones 
cosmologicas hidrodinamicas son utiles para predecir los niveles de emision de rayos gamma 
en objetos virializados a gran escala (Keshet et al. 2003; Pfrommer et al. 2007, 2008). Si 
embargo, los resultados de estos estudios no predicen una contribution significativa al 
fondo difuso de rayos gamma (Berrington & Dermer 2003). Estos resultados teoricos son 
consistentes con el hecho de que no se hayan detectado cumulos de galaxias en HE y 
VHE con los telescopios de rayos gamma que actualmente entan funcionando. Solo han 
sido reportadas evidencias marginales de emision de altas energias del cumulo Abell 1758 
dentro de la caja de error de la fuente 3 EG J1337+5029 (Fegan et al. 2005). Los cumulos 
de galaxias Perseus y Abell 2029 fueron observados por Perkins y colaboradores (2006) 
con el telescopio Cherenkov Whipple, aunque no encontraron evidencias de fuentes de 
rayos gamma puntuales ni extendidas en el rango de los TeV. Por otro lado, observaciones 
recientes de los cumulos Coma y Abell 496 llavadas a cabo con HESS (Domainko et al 
2007) no detectaron emision significativa en rayos gamma a partir de observaciones con 
tiempos de exposition moderados, de ~ 10 — 20 horas. 

Todas las estimaciones hechas al momento de la posible emision de rayos gamma de 
los cumulos de galaxias han sido realizadas para fuentes que no son reales sino simuladas 
numericamente. Por esto, en el presente capftulo estudiaremos los procesos no termicos que 
tienen lugar en los bordes del cumulo Abell 3376, en donde se han detectado en frecuencias 
radio dos remanentes gigantes y simetricos. En lo que sigue, analizaremos tanto el contenido 
de particulas relativistas en los remanentes de este cumulo como asi tambien la emision no 
termica producida por las particulas aceleradas, concentrandonos en el analisis del espectro 
de rayos gamma. 

7.2. El cumulo Abell 3376 

El cumulo de galaxias Abell 3376 esta ubicado a una distancia d ~ 187 Mpc, con un 
corrimiento al rojo (redshift) z ~ 0.046 si consideramos una cosmologia estandar de ma- 
teria oscura frfa ACDM (fi m =0.3, fi A =0.7, H G = 100 h^ 1 km s" 1 Mpc" 1 , con h = 0.7). La 
masa del cumulo ha sido estimada aplicando el teorema del virial a las galaxias miembro 
del mismo. Asumiendo que la distribution de masa de la fuente se comporta de la misma 
manera que la distribution de masa de las galaxias miembro (Girardi et al. 1998), obtene- 
mos una masa virial MA3376 ~ 3.64 x 10 14 hr x M . El correspondiente radio virial resulta 
#A3376 ~ I Ah- 1 Mpc. 

El cumulo Abell 3376 ha sido detectado por los satelites ROSAT y VMM-Newton en 
rayos X siendo esto una evidencia de la actividad de fusion de subcumulos. La luminosidad 
en X de esta fuente es L x (0A — 2.4 keV) ~ 2.48 x 10 44 erg s _1 , producida termicamente 
por el ICM caliente, que tiene una temperatura promedio T x ~ 5.8 x 10 7 K (?s 5 keV). 

Observaciones en radio llevadas a cabo por Bagchi y colaboradores (2006) con el arreglo 
de telescopios VLA muestran la presencia de dos estructuras gigantes en el borde del 
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Figura 7.2: Composition de la emision en frecuencias radio (contornos amarillos) y en 
rayos X (estructura central) del cumulo Abell 3376. Con amarillo se dibuja una elipse 
dentro de la cual estan la emision de rayos X (en el centro) y los remanentes (en el borde) 
(Bagchi et al. 2006). 
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Figura 7.3: Mapas compuestos en frecuencias radio y opticas de los remanentes este (iz- 
quierda) y oeste (derecha) del cumulo Abell 3376. En rojo se muestra la emision en radio 
y en azul la radiation detectada en optico (Bagchi et al. 2006). 
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cumulo, a una distancia ~ 1 Mpc del centre El flujo detectado en radio a v = 1.4 GHz es 
S u = 302 mJy, que corresponde a una luminosidad de 2.1 x 10 40 ergs _1 . Estas estructuras 
tienen las caracterfsticas tfpicas de los radio remanentes (relics) (Giovannini & Ferreti 
2004). Ellas se ajustan bien con una elipse proyectada en el piano del cielo de semiejes 
menor y mayor ~ 1.6 y ~ 2 Mpc, respectivamente, como se muestra en la Figura 17.21 
Adoptando una profundidad de la lmea de la visual en la elipse de ~ 270 kpc (Bagchi et al. 
2006), el volumen del elipsoide tridimensional resulta ~ 0.45 Mpc 3 . De la figura mencionada 
anteriormente inferimos que solo el 20 % de este volumen corresponde a los remanentes en 
radio, es decir, V r ~ 0.09 Mpc 3 . En la Figura [7731 se muestra con mas detalle la estructura 
y tamano de cada uno de los remanentes. 

Algunos parametros del cumulo que son relevantes para los propositos de nuestro es- 
tudio, como la densidad del gas en los remanentes, n r , y la velocidad del choque, t> c h, no 
son provistos por las observaciones y para hallarlos nos valemos de simulaciones numeri- 
cas. Considerando una simulacion hidrodinamica de iV-cuerpos/SPH (Smoothed Particle 
Hydrodynamics) para un cumulo cuya masa virial es M vir ~ 1.3 x 10 15 h^ 1 M (~ MA3376), 
del analisis de los resultados de esta simulacion obtenemos que n r = 2 x 10~ 5 cm -3 y 
f c h = 1000 km s -1 . Teniendo en cuenta que la temperatura del ICM, donde el choque se 
esta propagando, es T ICM ~ 0.1T r (Hoeft et al. 2004) el numero de Mach resulta M ~ 4.2 
(Gabici & Blasi 2003), donde T r es la temperatura del ICM en la region de los remanentes. 
En el resto del caprtulo, expresaremos los resultados numericos adoptando h = 0.7. En 
la Tabla 17.11 listamos los valores de los parametros (observados y estimados) de la fuente 
Abell 3376. 



Parametro valor 



Corrimiento al rojo 


z = 


: 0.046 


Distancia 


d = 


: 187 Mpc 


Masa 


M A3 376 ~ 5.2 x 1O 14 M 


Tamano 


#A3376 ~ 2 MpC 


Volumen de los remanentes 


V t ~ 0.09 Mpc 3 


Densidad 


n r - 


= 2 x 10~ 5 cm- 3 


Velocidad de los choques 


Vch 


= 1000 km s" 1 


Luminosidad en X 


Lx 


= 2.48 x 10 44 erg s" 1 


Luminosidad en radio 


2.1 


x 10 40 ergs _1 


Temperatura del gas 


T x 


« 5.8 x 10 7 K 



Tabla 7.1: Parametros del cumulo de galaxias Abell 3376. Los valores estan dados consi- 
derando h = 0.7. 
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Figura 7.4: Mapas de densidad (izquierda) y de velocidad (derecha) de una simulation 
de un cumulo de galaxias cuya masa virial es ~ 1.3 x 1O 15 /i _1 M . Del analisis de estas 
figuras podemos obtener una estimation de los valores de la densidad y de la velocidad de 
los choques a una distancia ~ 1 Mpc del centro (donde se localizan los remanentes en la 
fuente Abell 3376). Hallamos n r ~ 2 x 10~ 5 cm~ 3 y v ch ~ 1000 km s _1 . 



7.3. Poblacion de particulas relativistas 

Las ondas de choque que se generan (por acrecion y fusion) durante la formation y 
evolution de los cumulos de galaxias son la fuente principal de la termalizacion del ICM y 
de la aceleracion de particulas (Pfrommer et al. 2006). La actividad de las radiogalaxias 
inmersas en los cumulos tambien contribuye a la poblacion de particulas relativistas, in- 
yectando un plasma fosil que se detecta como cavidades en los mapas de brillo superficial 
de rayos X (Churazov et al. 2000). En el cumulo Abell 3376 se han detectado dos radio- 
galaxias. Una de ellas es la fuente MRC 0600-399, asociada con el segundo miembro mas 
luminoso del cumulo, y la otra radiofuente es originada posiblemente por una galaxia elipti- 
ca. Estas dos radiogalaxias se encuentran dentro de la region central del cumulo, donde la 
emision termica es detectada. Por esto, consideramos que las estructuras que se detectan 
en frecuencias radio en el borde del cumulo no estan asociadas a estas radiogalaxias (que 
son fuentes puntuales). 

La morfologia de la emision en frecuencias radio y en rayos X observada en el cumulo 
Abell 3376 sugiere que este esta sufriendo una fusion con otro cumulo de masa similar. 
Como muestran los calculos numericos que simulan la fusion de dos cumulos de galaxias 
(Hoeft et al. 2004), las ondas de choque se propagan en ambas direcciones a lo largo de 
la lfnea que conecta los centros de ambos cumulos con los remanentes casi exclusivamente 
observados en la localization de los frentes de choque. En los casos en los cuales solo se 
detecta un unico remanente, probablemente la fusion haya ocurrido entre dos cumulos de 
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Figura 7.5: Esquema de la fusion de dos cumulos {cores). En la figura de arriba se muestra 
el sistema antes de la fusion y en la de abajo despues. Como puede observarse en el ultimo 
grafico, luego de la colision queda un micleo compuesto por la fusion de ambos cumulos y 
dos choques que se propagan en direcciones opuestas y hacia afuera del sistema (Markevitch 
et al. 2000). 



masas muy diferentes y solo el choque que se propaga a traves del cumulo mas masivo 
perdura en el tiempo y se detecta. El tiempo de vida de estas estructuras suele estimarse 
como r vic i a ~ Rr/vch, donde R r es la distancia recorrida por los choques. En el caso de la 
fuente que estamos estudiando, R x ~ 1 Mpc y t> ch ~ 1000 km s _1 , con lo cual r vida resulta 
~ 1 Gyr. Notamos que estos choques no llegan a un estado estacionario ya que son el 
producto de un evento puntual en el tiempo: una fusion. La energia que llevan los choques 
de fusion se disipa a medida que estos se propagan por el ICM. 

Teniendo en cuenta las evidencias observacionales y los resultados numericos, supone- 
mos que el contenido de particulas relativistas en los remanentes es debido a la aceleracion 
en los choques de fusion, despreciando la posible contribution de las radiogalaxias y de 
los choques de acrecion. En las siguientes subsecciones describiremos los procesos de ace- 
leracion y enfriamiento tanto de electrones (de los cuales tenemos evidencia observacional) 
como de protones, que determinan las distribuciones de energia de estas particulas y su 
subsiguiente evolution. 

7.3.1. Aceleracion de particulas y perdidas radiativas 

La emision difusa en radio producida en la localization de los remanentes es evidencia 
de la presencia de electrones relativistas en esa region. Por otro lado, aunque no tengamos 
evidencia observacional de la presencia de protones en los remanentes, estas particulas 
pueden ser aceleradas tambien de la misma manera que los electrones. Para estimar el 
tiempo en el cual estas particulas se aceleran hasta una energia dada necesitamos conocer 
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Tabla 7.2: Campo magnetico, B v , y constantes de normalizacion de la distribution de 
electrones primarios, k' ei , y protones, k' p , acelerados en los remanentes para los tres casos 
considerados y caracterizados por el parametro a. 



a 


B x 
[G] 


[erg Pe i _1 cm -3 ] 


k p 

[erg Pp_1 cm~ 3 ] 





9 x 10~ y 


1.4 x 10~ i5 




1 


1.1 x 10~ 6 


1.1 x 10~ 15 


8.7 x 10~ 16 


100 


3.4 x 10~ 6 


3 x 10~ 16 


2.5 x 10~ 14 



la magnitud del campo magnetico B r en los remanentes/acelerador. 

Para estimar B r suponemos que la energia se reparte de la misma manera tanto en el 
campo magnetico, Ub, como en las particulas relativistas, u Qt : 

B 2 

=u nt = u ei +u p + u e2 , (7.1) 

donde u ei y u p son las densidades de energia de los electrones y protones primarios, respec- 
tivamente, y u e2 la densidad de energia de los pares secundarios. De la misma manera que 
en el Capitulo HI para estimar u-q a partir de ( 17.11) necesitamos conocer las distribuciones 
de particulas, es decir, las constantes de normalizacion k[, para i = ei,p,e2, y los indices 
p\ de las potencias (rii(Ei) = k' i E i Pl ). Para esto procedemos de la misma manera que se 
ha expuesto en la Seccion 14.3. 2\ fijando u v = au ei y considerando tres valores para el 
parametro libre a: 0, 1 y 100; ademas de tener en cuenta que u e , 2 = f e ±u p . 

Para que la produccion de pares secundarios sea eficiente, es necesario un medio con una 
densidad alta. Sin embargo, esta condicion no se satisface en las regiones perifericas de los 
cumulos, las cuales estan caracterizadas por densidades muy bajas del gas (~ 10 _5 cm~ 3 ). 
Por esto, la contribucion de u e2 en la ecuacion fl7.ll) es despreciable y no la tenemos en 
cuenta en los calculos de B T . Luego, aplicando la restriccion de que el flujo observado en 
radio de la fuente Abell 3376 es S(u = 1.4 GHz) = 302 mJy, podemos estimar el valor de B T 
y de las constantes k\ asumiendo que p 1 = p' p = 2.1. Los resultados obtenidos se muestran 
en la Tabla [7T2l Considerando estos valores obtenemos que p e2 = 2.08 y el cociente entre 
las densidades de energia de los pares e ± y de los protones es muy chico, u e2 /u p ~ 10~ 5 , 
lo que justifica nuestra suposicion previa sobre la contribucion de los pares secundarios al 
flujo detectado en frecuencias radio. 

Asumiendo que la velocidad del choque es = 1000 km s -1 , el tiempo de aceleracion 
result a 

Las perdidas por escape no son relevantes ya que las escalas espaciales que deben recorrer las 
particulas para escapar de los remanentes son muy grandes. Siendo el ancho del remanente 
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~ 2l T ~ 0.6 Mpc (ver la Figura I7T3]) tenemos que el tiempo de conveccion a una velocidad 
~ f c h/4 es muy largo, 

tconv ~ ~ 7 x 10 16 s. (7.3) 

fch/4 

Para estimar el tiempo de difusion consideramos que la escala espacial minima que deben 
recorrer las particulas para difundir y escapar del remanente es l T ~ 0.3 Mpc y hallamos 
que 

tdif ~ 2 x 10 22 ( p- UA V 1 8 . (7.4) 



Perdidas leptonicas 



1 /iG / \ler 



Las perdidas radiativas mas importantes que sufriran los electrones en los remanentes 
de fusiones son por radiacion sincrotron, por Bremsstrahlung relativista y por dispersiones 
IC. Las primeras estan caracterizadas por un tiempo de enfriamiento (ver ( 13. 15ft ) 

f -~ 4 - ixi ° i4 (3" 2 fe)" ls - (7 - 5) 

Para las perdidas por IC consideramos el campo de fotones provistos por el CMB como 
asf tambien los rayos X emitidos termicamente por el ICM caliente. En el primer caso, la 
densidad de energia del CMB es uqmb = 1-2 x 10 _13 ergcm -3 al redshift z correspondiente 
al cumulo Abell 3376. En el ultimo caso ux = L x (4 7r R\ c)~ x (1 + z) A , donde Rx es el radio 
de la region emisora de los fotones X y es ~ 0.15 — 0.3-R v ; r (Balestra 2007), con lo cual en 
la fuente que estamos estudiando resulta Rx ~ 0.3 Mpc, dando «x ~ 1-3 x 10 -15 erg cm -3 . 
Debido a que ux <C mcmb solo consideraremos los fotones del CMB para estimar las 
perdidas por IC. Teniendo en cuenta que la energia caracteristica de los fotones del CMB 
es -Eph^ 5 ~ 1.9 x 10~ 4 eV, solo aquellos electrones con energfas mayores que 1.6 x 10 15 eV 
van a interactuar en el regimen de KN. El tiempo de enfriamiento por IC con los fotones 
del CMB resulta 

ilc = 15 .25 }) +S3v \, {l + lZy2) a . , (7.6) 
ln(l + 0.2 y) (1 + 0.5 y + 1.3 y 2 ) V ' 

donde y = E^EJ (5.1 x 10 5 eV) 2 . 

Finalmente, considerando que la densidad del medio en la localizacion de los remanentes 
es n T = 2 x 10 _5 cm -3 , el tiempo de enfriamiento por Bremsstrahlung relativista es muy 
grande, tBrem ~ 2 x 10 22 s, con lo cual este proceso no sera relevante. En la Figura 17.61 se 
muestran los tiempos de enfriamiento calculados anteriormente juntamente con los tiempos 
de aceleracion y de escape (conveccion) para el caso a — 1. Como puede observarse, la 
energia maxima de los electrones primarios queda determinada por la igualdad t ac = tic Y 
da un valor E™ x ~9x 10 13 eV para diferentes valores de a, como se muetra en la segunda 
columna de la Tabla [7731 Para los electrones con energfas ~ -E™ ax , encontramos que el radio 
de giro de estas particulas es r g ~ 0.1 pc (<C l T ), con lo cual estan contenidas dentro de la 
region de aceleracion. 
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Figura 7.6: Tiempos de aceleracion y de enfriamiento para los electrones primarios. Esta 
figura corresponde al caso con a = 1. Al resto de los casos les corresponden graficos 
similares, ya que la linica diferencia esta en el valor de B r (que afecta a t ac , t sin y tdif) y 
este no cambia demasiado (B T ~ G para los 3 valores de a considerados) . 



Para estimar el estado evolutivo de la poblacion de electrones relativistas acelerados e 
inyectados en los remanentes, consideramos la ecuacion (12. 64ft con un tiempo de vida de la 
fuente r v ida ~ 1 Gyr. Los resultados se muestran en la Figura [7T71 para los casos a = 1 y 100 
y para diferentes tiempos de inyeccion t. Como se muestra en estos graficos, los espectros 
presentan un quiebre en E q ~ 5 y 1 GeV, para los casos a = 1 y 100, respectivamente. Este 
quiebre se debe a que para E ei < E q , las perdidas de energfa por escape son mas eficientes 
que el enfriamiento radiativo. Finalmente, dado que t conv ~ r vi d a , podemos considerar que 
el espectro de los electrones relativistas inyectados en los remanentes alcanza el estado 
presente despues de un tiempo ~ r vida . 



Tabla 7.3: Las energfas maximas obtenidas para los electrones primarios, E™ ax } y protones, 
E'™ ax , acelerados en los remanentes de los choques de fusion del cumulo Abell 3376. Las 
energias maximas E™ ax alcanzadas por los pares e son tambien presentadas. 



a 


[eV] 


[eV] 


[eV] 





9 x 10 13 






1 


9.3 x 10 13 


5.0 x 10 17 


4 x 10 16 


100 


8.8 x 10 13 


1.3 x 10 18 


10 17 
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log(E [eV]) log(E [eV]) 

Figura 7.7: Distribuciones de energia (x_E 2 J de los electro nes primarios para los casos 
a = 1 (derecha) y 100 (izquierda). Se muestran los estados evolutivos a diferentes tiempos 
de vida de los remanentes t r : 0.125 Gyr (lmea punteada), 0.25 Gyr (lmea de puntos y 
rayas), 0.5 Gyr (lmea rayada) y 1 Gyr (lmea llena). El estado presente se alcanza para 
edades ~ r vi d a . 



Perdidas hadronicas 

El tiempo de enfriamiento debido a las colisiones inelasticas pp resulta muy largo, 
t pp ~ 10 20 s, debido a las bajas densidades del medio en los bordes del cumulo. Por esto, la 
produccion de pares no sera muy eficiente en los remanentes. Considerando las perdidas 
por pp la energia maxima que se obtiene es ~ 2 x 10 21 y 5 x 10 21 eV para los casos con 
a = 1 y 100, respectivamente. Sin embargo, para energias E p > 5 x 10 18 eV la produccion 
de pares por la interaction proton-foton del CMB es relevante, y para energias mas alia de 
~ 5 x 10 19 eV las perdidas por produccion de piones son las dominantes (Berezinsky & 
Grigorieva, 1988; Kelner et al. 2008). Sin embargo, no es necesario tener en cuenta estos 
procesos ya que la energia maxima de los protones queda determinada por restricciones 
adicionales: el tamano del acelerador y el tiempo de vida de la fuente. 

Los valores de las energias maximas hallados tanto para los electrones como para los 
protones son validos siempre y cuando estos permanezcan dentro de la region de acele- 
racion cuyo tamano definimos como Z a . Por esto, las particulas deben satisfacer el reque- 
rimiento r g < Z a . Asumiendo que las particulas se aceleran en los choques trazados por 
los remanentes observados en frecuencias radio, adoptamos para el tamano del acelerador 
/ a ~ / r ~ 0.3 Mpc. Debido a que el valor de la energia maxima de los protones obtenido 
en el parrafo anterior a partir de las perdidas radiativas corresponde a r g > Z r , determi- 
nants la energia maxima imponiendo que r g = Z r , con lo cual obtenemos E™ 3 ^ ~ 10 20 y 
2.6 x 10 20 eV, para a = 1 y 100, respectivamente. Sin embargo, el tiempo requerido para 
que los protones alcancen esta energia mediante el mecanismo de aceleracion de Fermi es 
mayor que el tiempo de vida de los remanentes, r vi d a ~ 1 Gyr. Asi, la verdadera energia 
maxima hasta la cual pueden acelerarse los protones en los remanentes se obtiene igualando 
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Figura 7.8: Distribuciones de energia (x E 2 e2 ) de los pares secundarios para los casos a = 1 
(derecha) y 100 (izquierda). Se muestran los estados evolutivos a diferentes tiempos de 
vida de los remanentes t T : 0.125 Gyr (linea punteada), 0.25 Gyr (linea de puntos y rayas), 
0.5 Gyr (linea rayada) y 1 Gyr (linea llena). El tiempo maximo de evolucion es ~ T v i da . 



^acc y TVidai dando _E™ ax ~ 5 x 10 17 y 1.3 x 10 18 eV para a = 1 y 100, respectivamente, 
como se muestra en la Tabla [7731 

La distribucion de energia de los protones acelerados en los remanentes queda deter- 
minada por el tiempo de vida de los liltimos r vida ~ 1 Gyr y debido a que r vida ~ t conv 
podemos considerar que n p (E p ) alcanza el estado presente despues de un tiempo ~ r vida . 
Por esto, 

n p (E p ) ~ q p (E p ) r vida . (7.7) 

Respecto de la produccion de pares e 1 * 1 , si bien la contribucion de estos a la SED 
sera despreciable, la calcularemos por completitud. Para esto necesitamos conocer como es 
la distribucion energetica de estas particulas. La evolucion del espectro de inyeccion de los 
leptones secundarios en los remanentes queda determinada, al igual que en el caso de los 
electrones primarios, por las perdidas radiativas (IC y sincrotron) para E e2 > E q y por el 
tiempo de escape para E e2 < E q , donde E q ~ 5 y 1 GeV, para a = 1 y 100, respectivamente. 
El espectro alcanza el estado de interes, como se muestra en la Figura I7.8[ en un tiempo 

~ Tvida- 



7.4. Distribuciones espectrales de energia 

Conociendo las distribuciones de energia n i: podemos calcular la emision que producen 
las poblaciones de particulas relativistas. Las luminosidades producidas en los remanentes 
estan dadas por la ecuacion 

EphLph = -Ep h jph(-Ep h ) V T , (7.8) 
donde V r ~ 0.1 Mpc 3 es el volumen de la region emisora, es decir, de los remanentes. 
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Figura 7.9: Distribution espectral de energia para el caso puramente leptonico, a = 0. La 
lfnea negra representa la suma de las tres contribuciones (sincrotron, IC y Bremsstrahlung 
relativista) a la SED total. La luminosidad observada con VLA a la frecuencia v = 1.4 GHz 
del ciimulo Abel 3376 se representa con un cuadrado rojo. 



7.4.1. Interacciones leptonicas 

Las emisividades diferenciales j p h(-Eph) producidas por los leptones por radiation sin- 
crotron, dispersiones IC y Bremsstrahlung relativista son calculadas usando las formulas 
dadas en el Capftulo [3j Las luminosidades espetificas son halladas a traves de la ecua- 
cion ( I7.8p . En los casos a = 1 y 100, en los cuales se asume que existe una poblacion de 
protones relativistas, las contribuciones de los pares secundarios a las SEDs tambien se 
estiman por completitud. 

Los resultados de las luminosidades espetificas producidas en el caso donde solo se 
aceleran electrones, caracterizado por el parametro a = 0, se muestran en la Figura 17.91 
Podemos ver que las interacciones IC son el proceso radiativo mas importante, con una 
luminosidad Lie ~ 9.1 x 10 41 ergs _1 , a energias E p ^ > 0.1 MeV y con un cut-off en 
~ 10 TeV. La luminosidad producida por Bremsstrahlung relativista es despreciable, en 
concordancia con los resultados mostrados en la Figura 17.61 

7.4.2. Interacciones hadronicas 

La emisividad correspondiente a las interacciones pp se calcula considerando las ecuacio- 
nes (I3.39P y (13.411) . y luego la luminosidad espeeffica se calcula a traves de la ecuacion (\7.8L 
Los resultados de nuestros calculos se muestran en la Figura 17.101 

En el caso a = 1, de la misma manera que como ocurre en el caso a = 0, la SED 
esta dominada por las interacciones IC, con una luminosidad Lie ~ 7.4 x 10 41 erg s _1 . La 
envision a energias mayores que ~ 1 GeV es producida por el decaimiento de piones neutros, 
alcanzando una luminosidad L pp ~ 1.6 x 10 38 erg s _1 con un cut-off en E p h ~ 10 17 eV. Sin 
embargo, su contribution a la SED es evidente para energias > 10 TeV, como se muestra 
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Figura 7.10: Distribuciones espectrales de energia para los casos con a = 1 (izquierda) y 
100 (derecha), en los cuales consideramos tanto una poblacion leptonica como hadronica 
de particulas primarias y la poblacion de pares e ± producida. La emision producida por 
los pares secundarios es menor que la correspondiente a los electrones primarios. Como 
mencionamos en la Figura 17.91 la linea negra indica la SED total (suma de todas las 
contribuciones, leptonicas y hadronica) y el cuadrado rojo indica la luminosidad observada 
con VLA. 



en la Figura I7TTU1 (izquierda) . 

En el caso correspondiente al valor a = 100, tenemos que u p > u ei , pero esto no implica 
un flujo significativo de fotones producidos en las interacciones pp ya que la densidad de 
particulas bianco es constante y muy baja en la localization de los remanentes. Por esto, el 
espectro de rayos gamma producidos por el decaimiento de los 7r° no domina la SED, como 
puede observarse en la Figura [7TTU1 (derecha). La luminosidad L pp es ~ 4.2 x 10 39 ergs -1 , 
siendo apenas mayor que la correspondiente al caso con a — 1. Por otro lado, contrariamente 
a lo que ocurre en los casos con a — y a — 1, la emision por IC es menor que la 
radiacion sincrotron, con luminosidades Ljq ~ 7.1 x 10 40 ergs _1 y L sin ~ 3.8 x 10 41 ergs~ 1 , 
respectivamente. El hecho de que la radiacion por dispersiones IC sea reducida en el caso 
con a = 100 puede explicarse como una consecuencia de la equiparticion asumida entre 
la densidad de energia magnetica y de las particulas relativistas. Como puede verse en la 
Tabla I7T2"| el campo magnetico se hace mas intenso cuando a crece. Esto reduce la cantidad 
de energia de los electrones que se radia por interacciones IC en comparacion con la que 
se pierde a traves del proceso sincrotron. Ademas, la densidad de energia en electrones 
relativistas tambien se reduce para poder ajustar el espectro con el flujo observado con un 
valor de B r mayor. Esto ultimo hace que Ljq sea menor en terminos absolutos en el caso 
a = 100. 

Respecto de la emision producida por los pares e ± , encontramos que sus contribuciones 
a la SED por los diferentes procesos expuestos en la seccion 17.3.11 son mayores en el caso 
con a = 100 que con a = 1, pero aun son irrelevantes. Esto es consistente con la baja 
eficiencia de las interacciones pp que tienen lugar en la periferia del cumulo. 



134 



Cumulos de galaxias 



7.5. Discusion 

En este capftulo hemos presentado el estudio realizado sobre el cumulo de galaxias 
Abell 3376. Esta fuente presenta fuertes evidencias de una actividad de fusion, caracterizada 
por grandes estructuras detectadas en frecuencias radio en la periferia del cumulo y que han 
sido asociadas con remanentes (Bagchi et al. 2006). Hemos modelado la emision producida 
en estos remanentes en un amplio rango de frecuencias, desde radio hasta rayos gamma, 
considerando la contribution de diferentes procesos radiativos y diferentes relaciones entre 
las densidades de energia de los protones y los electrones relativistas. 

Las perdidas leptonicas mas importantes que afectan la distribution de parti'culas son 
la radiation sincrotron y las dispersiones IC con fotones del CMB. Los protones relativistas 
se enfrian por interacciones pp con el ICM, produciendo rayos gamma e inyectando una 
poblacion de pares e ± los cuales se enfrian por los mismos procesos que los electrones 
primarios. Los parametros involucrados en estos procesos no termicos, como el campo 
magnetico y las distribuciones de parriculas relativistas, se han estimado a traves del ffujo 
observado en radio de la fuente y de suponer equiparticion de la energia magnetica y no 
termica. Ademas, las observaciones en rayos X nos permiten extraer information sobre el 
volumen de los remanentes. Por otro lado, los resultados de simulaciones numericas nos 
proveen de valores tipicos para la densidad del medio en los bordes del cumulo, como 
asi tambien de la velocidad de los choques de fusion. 

Las SEDs mostradas en las Figuras I7T91 y 17. 101 se han calculado considerando la contri- 
bution de diferentes procesos radiativos no termicos en un amplio rango de energias, desde 
radio hasta rayos gamma. En este capftulo hemos estudiado particularmente la detectabi- 
lidad del cumulo Abell 3376 en HE y VHE. Encontramos que los procesos radiativos mas 
importantes que contribuyen a la emision de rayos gamma son las interacciones IC y el 
decaimiento de 7r°, siendo la luminosidad de este ultimo mucho menor que la alcanzada 
por el primero, aiin en el caso en el cual u p = 100 u ei . La luminosidad por pp obtenida 
en este caso es L pp ~ 4 x 10 39 erg s _1 . La luminosidad por dispersiones IC de fotones del 
CMB es Lie ~ 7.4 x 10 41 erg s" 1 . Asi, los rayos gamma emitidos en el cumulo Abell 3376, 
y bajo las suposiciones de nuestro modelo, podrian ser detectables con Fermi y HESS, que 
operan en las bandas (aproximadas) 100 MeV-100 GeV y 0.1-10 TeV, respectivamente. El 
futuro telescopio CTA podria permitir incluso la detection de la componente hadronica a 
mas alt as energias. 

Finalmente, la proximidad del cumulo Abell 3376 y su posible contenido alto de parti'cu- 
las relativistas hacen esta fuente interesante para ser observada tanto con instrumentos que 
actualmente estan funcionando como asi tambien con futuros telescopios Cherenckov ya 
planeados. Notamos que hasta el di'a de hoy no se han detectado cumulos de galaxias en ra- 
yos gamma y que la fuente Abell 3376 es uno de los mejores candidatos para ser observado 
en esta banda de energia. 



Capitulo 8 

Conclusiones y perspectivas 



A lo largo de esta tesis hemos estudiado los procesos no termicos que tienen lugar en 
diferentes fuentes astrofisicas, tanto galacticas (YSOs, MQs) como extragalacticas (AGNs, 
ciimulos de galaxias). A continuation resumimos las investigaciones realizadas en cada tipo 
de fuente y luego damos las conclusiones generates de la tesis. 

■ YSOs: Algunas estrellas masivas en formacion presentan jets y cuando estos son 
frenados por el material de la nube molecular en la cual se hallan embebidos, se 
producen choques terminales fuertes que pueden acelerar particulas hasta energias 
relativistas. Estas particulas emiten en frecuencias radio por mecanismo sincrotron 
y es asi como se detectan lobulos no-termicos en las regiones terminales de los jets 
que emanan de protoestrellas de gran masa. 

Hemos modelado uno de los lobulos que forman parte el sistema triple en radio 
asociado a la fuente IRAS 16547-4247, la protoestrella mas luminosa detectada a la 
fecha (L* ~ 6.2 x 10 4 L©). Suponiendo equiparticion entre las densidades de energia 
magnetica y de las particulas relativistas y valiendonos de los flujos observados en 
frecuencias radio hemos estimado el campo magnetico y las distribuciones espectrales 
de energia para diferentes relaciones entre las densidades de energia de los protones y 
electrones acelerados en el lobulo. Luego, calculamos las SEDs considerando diferentes 
procesos radiativos no termicos. 

Los procesos dominantes resultan ser Bremsstrahlung relativista y las interacciones 
pp, ya que la densidad de la nube molecular es muy grande, ~ 10 5 cm -3 . Los niveles 
de emision predichos son ~ 10 32 — 10 33 erg s _1 en altas y muy altas energias, con 
lo cual la fuente seria detectable con el satelite Fermi y con el futuro arreglo de 
telescopios Cherenkov CTA. 

Las protoestrellas de gran masa, de ser detectadas en altas energias, serian un nuevo 
tipo de fuentes de rayos gamma. De esta manera, se abre una nueva ventana del 
espectro electromagnetico a traves de la cual se puede obtener information sobre 
como es el proceso de formacion de las estrellas tempranas. Este proceso aun hoy no 
esta claramente establecido. Especificamente, se podran estudiar con nuevas herra- 
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mientas las propiedades fisicas (como la densidad y campo magnetico) de los jets, 
asi como las del medio cincundante. 

■ MQs: Existen evidencias observacionales de que los vientos de las estrellas de gran 
masa no son homogeneos sino que tienen una estructura porosa. En los HMMQs, 
algunos clumps del viento de la estrella companera pueden llegar hasta los jets ge- 
nerados por el objeto compacto. Debido a la interaccion del clump con el material 
del jet se producen dos choques: uno en el jet y otro en el clump. En el primero se 
pueden acelerar particulas que luego radiaran localmente por diferentes procesos no 
termicos, y el segundo calienta el material del clump, que radia termicamente. Las 
particulas aceleradas en el bow shock se enfrian tanto en el jet como en el clump, ya 
que las mas energeticas se difunden hasta alii. Hemos calculado los procesos dinami- 
cos y radiativos mas importantes, considerando diferentes valores para el tamano de 
los clumps y para el campo magnetico en la region del jet donde ocurre la interaccion 
(~ a mq /2). 

De acuerdo a un estudio de la dinamica de la interaccion jet-clump, hemos podido es- 
timar la velocidad del choque en el clump y predecimos que la nube no sera destruida 
antes de que las particulas aceleradas en el bow shock se enfrien significativamente. 
Considerando diferentes valores para el tamano de los clumps (R c = 10 10 y 10 11 cm) 
y del campo magnetico en la region de interaccion (B^ = 150 y 1 G), hemos cal- 
culado las SEDs, prediciendo emision significativa a lo largo de todo el espectro 
electromagnetico. Las luminosidades mas altas obtenidas en rayos gamma son las 
producidas por interacciones IC con los fotones emitidos por la estrella companera, 
Lie ~ 10 35 erg s -1 , en HE, y en VHE se han alcanzado luminosidades tan altas como 
L, pp ~ 10 32 erg s _1 producidas por colisiones pp en el clump. 

De acuerdo al factor de llenado de clumps en el viento que hemos supuesto (/ ~ 
0.005), el numero de clumps que simultaneamente pueden interactuar con el jet es 
~ 350 si -R c = 10 10 cm y ~ 0.5 si R c = 10 11 cm. Luego, la emision sera estacionaria 
en el primer caso, con una luminosidad total ~ 350 veces mas alta que la de una 
interaccion simple, con una escala de tiempo de los flares de ~ 1 hora, y recurrencia 
en el plazo de varias horas. De esta manera, las interacciones jet-clump son una 
posible explication a la production de los flares observados en algunas binarias de 
rayos gamma. 

Los niveles de emision alcanzados son detectables por intrumentos como Fermi en 
HE y los telescopios Cherenkov de nueva generation en VHE. Esta emision, de ser 
detectada, nos proveeria importante information sobre los jets de los HMMQs, como 
asi tambien de los vientos de las estrellas de gran masa (e.g. Cygnus X-l, Cygnus X- 
3). 

■ AGNs: Un estudio similar al realizado sobre los HMMQs ha sido desarrollado en el 
escenario de los AGNs, considerando que nubes de la BLR pueden penetrar en los 
jets que emanan de las cercanias del SMBH. Debido a la interaccion de una nube 
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con uno de los jets se produce un bow shock fuerte en el jet, donde pueden acelerarse 
particulas eficientemente, y otro choque mas debil en la nube. Como en el caso de los 
HMMQs, las particulas aceleradas en el bow shock se enfrfan tanto en el jet como en 
la nube, ya que las mas energeticas se difunden hasta alK. 

Los procesos dinamicos que se desarrollan permiten que la nube entre y permanezca 
entera suficiente tiempo como para que las particulas relativistas aceleradas en el 
bow shock puedan radiar una fraction significativa de su energfa. Las luminosidades 
obtenidas por la interaccion de una sola nube con el jet son relativamente bajas, 
~ 10 38 erg s _1 , pero si muchas nubes se encuentran simultaneamente interactuando 
con ambos jets entonces la contribution total es significativa y detectable en fuentes 
no alineadas (e.g. radiogalaxias) . Debido a que el emisor esta practicamente quieto, la 
radiation no es amplificada por efecto Doppler como acurre en la emision producida 
en los jets de los blazares. La interaccion de nubes de la BLR con jets podrfa ser un 
mecanismo que explique la emision observada en algunos AGNs. 

Finalmente, hemos aplicado nuestro modelo a galaxias FR I y FR II. En las primeras, 
si bien no hay detecciones claras de la presencia de una BLR, es posible suponer una 
poblacion de nubes oscuras aunque no podemos estimar el mimero de ellas. Por 
esto, consideramos aqui solo la interaccion de una nube con el jet. Debido a que 
estas fuentes son cercanas (Cen A, M 87) la interaccion de una sola nube de tamaiio 
R n ~ 10 14 cm podria producir flares detectables. Sin embargo, en las galaxias FR II 
la cantidad de nubes es grande, y la interaccion de muchas de ellas con ambos jets 
del AGN puede producir emision detectable si los jets no estan alineados con la lfnea 
de la visual. 

■ Cumulos de galaxias: Los cumulos de galaxias son candidatos a ser fuentes de rayos 
gamma, ya que hay evidencias de la existencia de particulas relativistas que emi- 
ten en frecuencias radio. En particular, en esta tesis hemos estudiado la production 
de rayos gamma en los remanentes (relics) detectados en el borde de la fuente cer- 
cana Abell 3376. Los remanentes son trazadores de choques de fusion. El cumulo 
Abell 3376 presenta grandes estructuras en frecuencias radio, y que han sido asocia- 
das a remanentes de choques producidos por la fusion entre dos cumulos de galaxias. 
Estos remanentes sugieren un rico contenido de particulas relativistas del ICM en 
estas estructuras. Suponemos que la aceleracion de particulas tiene lugar en los hor- 
des de los cumulos donde se detectan los remanentes radio. Estas particulas luego se 
enfrian por diferentes procesos no termicos. 

Ademas de los electrones primarios y protones acelerados en los choques, considera- 
mos tambien una poblacion de pares secundarios generados por interacciones pp de 
los protones relativistas con el ICM. Los parametros involucrados en los procesos no 
termicos han sido estimados considerando los datos observacionales como asi tam- 
bien suposiciones como la de equiparticion de la energia magnetica y no termica. 
Por otro lado, hemos hecho uso de simulaciones numericas para obtener aquellos 
parametros del medio (velocidad del choque y densidad del ICM) no provistos por 
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las observaciones. 

La emision no termica mas intensa producida en los bordes de este ciimulo es debida 
a las interacciones IC de electrones primarios con los fotones del CMB, alcanzando 
luminosidades ~ 10 41 erg s" 1 en HE y VHE. Si bien al dfa de hoy no se han detectado 
cumulos de galaxias en rayos gamma, la fuente Abell 3376 es un buen candidato para 
ser detectado por medio de exposiciones prolongadas de acuerdo a los resultados de 
nuestro modelo y a su proximidad. 

Las ondas de choque se producen en diferentes tipos de fuentes astrofisicas y por diversos 
mecanismos. Los frentes de choque asociados a estas ondas tienen tamanos que van desde 
~ 10~ 2 pc en los jets de los YSOs hasta ~ 1 Mpc en los cumulos de galaxias. Sin embargo, 
las velocidades de estos choques no son proporcionales a su tamaho, ya que en YSOs y 
en cumulos de galaxias i> c h ~ 1000 km s" 1 mientras que en los jets relativistas de los 
MQs y AGNs t> c h ~ c. En todos estos choques se pueden acelerar particulas cargadas 
hasta velocidades relativistas. El tamano de los frentes de choque influye en las perdidas 
de escape, que se tornan muy lentas si aquellos son muy grandes. El tamano tambien es 
relevante para las densidades de energia de los campos ambientales. Si las densidades de 
energia magnetica y de particulas (fotones y materia) son bajas, entonces las perdidas 
radiativas tampoco son muy eficientes y las energias maximas seran altas, > 1 PeV. Sin 
embargo, si las densidades de energia de los campos ambientales son altas, es posible 
alcanzar energias maximas de ~ 1-10 TeV en el regimen de saturacion. 

La competencia entre los tiempos dinamicos y radiativos determina la eficiencia radia- 
tiva de la fuente. Si t ra( j < tdm en algiin rango de energias de la distribution de particulas 
relativistas, entonces las particulas con esas energias radiaran significativamente y se dice 
que han alcanzado el regimen de saturacion. En todos los escenarios explorados en esta 
tesis, el regimen de saturacion es alcanzado por los electrones en algiin rango de energia y 
para algiin conjunto de parametros considerado. Sin embargo, solo en aquellos casos en los 
cuales las particulas mas energeticas son las que saturan, la emision de la fuente resulta 
intensa y detectable. De ser confirmados los resultados hallados en cada tipo de fuente 
que hemos estudiado, se aportarian nuevos conocimientos en cada campo de investigation 
respectivo. 

En lo que refiere a las estrellas de gran masa, por un lado modelamos la emision de 
la protoestrella asociada a la fuente IRAS 16547-4247 obteniendo flujos detectables en 
rayos gamma, lo cual solo habia sido someramente sugerido anteriormente por Henriksen 
y colaboradores (1991), y siendo asi estos objetos un nuevo tipo de posibles de fuentes 
de rayos gamma. Esto ha sido estadisticamente comprobado a traves de la muy buena 
correlation encontrada entre las fuentes detectadas por Fermi y la localization en el piano 
del cielo de protoestrellas de gran masa (Munar y colaboradores, comunicacion personal). 

Por otro lado, las propiedades de las inhomegeneidades de los vientos de las estrellas 
tempranas pueden conocerse a traves de observaciones en rayos gamma, ya que las in- 
teracciones de estas con los jets de los MQs producirian efectos detectables. Para poder 
determinar la microestructura de la variabilidad es necesario el uso de telescopios con gran 
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resolution temporal (gran sensibilidad) . Un instrumento apropiado para esto sera el arreglo 
de telescopios Cherenkov CTA. 

Las interacciones de obstaculos con jets nos proveen un nuevo mecanismo de radiation, 
ya que ademas de acelerarse las particulas en los bow shocks que se forman en los jets, 
el obstaculo sirve de bianco para las interacciones de estas particulas relativistas. En el 
campo de las binarias de rayos gamma de alta masa, el proceso descripto puede ocurrir si 
los clumps del viento de la estrella companera llegan hasta el jet generado por el agujero 
negro y asi podrian producirse los flares detectados en algunas de estas fuentes. En el 
campo de las galaxias activas, las interacciones de nubes de la BLR con la base de los 
jets podrian explicar la emision esporadica observada en galaxias FR II cercanas y la 
emision estacionaria de los AGNs mas lejanos que no son blazares. Como trabajo a futuro 
en este campo, haremos una aplicacion del modelo expuesto en el Capitulo E] a la fuente 
3C 120, de la cual contamos con datos observacionales de alta resolution espacial. De 
estos datos es posible extraer information sobre la microestructura del jet, que presenta 
inhomogeneidades de ~ cm (Gomez et al. 2008). Como una consecuencia natural de los 
modelos desarrollados de interacciones de jets con obstaculos, nos proponemos estudiar 
que sucede cuando una estrella masiva interactua con el jet de un AGN. Si bien aqui el 
choque en el jet se produciria por la colision entre el viento de la estrella y el material del 
jet, gran parte del modelo ya desarollado puede aplicarse a este nuevo escenario. 

En el topico de los cumulos de galaxias ofrecemos un candidato concreto para ser 
observado con los telescopios de rayos gamma. A traves de una modelizacion de la fuente 
Abell 3376, estimamos que el flujo emitido por esta seria detectable por los instrumentos 
actuales que observan en HE como asi tambien por los futuros telescopios ya planeados. 
Sin embargo, el satelite Fermi no ha detectado ningiin cumulo de galaxias en tres anos de 
observation. Esto muestra que alguno de los valores de los parametros supuestos por los 
modelos actuales esta siendo sobreestimado. Debido a que el efecto mas importante es el 
IC sobre los fotones del CMB, y el valor de u cm b esta bien determinado, es licito pensar que 
quizas la poblacion de particulas relativistas no es tan importante. Esto podria ocurrir si 
los choques de fusion no fuesen tan intensos como resulta de las simulaciones numericas y 
la velocidad de estos fuese menor que 1000 km s -1 . Como trabajo a futuro en este campo, 
nos proponemos hacer un analisis del espacio de valores de los parametros considerados en 
la modelizacion de la emision de rayos gamma de los cumulos de galaxias. 

Finalmente, de ser detectadas las fuentes propuestas en esta tesis y de comprobarse los 
resultados obtenidos al estudiar las interacciones de jets con obstaculos, se increment arian 
las clases de emisores de rayos gamma y se podrian explicar la emisiones esporadicas 
producidas en algunas fuentes, mediante mecanismos no explorados anteriormente. 
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Apendice A 

Lista de acronimos 



Debido al uso generalizado de las siglas inglesas en la jerga astronomica referida al tenia 
de esta tesis, se han utilizado a lo largo de este trabajo una gran cantidad de acronimos 
que se refleren a palabras en la lengua mencionada. 

AGN: Active Galactic Nucleus 

ATCA: Australia Telescope Compact Array 

BH: Black Hole 

BLR: Broad Line Region 

CDM: Cold Dark Matter 

CMB: Cosmic Microwave Background 

CTA: Cherenkov Telescope Array 

EC: External Compton 

FR: Faranoff-Rayleigh 

FSRQ: flat spectrum radio-quasar 

GRB: Gamma-Ray Burst 

HE: High Energy 

HH: Herbig-Haro 

HMMQ: High Mass Microquasar 

IC: Inverse Compton 

ICM: Intra Cluster Medium IR: Infra Red 

KH: Kelvin-Helmholtz 

MHD: Magnetohydrodynamics 

MQ: Microquasar 

RT: Rayleigh-Taylor 

SED: Spectral Energy Distribution 

SMBH: Super-Massive Black Hole 

SPH: Smoothed Particle Hydrodynamics 

SSC: Synchrotron Self Compton 

VHE: Very High Energy 

VLA: Very Large Array 
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Lista de acronimos 



WR: Wolf-Rayet 

YSO: Young Stellar Object 
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